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I.  INTRODUCTION 


A  recurring  problem  in  vulnerability  analysis  is  the  requirement  to  predict  the  damage 
to  structural  panels  produced  by  the  detonation  of  high  explosive  charges  in  close  proxim¬ 
ity  to  the  target  panels.  Examples  of  problems  where  this  requirement  arises  are  (a)  the 
blast  effect  of  shaped  charge  warheads  against  lightly  armored  vehicles  and  (b)  the  effect 
of  small  caliber  high  explosive  shells  which  detonate  close  to  aircraft  panels.  Although 
both  these  examples  involve  additional  damaging  effects  (e.g.,  the  shaped  charge  jet,  and 
shell  casing  fragments), this  report  is  concerned  only  with  a  rational  analysis  of  the 
response  induced  by  the  air  blast,  leaving  the  combined  effects  of  all  lethal  mechanisms 
for  subsequent  consideration. 

While  the  methodology  for  blast  loading  prediction  is  far  from  satisfactory  at  present,  in 
the  sequel  it  will  be  assumed  that  the  blast  pressure  on  the  target  surface  is  a  known  func¬ 
tion  p  (r,  6 ,  t).*  Since  we  will  be  concerned  with  the  blast  from  conventional  chemical 
explosives  the  pressure  will  be  of  high  intensity  and  short  duration,  resulting  in  the 
delivery  of  a  significant  impulse  to  the  target.  However,  the  peak  pressure  will  be  assumed 
insufficient  to  produce  spallation  from  the  back  side  of  the  target  panel,  thus  excluding 
consideration  of  explosives  in  contact  with  the  target.  This  is  not  a  serious  restriction, 
since  stand-off  explosions  can  produce  catastrophic  damage  to  structural  panels. 

It  is  now  necessary  to  adopt  a  structural  response  analysis  tool  which  has  attributes 
sufficient  for  adequate  modeling  of  the  physical  phenomena  expected  to  occur.  Specifically 
the  analysis  should  be  capable  of  treating  finite  amplitude  elastoplastic  response  of  shell 
structures  having  a  variety  of  physical  edge  conditions,  which  are  subjected  to  transient 
distributions  of  surface  pressures.  It  should  have  flexibility  regarding  material  constitutive 
representation,  including  strain  hardening  and  strain-rate  dependence.  It  should  also  be 
possible  to  introduce  various  material  failure  models  into  this  analysis  and,  preferably,  to 
perform  some  post-failure  calculations.  Further,  since  the  blast  load  will  initially  appear  at 
some  interior  point  on  the  panel  and  then  spread  rapidly  to  cover  the  entire  panel,  the 
analysis  tool  should  properly  account  for  the  propagation  of  shear  waves  as  well  as  flexural 
and  membrane  waves.  Finally,  it  is  desired  to  avoid  a  general  three-dimensional  analysis,  if 
possible,  for  reasons  of  computational  economy. 

Taking  the  foregoing  considerations  into  account,  a  variety  of  available  finite  element 
and  finite  difference  computer  programs  were  reviewed  and  it  was  decided  that  the 
PETROS  4  code1,2  provided  the  best  point  of  departure  for  meeting  these  requirements. 


See  Nomenclature,  p.  41  .for  definition  of  symbols. 

*  S.  D.  Pirotin,  L.  Mori  no,  E.  A.  Witmer,  and  J.  W.  Leech,  " Finite-Difference  Analysis  for 
Predicting  Large  Elastic-Plastic  Transient  Deformations  of  Variable-Thickness  Kirchhoff,  Soft 
Borded  Thin,  and  Transverse-Shear  Deformable  Thicker  Shells US  Army  Ballistic  Research 
Laboratory  Contract  Report  No.  315,  September  1976.  AD  B013924L 

2  S.  D  Pirotin,  B.  A.  Berg,  and  E.  A.  Witmer,  ” PETROS  4:  New  Developments  and  Program 
Manual  for  the  Finite-Difference  Calculation  of  Large  Elastic-Plastic,  and/or  Viscoelastic 

Transient  Deformations  of  Multilayer  Variable-Thickness  (7)  Thin  Hard-Bonded , 
(2)  Moderately -Thick  Hard-Bonded,  or  ( SJ  Thin  Soft-Bonded  Shells,"  US  Army 
Ballistic  Research  Laboratory  Contract  Report  No.  376,  September  7976. 


This  report  is  concerned  with  documentation  of  the  modifications  to  this  code  which  were 
required  to  achieve  the  desired  analysis  capability. 


II.  DEFICIENCIES  OBSERVED  IN  THE  ORIGINAL  PETROS  4  PROGRAM 

The  user  of  the  PETROS  4  code  has  first  to  choose  one  of  the  three  versions  referred 
to  in  the  title  of  Reference  2.  For  the  reasons  previously  stated  the  primary  attention  has 
been  given  to  the  moderately-thick  hard-bonded*  transverse  shear  deformable  option 
although  occasional  use  of  the  thin  hard-bonded  KirchhofT  shell  model  has  been  made  for 
comparison  purposes.  Another  choice  to  be  made  is  the  plasticity  theory  to  be  employed, 
which  is  selected  by  the  value  assigned  to  the  input  variable  ISTRES,  as  follows: 


ISTRES 


Plasticity  Theory 

Mechanical  sublayer  model^’*’^,  3-D  stress 

Prandtl-Reuss  model,  3-D  stress 

Mechanical  sublayer  model,  T3j  —  0  (j*=  1,2,3) 

This  option  was  recommended  for  the  KirchhofT  shell. 


Mechanical  sublayer  model,  T  —  0 

(This  stress  component  is  oriented  along  the  normal  to 

the  shell  reference  surface.) 


The  PETROS  4  code  can  treat  shells  composed  of  layers  qf  different  materials  (although  the 
immediate  interest  is  in  applications  involving  only  a  single  layer). 

H.  F.  Bohnenblust,  and  P.  Duwez,  " Some  Properties  of  a  Mechanical  Model  of  Plasticity," 
Journal  of  Applied  Mechanics,  Vol.  15,  No.  3,  September  1948,  pp.  222-225. 

G.  N.  White,  Jr.,  " Application  of  the  Theory  of  Perfectly  Plastic  Solids  to  Stress  Analysis  of 
Strain  Hardening  Solid,"  Graduate  Div.  of  Applied  Math.,  Brown  University  Tech  Report  51, 
August  1950. 

J.  F.  Besseling,  "A  Theory  of  Plastic  Flow  for  Anisotropic  Hardening  in  Plastic  Deformation  of 
an  Initially  Isotropic  Material ,"  Report  5410,  National  Aeronautical  Research  Institute, 
Amsterdam,  The  Netherlands,  1953. 


Another  input  quality  which  must  be  selected  is  INORML,  which  controls  the  manner 
in  which  the  variable  D3  is  calculated.  D3  is  the_oomponent  of  the  vector  D  in  the  direc¬ 
tion  of  the  normal  to  the  reference  surface  and  D  represent  the  three  additional  degrees- 
of-freedom  of  the  SHEAR  model  besides  those  of  the  Kirchhoff  model  at  each  mesh 
point.  The  options  for  INORML  are: 

INORML  D  Calculation 

0  The  cartesian  components  of  D  are  calculated 
using  three  equations  of  motion 

1  is  set  to  zero  after  D  is  calculated 

2  The  incremental  change  in  D3  is  calculated 
from  the  incremental  strain  dy/ 

using  the  elastoplastic  constitutive  relations; 
this  corresponds  to  a  thickness  change  which 
affects  the  stresses  at  the  next  time  step. 


A.  Stress  Calculation  Inconsistencies 

In  order  to  treat  the  problem  of  a  blast-loaded  structural  plate  it  would  appear  appropri¬ 
ate  to  use  the  SHEAR  version  of  the  PETROS  4  code  with  the  options  ISTRES  -  0, 
INORML  —  0  since  these  are  the  most  general  (least  restrictive)  choices  available.  This 
combination  of  options  has  been  employed  to  treat  the  following  physical  example: 

A  square  plate  of  rolled  homogeneous  steel  armor,  0.1905m  (7.50  in)  by  0.1905m 
(7.50  in)  by  9.53mm  (0.375  in)  thickness.  Young’s  modulus:  2.068GPa  (30x10°  psi) 
Poisson’s  ratio:  0.25.  The  uniaxial  strain-hardening  characteristics  of  this  material 
were  represented  in  the  mechanical  sublayer  model  by  the  following  stress-strain  coor¬ 
dinates  (connected  by  linear  segments): 

Coordinate  No.  _ Stress _  Strain 

1  1.048  GPa  (152000  psi)  0.005067 

2  1.145  GPa  (166000  psi)  0.0135 

3  1.248  GPa  (181000  psi)  0.0530 

4  1.675  GPa  (243000  psi)  0.2800 

The  boundary  conditions  imposed  on  all  four  edges  were  complete  fixity.  The  plate 
was  loaded  by  the  blast  from  a  0.907kg  (2  pound)  spherical  pentolite  charge  detonated 
63.5mm  (2.5  in)  above  its  midpoint 


Results  of  PETROS  4  calculations  of  transient  stress  components  using  the  cited  options 
are  shown  in  Figure  1.  The  location  of  these  stresses  is  a  point  9.53mm  (0.375  in)  from 
the  plate  midpoint  along  a  line  perpendicular  to  an  edge  and  0.661mm  (0.026  in)  above 
the  lower  surface.  One  sees  that  all  three  normal  stress  components  have  become  very 
large  (as  has  the  mean  or  hydrostatic  stress  rs)  by  the  end  of  this  short  run  and  it  may  be 
readily  verified  that  the  elastoplastic  stress  components  satisfy  the  von  Mises  yield  function 
and  the  associated  flow  rule.  On  the  other  hand  the  prescribed  blast  overpressure  at  this 
location  jumps  to  a  peak  value  of  834MPa  (120900  psi)  and  decays  exponentially  to 
621  MPa  (90000  psi)  by  the  end  of  the  run  (while  the  pressure  on  the  lower  surface 
remains  at  zero).  Consequently,  one  would  expect  that  the  mean  value  of  the  through- 
thickness  stress  component  T3  would  be  negative  (compressive)  and  that  any  tensile 
excursion  would  be  small.  For  this  reason  it  is  felt  that  the  normal  stresses  displayed  in 
Figure  1  are  exceedingly  suspect. 

Before  speculating  on  the  cause  of  this  behavior  let  us  compare  solutions  of  the  same 

physical  problem  obtained  by  use  of  other  options  of  the  PETROS  4  code.  Using  the  Kir- 

chhoff  model  with  ISTRES  —  2,  INORML  *-  0  the  distinctly  different  and  more  plausible 

results  shown  in  Figure  2  were  derived.  Rather  than  a  runaway  increase,  the  normal 
1  2 

stresses  r,  and  T2  appear  to  be  reaching  a  maximum  at  stress  levels  which  armor  plate 
may  sustain.  However,  this  Kirchhoflf  solution  has  the  drawbacks  that  (Da  two- 
dimensional  constitutive  relation  is  employed  (r3  =  0)  so  that  the  boundary  condition  on 
the  upper  surface  cannot  be  satisfied  and  (2)  the  early  time  solution  may  be  inaccurate 
since  transverse  shear  deformation  is  neglected. 

Now  consider  solutions  of  the  same  problem  derived  using  the  SHEAR  model  with 
ISTRES  —  3  where  the  three-dimensional  constitutive  relation  is  utilized  subject  to  the 
constraint  T3  =  0.  For  INORML  —  0  the  predicted  stresses  are  displayed  in  Figure  3. 
These  stresses  appear  entirely  plausible;  however,  experience  with  longer  computer  runs 
using  this  option  combination  has  revealed  a  tendency  to  unchecked  growth  in  magnitude 
of  the  variable  D3  and  the  associated  through-thickness  strain  component  y^.  For 
INORML  —  1  the  calculated  stresses  are  plotted  in  Figure  4.  These  results,  while 
different  from  the  preceding,  are  also  plausible.  In  this  case  transverse  shear  deformation 
is  permitted  but  the  "breathing"  deformation  mode  is  inhibited  by  the  non-  physical  con¬ 
straint  y3=  0.  Finally,  the  solution  for  INORML  —  2  is  presented  in  Figure  5.  It  had. 
been  expected  that  this  option  combination  would  provide  the  "best"  predictions  since  both 
shear  deformation  and  "breathing"  are  permitted.  However,  in  this  and  all  other  runs 
made  with  this  option  unstable  results,  including  negative  plastic  work,  were  obtained.  It 
must  be  concluded  that  there  exists  an  error  in  either  the  formulation  or  the  coding  for 
this  option. 

It  is  apparent  that  each  option  combination  leads  to  a  different  solution  for  the  stresses. 
The  solutions  shown  in  Figs.  1  and  5  are  obviously  unsatisfactory.  The  Kirchhoff  shell 
analysis  of  Figure  2  is  correct  within  the  limitations  of  classical  thin  shell  theory  but  a 
more  refined  analysis  including  transverse  shear  deformation  is  desired.  The  differences 
between  the  results  of  Figures  3  and  4  are  significant;  in  sequel  these  differences  will  be 
explained  and  a  more  satisfactory  formulation  derived. 
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B.  Unstable  Growth  of  D3 


As  previously  noted,  the  three  independent  components  of  the  vector  D  represent  the 
additional  degrees-of-freedom  possessed  by  the  SHEAR  option  of  the  PETROS  4  code. 
With  the  INORML  »■  0  option  the  Cartesian  components  of  D  are  determined  by  three 
equations  of  motion  derived  using  a  variational  principle.  However,  it  appears  that  this 
formulation  provides  no  material  stiffness-based  restoring  force  to  oppose  changes  in  the 
magnitude  of  D3,  the  component  of  D  in  the  direction  of  the  normal  to  the  reference  sur¬ 
face.  As  a  consequence  solutions  using  this  option  tend  to  exhibit  a  monotonic  growth  or 
decrease  in  the  magnitude  of  D3  and  of  the  directly  dependent  through-thickness  strain 
see  Figure  6.  Corrective  measures  to  circumvent  this  defect  will  be  presented  in 
Chapter  IV. 


C.  Constitutive  Relation  for  I  ST  RES  —  3 

In  the  course  of  checking  the  PETROS  4  code  it  was  discovered  that  values  of  the  nor¬ 
mal  components  of  trial  stresses  (TR  (I,  I)  in  the  code)  were  being  evaluated  incorrectly  in 
the  STRESS  subroutine  for  the  option  ISTRES  —  3.  Specifically,  the  non-zero  value  of 
DGAMMX  (3,  3)  was  being  included  in  the  calculation  of  DGAMMA,  which  is  inap¬ 
propriate  when  the  stress-strain  relation  is  constrained  by  the  condition  T3s0  .  Also, 
non-zero  values  of  corrector  stress  TC  (3,3)  (as  well  as  TR  (3,3)  and  TM  (3,3))  were 
being  used  in  the  calculation  of  X  and  (Tj')^  ,  causing  an  additional  error  in  the  elasto- 
plastic  stress  evaluation.  The  fact  that  later  in  the  cycle  T  was  set  equal  to  zero  did  noth¬ 
ing  to  remedy  the  errors  introduced  into  the  other  component  of  Tu.  Once  this  problem 
was  detected,  appropriate  corrections  were  readily  made  to  the  STRESS  subroutine. 


D.  Effect  of  Through-Thickness  Normal  Stresses 

For  most  shell  structures  subjected  to  surface  loadings  the  magnitude  of  the  through- 
thickness  normal  stress  T3  is  negligible  in  comparison  to  induced  flexural  and  membrane, 
stresses  appearing  as  the  components  Tj  and  T2  .  However,  for  the  presently  contem¬ 
plated  application  the  blast  pressure-induced  values  of  T3  during  the  early  portion  of  the. 
loading  are  of  the  same  order  as  the  other  normal  stresses  and  deserve  to  be  taken  into 
account  when  applying  the  constitutive  relations.  This  raises  the  question  as  to  whether 
such  problems  can  be  treated  in  a  rational  manner  without  resorting  to  a  complete  three- 
dimensional  analysis.  It  should  be  noted  that,  of  the  options  available  with  the  PETROS  4- 
code,  only  the  ISTRES  —  0  option  does  not  set  T3  =0  .  Unfortunately,  as  shown  in  Figure 
1,  solutions  obtained  using  this  option  predict  unreasonably  large  tensile  values  of  r3 
rather  than  the  compressive  stresses  which  would  be  expected  to  result  from  surface  pres¬ 
sure  loading.  An  alternative  formulation  for  incorporating  the  effects  of  through-thickness 
normal  stresses  will  be  given  in  Chapter  III. 


E.  Omission  of  Surface  Traction  Effects 


In  the  theoretical  formulation  report'  for  the  PETROS  4  code  the  effect  of  surface  trac¬ 
tions  was  embodied  in  the  equations  of  motion  by  terms  designated  E(kn).  Later  in  the 
same  report  it  was  argued  that  the  terms  E<k )  and  E^j  could  be  neglected  for  thin  shells. 
However  in  Appendix  D  of  Reference  1,  where  the  equations  of  motion  for  the  SHEAR 
(moderately  thick  shell)  equations  are  presented,  the  term  Eq)  is  retained  (as  it  should 
be).  It  was  discovered  that  this  term  was  not  included  in  the  finite  difference  equations  of 
motion  in  subroutine  EQUIL2  of  the  PETROS  4  code  which  are  used  to  calculate  the  com¬ 
ponents  of  D. 


F.  Reconstitution  of  Mixed  Tensor  Stresses 

In  the  cyclic  time  marching  solution  procedure  employed  by  the  PETROS  4  code  the 
values  of  the  unsymmetric  mixed  tensor  stress  components  Tk  in  each  sublayer  at  the  pre¬ 
vious  time  step  are  needed  in  the  calculation  of  elastoplastic  stresses  at  the  current  time 
step.  However,  in  an  apparent  effort  to  economize  on  use  of  computer  memory,  the  sym¬ 
metric  contra  variant  stress  tensor  components  ru  are  saved  rather  than  the  mixed  tensor 
components.  Thus  at  the  previous  time  step  the  calculations  TU  —  G^rjl  are  performed 
and,  when  stress  calculations  are  resumed  at  the  next  cycle,  the  values  of  Tk  are  reconsti¬ 
tuted  by  use  of  Tk  —  GkiTU  .  However,  in  the  interim  a  new  set  of  metric  tensors  has 
been  calculated  so  that  the  reconstituted  values  of  Tk  are  not  generally  idjntical  with  the 
values  determined  during  the  previous  cycle.  In  fact,  if  a  significant  geometry  change  has 
occurred  the  differences  may  be  appreciable.  It  is  feared  that,  for  long  computer  runs, 
these  differences  may  accumulate  to  cause  serious  departures  from  the  correct  solution. 


G.  Plastic  Work 

The  PETROS  4  code  calculates  the  total  plastic  work  performed  within  the  boundaries 
of  the  structure  as  one  of  the  ingredients  of  an  energy  balance  diagram  which  is  useful  for 
detecting  numerical  instabilities  and  for  determining  an  appropriate  time  to  terminate  the 
solution.  The  other  ingredients  are  the  total  kinetic  energy,  total  elastic  strain  energy,  and 
the  total  work  done  on  the  structure  by  external  loads.  For  conservation  of  energy  the 
sum  of  the  kinetic  energy,  strain  energy,  and  plastic  work  should  not  exceed  the  external 
work,  except  possibly  for  a  small  discretization  error.  However,  for  the  blast  loaded  panels 
of  current  interest  the  energy  balance  diagram  of  Figure  7  is  typical.  By  a  process  of  elimi¬ 
nation  it  has  been  concluded  that  the  observed  discrepancy  is  due  to  an  error  in  the  plastic 
work,  either  in  the  finite  deformation  formulation  or  the  coding.  Fortunately,  the  compu¬ 
tation  of  plastic  work  is  an  auxiliary  calculation  which  has  no  effect  on  the  basic  solution 
process. 

The  next  three  chapters  are  devoted  to  modifications  to  the  PETROS  4  code  designed 
to  remedy  the  foregoing  deficiencies. 


III.  PRESCRIBED  VARIATION  OF  THROUGH-THICKNESS  NORMAL  STRESS 


The  issues  raised  in  Section  D  of  Chapter  II  concerning  the  effects  of  through-thickness 
normal  stresses  will  now  be  examined  in  more  detail. 

A.  Stress  Wave  Considerations 

In  order  to  assess  the  significance  of  the  T3  stress  component  on  elastoplastic  calcula¬ 
tions  it  is  necessary  to  have  some  information  as  to  the  manner  in  which  this  component 
varies  through  the  shell  thickness.  Therefore,  consideration  was  given  to  an  idealized  one¬ 
dimensional  problem  of  elastic  stress  wave  propagation  for  the  through-thickness  direction. 
This  is  not  as  restrictive  as  it  may  seem.  One  is  not  concerned  with  blast  pressures  great 
enough  to  induce  plasticity  in  the  first  pass  of  the  stress  wave  through  the  thickness  or  to 
cause  spallation  off  the  far  side  of  the  shell  because  it  is  known  that  rupture  of  the  shell 
will  occur  for  lower  blast  pressures.  The  analysis  which  follows  is  for  only  slightly  more 
than  two  wave  passes  through  a  plate  and  it  is  known  that  plasticity  is  not  induced  in  the 
plate  (due  to  flexure  and  stretching)  until  much  later.  However,  it  will  be  possible  to  draw 
conclusions  which  will  also  apply  during  elastoplastic  response. 

In  the  idealized  problem  it  was  assumed  that  the  upper  surface  of  a  plate  was  subjected 
to  a  uniformly  distributed  blast  pressure  which  jumped  to  a  value  pD  followed  by  an 
exponential  decay.  The  lower  surface  of  the  plate  was  assumed  stress-free.  The  solid  line 
in  Figure  8  is  the  traveling  wave  solution  for  the  through-thickness  normal  stress  at  the 
Gauss  point  closest  to  the  loaded  surface.  On  the  other  hand,  if  one  assumes  a  linear  vari¬ 
ation  of  through-thickness  stress  from  -p(t)  at  the  upper  surface  to  zero  at  the  lower  sur¬ 
face  the  stress  at  each  Gauss  point  can  be  calculated.  In  this  manner  the  dashed  curve  in 
Figure  8  was  obtained  (this  curve  is  actually  the  "variable  mean"  of  the  traveling  wave 
solution).  Similar  results  can  be  obtained  at  the  other  Gauss  points.  The  one-dimensional 
traveling  wave  analysis  is  only  applicable  for  a  uniform  blast  pressure  which  is  not  the 
actual  distribution  for  an  explosive  charge  detonated  near  a  plate.  Inasmuch  as  it  is 
desired  to  avoid  a  general  three-dimensional  response  analysis  it  is  felt  that  the  linear  vari¬ 
ation  of  through-thickness  normal  stress  represents  a  reasonably  satisfactory  approximate 
basis  for  defining  this  component  of  the  stress  tensor  in  subsequent  calculations;  certainly 
this  is  better  than  assuming  T3  =  0. 


TRAVELING  WAVE 
SOLUTION - 


LINEAR  STRESS  VARIATION 
THROUGH  THICKNESS 
(INSTANTANEOUS 
PROPAGATION) 


B.  Modified  Constitutive  Relation 


At  the  upper  surface  of  the  plate  (or  shell)  the  boundary  conditions*  are 

rV.y.t) — G"(f“.y,t>  P  (f.y.O  .  t31-t3!-0  (I)" 

while  on  the  lower  surface  ([  —  -  h/2)  all  three  stress  components  vanish.  The  mixed  ten¬ 
sor  through-thickness  component  at  the  upper  surface  is  given  by 

r3V,y,t)  -  G3j(f*,y,t)  rV.y.O  t 

-G33(r,y,t)r33(fa,y,t)  (2) 

-  -G33(f“Al)  p(fttA° 


For  the  desired  linear  variation  of  T3  between  the  surface  boundary  conditions, 

T3V,£g,t)  -  —G33(f°,y,t)  G33^  y.t>  p(r.y,t>  <y+~)  <3> 

where  £g  are  the  locations  of  the  discrete  Gauss  points. 

In  the  PETROS  4  code  the  nine  incremental  stress  components  are  determined  from 
nine  strain  increments  Ayj  and  the  known  values  of  the  stress  components  (Tj)n  at  the 
lime  t,,  prior  to  the  incremental  change.  When  the  stress  component  r3  is  prescribed  as 
shown  in  equation  (3)  its  value  is  known  at  both  tn  and  t,^.,  so  that 

AT,3  -  (tj3)^,  -  (tj),,  (4) 

is  also  prescribed.  Thus  the  elastoplastic  constitutive  problem  is  shifted  to  determining 
eight  incremental  stress  components  and  an  incremental  strain  component  Ay|  which  gen¬ 
erally  differs  from  that  provided  by  the  ZETA  subroutine. 


These  are  the  physical  boundary  conditions.  The  displacement  model  embodied  in  the 
PETROS  4  code  does  not  provide  for  satisfaction  of  these  conditions. 

I 

Superscripts  and  subscripts  range  over  the  integers  as  follows:  greek  1,2:  latin  1,2,3. 

The  summation  convention  is  employed:  terms  having  a  repeated  index,  once  as  a  subscript 
and  once  as  a  superscript,  are  to  be  summed  over  the  range  of  that  index. 


The  generalized  Hooke  s  law  can  be  rearranged  to  obtain  the  following  expressions  for 
the  trial  stress  increments: 

ATI*  "  ' (Ayj'+vAyj2)  +  -j^AT* 

•jp 

Arl  "  +  AY2>  +  1Z^AT3  (s) 

atj  “  for  i9dj 

The  trial  stresses  at  time  t^,  are 

(r])n+i-  ^j)n+  Ar]  (6) 

To  determine  whether  plasticity  occurs  during  the  interval  t^  —  tn  —  At  the  values  of  the 
trial  stresses  are  substituted  in  the  plastic  potential  function 


T  TT  ,  T 

T  _  _i_J _ 1_  m 

^n+1  3  l^m 


-  2 


which  may  be  recognized  as  the  von  Mises  yield  function.  If  the  resulting  value  of 

T  Tj 

^irH  ^  0  the  stress  change  is  elastic  and  the  trial  stresses  (Tj)n'H  become  the  actual 

stresses  For  </>„+!  >  0  plasticity  occurs  and  the  total  strain  increments  are  com¬ 

posed  of 

Ay/  -  Ay/  +  Ay/  (8) 

i.e.,  elastic  and  plastic  parts.  The  plastic  strain  increments  are  obtained  from  the  flow  rule 

4Tj)„-ykl*iAX  <»> 

where  AX  is  a  scalar  multiplier.  Even  in  the  presence  of  plasticity  the  stress  increments  are 
related  to  the  elastic  strain  increments  by  Hooke’s  law: 


Then,  substituting 

Ay]  -  Ay/  -  Ay/ 

and  the  plastic  strain  increments  from  equation  (9),  one  obtains 


where  Alj  is  given  by  equations  (5), 


i  f  l)  1  1 — 2v  (  ml  f  3|  i 

Tj|n"lTi)n  \-p  3  lT"»In"*"'VlT3ln 


2EAX  3  T  3  C  3 

and  X  —  -j— Since  Ar3  is  prescribed  it  is  arbitrary  whether  AT3  or  (AT3)n  is  defined 

as  long  as  equation  (12)  is  satisfied.  However,  it  is  preferable  to  set 
T  j  3  C  3 

AT3-AT3  ,(AT3)n  —  0  so  that  when  the  test  for  plasticity  (equation  (7))  is  applied, 

/tA 

It  j  is  already  known.  The  stresses  at  the  end  of  the  time  interval  are 

h'L  -  (tj)„ + iT-  -  (tj'L + at;  -  x  Tj|n  -  r;]ri+i  -  *|r;jn  <u) 

In  the  mechanical  sublayer  constitutive  model  each  hypothetical  sublayer  is  treated  as 
an  elastic,  perfectly  plastic  material  having  a  distinct  yield  stress.  Consequently  the  condi¬ 
tion  4>n+i  *“  0  is  imposed  in  order  to  determine  the  parameter  X  for  each  sublayer 
experiencing  plasticity.  When  the  stresses  given  by  equation  (14)  are  subjected  to  this 
condition  a  quadratic  equation  of  the  form 

AX2-2BX+C-0  (15) 


results,  where 


c| (c|  ,  C  ]2 

1  j  1  m 

N«T'n“lK, 


T  C,  ,  T  C 

_  i  J  I  m  r 

B  -  lTJ  n+i  Ti  „  ~  T  T' 


From  equation  (15), 


T  ,  T 

rJ  -  —  tW 
rt+l  ■  JrH-i  3  m  n+1 


—  — 2 


Once  a  real,  positive  value  of  X  has  been  obtained  the  stresses  (Tj)^  can  be  calculated 
by  use  of  equation  (14).  Huffing  lo  n^  has  presented  a  procedure  for  dealing  with  complex 
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roots  should  they  arise  by  subdividing  At  for  purposes  of  stress  evaluation  only.  This  pro¬ 
cedure  was  already  incorporated  in  the  PETROS  4  code.  However,  it  has  been  found  that 
under  certain  circumstances  even  when  a  real  root  X  is  obtained  without  subdividing  Al, 
inaccurate  or  oscillatory  stresses  may  result.  This  difficulty  is  associated  with  an  exces¬ 
sively  large  excursion  of  the  trial  stress  vector  (t)  outside  the  yield  surface  in  stress  space. 
A  technique  for  coping  with  this  problem  devised  by  Huffington  (see  Appendix  B  of 
Reference  7)  has  recently  been  introduced  into  the  PETROS  4  code.  It  entails  calculating 
an  integer  L  which  defines  the  number  of  subdivisions  of  the  time  step  At,  where 


L-INT 


YLDFAC 


J 


+  <* 


)  /  cry2  —  1 


+  1 


(18) 


YLDFAC  is  a  parameter  which  varies  the  accuracy  of  the  stress  evaluation;  it  ranges  from 
0  (no  subdivision  of  At)  to  »  (differential  subintervals).  Normally  YLDFAC  “  1  is  used. 

The  foregoing  formulation,  entailing  the  prescribed  linear  variauon  of  Tj  through  the 
thickness,  has  been  incorporated  in  the  PETROS  4  code  as  option  ISTRES  —  4. 


C.  Through-Thickness  Strain  Calculation 

The  strain  increment  component  Ay^  may  be  determined  by  use  of 

“  t[4tj  -  *'[4t‘  +  4Tj)]  +  -  r!  -  Tj)nx 


(19) 


once  X  has  been  evaluated.  If  plasticity  is  occurring  at  a  Gauss  j)oint  the  value  of  Ay^  will 
generally  be  different  for  each  sublayer.  A  weighted  average  Ay|  for  the  Gauss  point  can 
be  obtained  through  multiplying  the  sublayer  Ay|  strains  by  the  same  weighting  factors 
employed  with  the  mechanical  sublayer  model  and  summing.  Alternatively,  after  the  total 
elastoplastic  stresses  Tj  at  a  Gauss  point  have  been  determined  the  value  of  Ay3  for  the 
Gauss  point  can  be  calculated  by  imposing  the  condition  of  plastic  incompressibility: 

ATk 

a>3  “1k  (20) 


Generally,  a  different  value  of  Ay|  will  result  for  each  Gauss  point  rather  than  the  con¬ 
stant  value  of  Ay^  through  the  thickness  which  the  PETROS  4  displacement  function 
admits.  In  the  next  chapter  a  compromise  resolution  of  this  discrepancy  will  be  presented. 
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IV.  MODIFICATION  OF  THE  D3  DEGREF.-OF- FREEDOM 

The  problem  of  unstable  growth  of  the  variable  D3(f°,t)  when  the  INORML  -  0 
option  is  employed  was  discussed  in  Chapter  II.  This  variable  corresponds  to  just  one  of 
the  six  degrees-of-freedom  which  appear  in  the  PETROS  4  displacement  function,  which 
has  the  following  forms: 


uUV,t> 


u0(r,t)  4-  £[N(£a,t)  -  n(r,0)]  4-  4D(f“.t) 

(21a) 

v°AQ  4-  wN  4-  £(N  -  nj  4-  £(DaAa  4-  I>iN> 

(21b) 

Y<1)k7k  +  ;Nkik  4-  CY<2)kTk  -  (rG  +  £n) 

(21c) 

The  first  form  is  a  vectorial  one,  showing  the  independent  variables  upon  which  these  vec¬ 
tors  depend.  Equation  (21b)  designates  the  components  of  the  vector  quantities  in  the 
directions  of  the  basis  vectors  Aa  lying  in  the  deformed  reference  surface  and  the  normal 
vector  N  perpendicular  to  this  surface.  The  third  form  represents  the  rectangular  cartesian 
component  version  of  the  same  displacement  functioa  Note  that 

L>Ua,l)  -  D“AQ  +  D3N  -  Y<mTk  (22) 

When  the  strain-displacement  relations  are  applied  to  the  displacement  model  of  equa¬ 
tion  (21 )  it  is  found  that  Ayj3  and  y 33  cannot  vary  with  the  through-thickness  coordinate 
£.  This  result  is  of  course  in  direct  contradiction  to  the  conclusion  regarding  the  variability 
of  strain  components  through  the  thickness  reached  at  the  end  of  the  preceding  chapter, 
where  these  quantities  were  evaluated  using  the  constitutive  relations.  Clearly,  a  generali¬ 
zation  of  the  displacement  function  to  permit  modeling  the  variation  of  y 33  with  {  would 
be  desirable.  However,  this  would  entail  an  extensive  reformulation  for  the  PETROS  4 
code  with  the  addition  of  other  degrees-of-freedom  and  increased  storage  and  computing 
requirements.  Since  this  is  not  feasible  at  present  it  appears  that  the  best  one  can  do  is  to 
modify  the  value  of  D3(f°,t)  so  that  the  value  of  by*  obtained  by  differentiating  the  dis¬ 
placement  function  will  agree  with  some  mean  value  of  the  Ay3  strain  increments  at  the 
(a  location  which  are  obtained  using  the  constitutive  relations  (i.e.,  equations  (19)  or 
(20)). 


A.  Modifications  to  Cartesian  Components 

Consider  that  the  quantities  are  known  at  a  time  and  that  the  EQUIL2  subrou¬ 
tine  of  PETROS  4  has  produced  the  next  set  of  incremental  changes  AY*2)j.  These  incre¬ 
mental  changes  also  satisfy  an  equation  similar  to  equation  (22): 

AD  -  AD°Aa  +  AD3^  -  AY^Tj  (23) 

After  modifications  to  certain  incremental  quantities  th::  following  equation  will  apply: 

m  _  m  _  m  _ 

AD  -  AD“  Aa  +  AD3  N  -  AY0*)  i; 


(24) 


where  the  terms  with_ovcT script  m  are  modified.  Note  that  the  terms  involving  An  are 
unmodified,  also  that  N  depends  only  on  Y<,)j  and  is  unaffected  by  changes  in  D.  Taking 
the  inner  product  of  both  sides  of  equations  (23)  and  (24)  with  the  unit  vector  "P  one 
obtains: 

AY<2)j  -  AD°Aa  •  ?  +  ADW  (25) 

m  _  m 

AY<2)J  -  AD°Atf  •  ?  +  AEPnJ  (26) 

Letting 

m 

tj(fa,0  -  AD3  -  AD3  (27) 

and  subtracting  equation  (25)  from  equation  (26)  gives 
m 

AY<2)*  -  AY<2)j  +  i} NJ  (28) 

The  modified  values  of  the  cartesian  components  of  D  at  the  next  time  step  are 
m  m 

“  (Y°)j)n  +  AY*2* 

-  (Y<2)j)n  +  AY*2*  +  ijN*  (jo) 

-  (Y^V,  +  tjN* 

Before  proceeding  further  a  formulation  for  tj  which  will  produce  the  desired  effect  is 
needed. 


The  derivation  of  an  expression  for  determining  ij  requires  examination  of  the  non¬ 
linear  incremental  strain-displacement  relations  employed  by  PETROS  4.  The  incremental 
through-thickness  strain  Ay  3  is  related  to  the  covariant  incremental  strains  by 

-  Ayjn.G"’3  (30) 


In  turn,  the  covariant  strain  increment  Ay 33  is  related  to  the  cartesian  components  of  the 
basis  vector  G3  and  its  incremental  change  by 

Ay33  -  jjAjJ  -  0.5  AJjAJJ  (31)* 


Since  these  are  cartesian  components  the  summation  convention  applies  even  though  the 
repeated  indices  are  superscripts. 
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For  the  SHEAR  option  (only), 

jj  -  NJ  4-  Y°>j  (32) 

AJj  -  AN*  +  AY°)j  (33) 

Using  these  relations  the  effect  of  the  previously  cited  modifications  can  readily  be  traced. 
As  noted  before,_the  surface  normal  depends  only  on  the  reference  surface  Y*1**  so  that 
modifications  to  D  have  no  effect  on  or  AN^. 

When  AY*2)j  is  modified  as  indicated  by  equation  (28)  the  components  of  the  basis  vec¬ 
tor  increment  AG3  are  affected  as  follows: 
m 

AJj  -  ANJ  +  AY<2)j  +  t,nJ  (34) 

-  AJj  +  ijNf 


Similarly, 


•  •• 

J  j  -  Jj  +  tjN1 


The  effect  on  Ay33  is,  by  use  of  equations  (31 ),  (34),  and  (35), 

m  m  m  ,  m  .  m 

Ay33-Jj  AJj -0.5  AJj  AJj 

-  (Jj  +  i]N0  (AJj  +  ijN0  -  0.5(AJj  +  ijN^)2 

-  AyM  +  Jjl^ij  +  0.5tj2 


Solving  this  quadratic  expression  for  rj , 


tj  -  -J jNJ  4- 


(jjNf  " 


2  Ay33  -  Ay33 


From  equation  (30)  it  follows  that 

Ay33  - 


Ay  3  -  Ay3,G13  -  Ay^G23 


a  ii 

and  a  similar  form  for  Ay33  .  Substituting  these  into  equation  (37): 


Ay  $  —  Ay3jG'3  —  AyjjC?23  Ay^  —  AyjiG13  —  Ay32G,J 


Up  to  this  point  no  approximations  have  been  made;  however,  since  the  quotient 
involving  modified  quantities  contains  several  unknowns  it  is  useful  to  assume 

m  m  m  m  m 

G33  ~  G33  f  Ay31G13  +  Ay32  G23  ~  Ay3,  G13  +  Ay32  G23 


On  this  basis  equation  (39)  reduces  to 


-JjNj  + 


JjN')2  - 


m 


Ay33  -  Ay^ 


G33 


Vi 


(40) 


which  will  provide  the  proper  yj  to  produce  the  desired  value  of 

C.  Implementation  of  the  Modification  to  D3 

The  normal  sequence  of  calculations  in  the  PETROS  4  code  is  indicated  in  the 
simplified  flow  chart  of  Figure  9  by  solid  lines.  Beginning  at  subroutine  EQUIL2,  where 
values  of  AY*2^  and  (Y*2,On+i  are  determined,  the  calculations  proceed  through  subrou¬ 
tines  GEOMET,  STRAIN,  and  ZETA  during  which  the  unmodified  geometric  quantities 
Nj,  AN*,  Jjj,  AJj(,  Gj<,  G*,  Ay^  and  Ayjj  *  are  computed.  In  subroutine  ZETA  there  is  a 
call  for  subroutine  STRESS  where  the  elastoplastic  stress  calculations  are  made,  after 
which  there  would  normally  be  a  return  to  subroutine  ZETA  for  calculation  of  force  and 
moment  resultants  at  tp+j.  However,  to  implement  the  desired  modification  to  D3  the  fol¬ 
lowing  changes  have  been  incorporated  into  PETROS  4. 

Once  the  elastoplastic  stress  increments  and  stresses  (T^n+i  have  been  determined  in 
subroutine  STRESS,  equation  (19)  is  used  to  calculate  the  strain  increment  Ay^  for  each 
sublayer.  The  Gauss  point  values  of  Ay3  are  then  determined  using  the  mechanical  sub¬ 
layer  coefficients.  After  returning  to  subroutine  ZETA  a  Gaussian  mean  value  of  Ay| 

through  the  thickness  is  computed  which  is  then  introduced  into  equation  (40)  as  Ay3, 
thus  completing  the  information  necessary  to  evaluate  tj  .  The  program  then  branches 
back  to  a  point  near  the  end  of  subroutine  EQUIL2  (see  dotted  path  on  Figure  9)  where 

AY(2)j  and  (Y^Op+i  are  calculated  using  equations  (28)  and  (29).  The  program  then 
proceeds  forward,  calculating  modified  values  of  j£,  AJ|[,  G^,  G*,  Ay^,  and  Ay^.  Clearly, 

,  m  i 

the  foregoing  procedure  could  be  continued  iteratively  to  cause  the  difference  Ay3  —  Ayj 


While  Ay33  is  independent  of  (  the  use  of  equation  (30)  introduces  a  very  slight  variation  of 
Ay 3  with  To  avoid  ambiguity  a  through-thickness  Gaussian  average  value  of  Ay3  is 
calculated  for  use  in  equation  (40).  A  similarly  averaged  value  of  G33  is  also  calculated  for 
this  purpose. 
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to  be  less  than  or  equal  to  some  small  quantity.  However,  experience  has  shown  that  a 
single  recycling  back  to  EQUIL2  for  the  cited  modification  is  sufficient  to  control  the 
growth  of  D3  and  provide  the  desired  compromise  value  of  y$.  This  single  recycling  step 
per  time  step  has  been  incorporated  in  PETROS  4  as  option  INORML  —  3. 


V.  ADDITIONAL  MISCELLANEOUS  PROGRAM  CHANGES 

A.  Addition  of  Surface  Traction  Terms 

The  omission  of  surface  traction  terms  in  the  PETROS  4  code  was  discussed  in  Section 
E  of  Chapter  II.  These  terms  were  introduced  in  a  general  manner  by  equation  (2.71b)  of 
the  theoretical  formulation  report.  The  specific  relations  needed  to  compute  values  of 
these  quantities  are:* 

tf,,— [V0G’'jikp]^J  (41) 

- [>/0  G31  jf  (  p]^  (42) 

The  evaluation  of  these  expressions  has  been  introduced  into  the  PETROS  4  code  through 
subroutine  S FORCE  The  finite  difference  equations  of  motion  in  subroutines  EQUIL  and 
BQUIL2  were  modified  to  incorporate  the  values  of  £(k,)  and  ,  respectively. 

It  was  hoped  that  the  addition  of  the  surface  traction  terms  would  eliminate  the  need, 
for  the  one-step  recycling  option.  Such  was  not  found  to  be  true  and  both  modifications 
are  required  for  a  satisfactory  solution,  at  least  when  the  ISTRES  —  4,  INORML  “  3 
option  combination  is  employed. 

B.  Symmetry  of  the  Contra  variant  Stress  Components 

From  equilibrium  considerations  the  contra  variant  stress  tensor  must  be  symmetric 
(in  the  absenoe  of  couple  stresses)  which  permits  the  storing  of  six  rather  than  nine  quan¬ 
tities  at  all  locations  where  values  of  this  tensor  must  be  saved  for  use  at  the  next  tjme 
step.  However,  the  calculations  in  the  STRESS  subroutine  which  provide  values  of  ru  do 
not  satisfy  this  requirement  exactly.  This  was  true  even  before  the  introduction  of  th| 
ISTRES  —  4  option,  for  which  the  problem  is  aggravated  since  the  prescribed  value  of  Tj 
generally  differs  from  that  which  would  be  calculated  using  the  symmetric  Ayu  strain 
increment  tensor.  After  the  elastoplastic  stresses  Tj,  have  been  determined  and  the  rela¬ 
tion  -  G*  Tj|  employed  the  resulting  stresses  TU  are  generally  not  equal  to  T*  for  is^j . 


*  I  am  indebted  to  my  colleague,  Dr.  J.  M.  Santiago,  Jr.,  for  providing  the  formulation  of  these 
expressions. 


In  the  original  version  of  PETROS  4  this  problem  was^dealt  with  by  ^electing  r*2,  t*3, 
and  T  as  the  correct  off-diagonal  terms  and  equating  T  ,  T  ,  and  r  to  these  quanu- 
ties,  respectively.  It  was  felt  that  this  procedure  could  bias  the  problem  solution  so  the 
program  was  modified  to  calculate  all  nine  components  TU  and  then  average  the  respective 
symmetrically  off-diagonal  components;  i.e., 

T**  -  T*)/2  for  i  j  (43) 

This  modification  was  found  to  have  a  slight  but  not  entirely  negligible  effect  upon  lengthy 
solutions. 


C.  Storage  of  Mixed  Tensor  Stresses 

The  problems  associated  with  the  reconstitution  of  the  Tj|  sublayer  stresses  each  time 
step  as  is  done  in  the  original  version  of  PETROS  4  were  discussed  in  Section  F  of 
Chapter  II.  A  revised  version  of  this  code  has  been  developed  in  which  the  stresses 
are  saved  for  use  at  the  next  time  step  rather  than  the  T**  stresses.  This  version  requires 
a  14%  increase  in  computer  memory  but  features  a  14%  reduction  in  running  time  for  a 
representative  length  run.  The  revised  version  is  preferred  because  (a)  computer  memory 
is  not  critical  today  and  any  reduction  in  running  time  is  appreciated,  (b)  the  cumulative 
deviations  from  the  true  solution  associated  with  reconstituting  the  TjJ  are  circumvented, 
and  (c)  this  version  provides  flexibility  for  future  applications  involving  material  failure. 
The  required  changes  to  the  code  are  primarily  confined  to  subroutine  STRESS,  which  is 
listed  in  the  Appendix. 


D.  Additions  to  Printed  Output 

The  format  for  listing  input  data  on  cards  for  the  original  version  of  the  PETROS  4 
code  is  presented  on  pp.  110-129  of  the  user’s  manual2.  In  sequel,  information  is  pro¬ 
vided  concerning  modifications  or  additions  to  input  data  controlling  various  options  for 
printed  output 

A  useful  feature  which  has  been  added  is  a  listing  at  the  end  of  a  run  of  the  maximum 
and  minimum  values  of  each  mixed  tensor  stress  component  along  with  the  locations  and. 
times  at  which  these  extreme  values  occur.  The  format  for  this  output  is  illustrated  in  Fig¬ 
ure  10.  For  some  applications  this  may  provide  all  the  desired  information  but  if  not  it 
directs  attenn'on  to  critical  locations  where  a  re-run  can  provide  detailed  printed  and  plot¬ 
ted  output.  Card  5  enters  the  values  of  fifteen  variables  in  format  (1515).  The  first  of 
these  variables  is  MAUXIL,  which  controls  the  printing  of  the  max/min  values:  0  —  no 
print,  l  —  print. 
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The  PETROS  4  code  controls  the  priming  of  groups  of  output  data  through  the  entries 
on  card  6a,  which  contains  the  variables 

KF,10UT(K)  (K-l.KF)  Format  (1615 ) 

KF  —  number  of  print  options  available  (currently  —  14)  and  IOUT(K)  is  the  cyclic  fre¬ 
quency  at  which  the  Kth  print  option  is  to  be  printed.  To  avoid  the  printing  of  the  Kth 
option,  set  IOUT(K)  to  an  integer  greater  than  the  final  time  step  —  ITIMEF.  The  first 
eleven  print  options  have  not  been  changed. 

Print  option  K— 12  has  been  modified  to  provide  a  rather  extensive  output  of  geometric 
and  stress  variables  which  is  useful  for  code  checking.  This  information,  a  portion  of 
which  is  illustrated  in  Figure  11,  is  provided  at  mesh  location  (IS|,IS2)  the  coordinates  of 
which  are  entered  on  card  27b  in  format  (2IS).  At  each  Gauss  point  the  following 
geometric  data  are  listed: 

J,(  =  GBASE(J,K),  AJj  S  DGBASE(J,K),  G*  =  G(J,K),  G  5  GTYPE, 

G*  3  GG(J,K),  Ay*  =  DGAM(J,K),  and  Ay|[  3  DGAMMX(J,K). 


This  is  followed  by  a  row  of  printing  which  gives  the  pressures 
p(£°,tn)  3  P(U,I2),p(fa,t,rn)  =  PPL(I1,I2)  and  the  value  of  G33  on  the  loaded  surface. 

Next,  the  code  lists  for  each  sublayer  associated  with  the  Gauss  point  the  arrays  of 

T 

(Tjj)„  3  TN(J,K)  and  Tj(  3  TR(J,K).  The  value  of  C  3  CZ  is  then  printed  as  well 
as  trY3  —  SIGMSQ.  If  C<0  the  stress  increment  in  the  sublayer  is  elastic  and 
(TkUi  =  TR(J,K).  For  C  >  0  the  stress  change  is  elastoplastic  and  the  following  infor¬ 
mation  is  printed:  A  3  AZ,B  3  BZ  ,  and  the  discriminant  B2  —  AC  3  DISCR.  This  is 
C  C 

followed  by  n  =  TCU.K),  X  —  HLAMDA  and  Tj*  3  TM(J,K).  Regardless 
of  whether  the  stress  state  is  elastic  or  plastic  the  code  then  prints 
"(REVISED)DGAMMX(3,3)  —  This  is  the  sublayer  strain  increment  Ay/  consistent 
with  the  constitutive  relations  which  is  determined  by  use  of  equation  (19).  The 
corresponding  value  of  Ay/  obtained  by  use  of  the  strain-displacement  relations  is  printed, 
with  the  rest  of  such  quantities  in  the  DGAMMX(J,K)  array  for  the  Gauss  point.  When 
the  preceding  material  has  been  printed  for  all  sublayers  at  one  Gauss  point  the. 


corresponding  material  is  listed  for  all  the  other  Gauss  points  at  the  selected  mesh  point. 

m  i 

At  the  end  of  this  output  group  the  value  of  Ay/  s  AG  AM  33  is  printed. 


Print  option  K— 13  has  been  added  which  provides  the  array  of  total  Gauss  point  mixed 
tensor  stresses  TAUF(J,K),  a  sample  of  which  is  shown  in  Figure  12.  The  frequency  of 
output  of  this  array  is  controlled  by  the  value  assigned  to  IOUT(13)  on  card  6a.  The 
locations)  at  which  these  stresses  are  printed  are  determined  by  entries  on  two  cards:  on 
card  27c  the  value  of  NUM  (—  number  of  points  at  which  mixed  tensor  stresses  are  to  be 
printed  (and  plotted))  in  format  (15)  and  on  card  27d  the  values  of  coordinate  pairs 
IPS|(1),1PS2(1)  in  format  (215)  for  I  —  l.NUM. 
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Figure  11.  Sample  of  Geometric  and  Stress  Output  Data 
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Figure  12.  Sample  Printout  of  Gauss  Point  Stresses 
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Figure  13.  Sample  Printout  of  Elastoplastic  Activity  Matrix 
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Print  option  K— 14  provides  the  matrix  of  integer  values  L  discussed  in  Section  B  of 
Chapter  III  (see  Figure  13).  These  integers  indicate  the  level  of  elastoplastic  activity 
currently  taking  plaee  at  each  mesh  location:  zero  indicates  elastic  behavior,  one  signifies 
normal  plastic  behavior,  and  any  integer  ^  2  specifies  rapid  plastic  flow  requiring  the  code 
to  subdivide  the  time  step  into  L  equal  substeps  for  purposes  of  stress  evaluation.  The 
code  will  provide  an  L  matrix  for  each  sublayer  at  each  Gauss  point  at  time  intervals  deter¬ 
mined  by  the  value  assigned  to  IOUT(14). 

E.  Additions  to  Plotted  Output 

The  PETROS  4  program  had  already  been  altered  to  couple  with  the  REPSIL  plotting 
program  (Appendix  D  of  Reference  7).  With  this  plotting  package  one  obtains  isometric 
plots  of  the  deformed  shell  surface  at  selected  time  intervals,  two  dimensional  plots  of  dis¬ 
placement  vs  time,  load  vs  time,  and  surface  strains  vs  time  as  well  as  an  energy  balance 
diagram.  The  plot  of  pressure  loading  vs  time  is  generated  for  mesh  location  (IPj,IP2)  , 
the  coordinates  of  which  are  entered  on  card  27a  in  format  (215). 

The  following  plotted  output  has  been  added: 

1.  The  cartesian  components  of  D  ,  i.e.,  Y*2)k  are  plotted  vs  time  at  the  same  mesh 
location  already  selected  for  . 

2.  Mixed  tensor  stress  components  Tj[  vs  time  (see  Figures  1-5)  are  plotted  for  each 
Gauss  point  through  the  thickness  for  the  mesh  point  location  selected  for  print  option 
K— 13. 

3.  Also,  for  the  same  location,  a  plot  of  the  through-thickness  strain  component  y|  vs 
time  has  been  added. 


VI.  CONCLUDING  REMARKS 

The  modifications  to  the  PETROS  4  code  discussed  in  the  foregoing  text  have  resulted 
in  an  improved  version  which  is  suitable  for  use  even  in  rather  exceptional  applications 
such  as  those  cited  in  the  Introduction,  where  the  hydrostatic  component  of  stress  is  a 
significant  fraction  of  the  largest  principal  stress.  The  concept  of  prescribed  through- 
thickness  normal  stress  is  considered  to  be  a  novel  approximate  procedure  for  taking 
account  of  this  stress  component  in  elastoplastic  stress  evaluations.  The  problem  of 
unstable  growth  of  through-thickness  strain  y/  has  been  successfully  circumvented  with 
the  introduction  of  the  INORML  —  3  (recycling)  option.  It  is  believed  that  the  SHEAR 
option  with  the  ISTRES  —  4,  INORML  —  3  combination  will  satisfactorily  treat  the  elas¬ 
toplastic  response  of  panels  subjected  to  blast  from  nearby  explosive  charges  and  serve  as  a 
point-of-departure  for  studies  of  panel  rupture. 


The  difficulty  with  the  ISTRES  —  0,  INORML  —  0  combination  cited  in  Chapter  II 
persists  and  is  not  alleviated  by  the  addition  of  surface  traction  terms  and  use  of  recycling. 
In  retrospect,  it  is  concluded  that  the  ISTRES  —  0  option  will  not,  in  general,  give  a  satis¬ 
factory  representation  for  stresses  in  thin  and  moderately  thick  shells  since  there  is  nothing 
in  the  basic  PETROS  4  formulation  to  enforce  the  stress  boundary  conditions  on  the  two 
shell  surfaces  (or  more  significantly,  at  the  Gauss  points  closest  to  the  surfaces).  By  con¬ 
trast  the  ISTRES  —  4  option  does  satisfy  the  surface  boundary  condition  on  the  normal 
stress  and  by  use  of  a  constrained  three-dimensional  constitutive  formulation  provides 
elastoplastic  stresses  which  appear  to  be  realistic. 

The  INORML  —  2  option  was  intended  to  take  account  of  an  average  thickness  change 
by  modifying  the  Y*2)^  variables  at  the  next  time  step.  While  the  problem  with  this  option 
(shown  in  Figure  5)  was  not  resolved  the  subject  appears  moot  since  the  new  INORML  — 
3  option  takes  account  of  Ay$  changes  in  the  current  time  step.  One  remaining  issue 
which  deserves  further  study  is  correction  or  improvement  of  the  formulation  for  the  plas¬ 
tic  work. 

The  cited  modifications  to  the  PETROS  4  program  have  affected  only  a  few  subroutines 
of  this  rather  lengthy  code;  a  listing  of  these  revised  subroutines  is  provided  in  the  Appen¬ 
dix.  Consideration  was  given  to  inclusion  of  an  application  of  the  modified  PETROS  4 
code  in  this  report  but,  owing  to  the  complexity  of  such  results,  it  is  preferred  to  pre¬ 
sent  these  as  a  separate  document. 
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NOMENCLATURE 


FORTRAN  NAME 

TEXT  SYMBOL 

DESCRIPTION 

Aa 

Deformed  reference  surface  base  vectors 

A(LA,LB) 

Covariant  components  of  metric  tensor  of  deformed 
reference  surface 

ABSCIS 

(IGAUSS,NGAUSL) 

Location  of  IGAUSSth  of  the  N  Gaussian  stations 
in  a  particular  layer;  interval  is  -1  to  +1. 

AOC(J,Il,T2) 

Y<i)j 

Acceleration  components 

ALPHA 

Coefficient  of  linear  thermal  expansion 

AGAM33 

m 

Ay/ 

Gaussian  average  of  SGAM33 

ANUM 

Status  of  material  (in  FMAT) 

AVEG33 

G33 

Gaussian  average  G33 

AVIS  (ILAYER) 

Viscoelastic  coefficient  of  ILAYER 

A Z 

A 

Coefficient  in  quadratic  equation 

B(LA,LB) 

Covariant  components  of  deformed  reference 
surface  curvature  tensor 

BM(LA,LB) 

Mixed  curvature  tensor  components  of  deformed 
reference  surface 

BSTIV(ILAYER) 

Elastic  modulus  coefficient  of  ILAYER 

BZ 

B 

Coefficient  in  quadratic  equation 

CAPQ1(LA,II,I2) 
CAPQ2(LA, 11,12) 
CAPQ3(LA,I1,I2) 

Generalized  force  resultant  tensor 

CAP2Q1(LA,I1,I2) 
CAP2Q2(LA, 11,12) 
CAP2Q3(LA, 11,12) 

Generalized  force  in  EQUIL2 

COEFHISB) 

Weighting  factors  of  the  mechanical 

sublayer  model 
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FORTRAN  NAME  TEXT  SYMBOL 


DESCRIPTION 


CONST 

Strain-rate  sensitivity  parameter 

CS(J,LA) 

See  ZETA  158-163 

CX(U) 

See  ZETA  358-375 

cz 

T 

Q^n+1 

Coefficient  in  quadratic  equation 

Cl 

Viscous  damping  parameter 

D(^a,t) 

Non-Kirchhoff  displacement  field 

Da,D3 

Tangential  and  normal  components  of 

D  in  basis  Gj 

D(J,  11,12) 

AY(1)j 

Incremental  change  in 

DA(LA,LB) 

Incremental  change  in  the  covariant  components  of 
the  metric  tensor  associated  with  the  deformed 
midsurface  of  the  shell 

DB(LA,LB) 

Incremental  change  in  the  corresponding 
curvature  tensor 

DD(J,LA) 

3AY(1)j 

dr 

DD2(J,LA) 

3AY(2)j 

d£a 

DEL 

See  ZETA  146,156 

DELBAR 

V 

See  ZETA  409 

DELNOR 

See  EQUIL2  174 

DELSN1 

See  EQUIL2  176,180 

DELSN2 

See  EQUIL2  177,181 

DELSN3 

See  EQUIL2  178,182 

DELTA(I,J) 

si 

Kronecker  delta 

DELTAP 

DELTAT 

At 

DGAM(I,J) 

A?ii 

DGAMAT 

DGAMA3 

DGAMMA 

DGAMMX(I,J) 

Ayj 

DGAM33 

&y} 

DGR\SE(I,J) 

AJJ 

DGTEMP 

DGM33 

A  y33 

DGOG 

DISCR 

DJR 

DN(J) 

AN^ 

DTAU33 

DTEMP 

DTMCI,J) 

DUM 

DX1 

DESCRIPTION 


Previous  time  increment 
Time  increment 

Covariant  components  of  strain  increment 

See  STRESS  263 

Average  Ay| 

See  STRESS  94,  95,  % 

Mixed  components  of  the  incremental 
strain  tensor 

Mixed  component  of  incremental 
strain  tensor  at  Gauss  point 

Cartesian  components  of  base  vector 
increment 

See  STRESS  93 

Gaussian  average  Ay^ 

See  ZETA  408 

Discriminant  of  quadratic  equation 

Saved  value  of  D  from  previous  time  step 

Incremental  change  in  component  of  surface  normal 

See  STRESS  123 

Temperature  increment 

Incremental  change  in  stress 

Intermediate  variable 

Increment  in  £ 1  coordinate 

Increment  in  f2  coordinate 


DX2 


FORTRAN  NAME 
DZU) 

DZAl(ll,I2) 

DZA2(I1,I2) 

DZB1(I1,I2) 

DZB2(I1,T2) 

D2(J,I1,I2) 

D33S 

EE 

EEP 

EL 

EPSL1(I1,I2) 

EPSL2(I1,I2) 

EPSU1(I1,I2) 

EPSU2(I1,I2) 

ES 

ETERM1 

ETERM2 

E1(J,I1,I2) 

E2(J, 11,12) 

FACTOR 

FMAS1 1(11,12) 
FMAS22(I1,12) 
FMAS23(I1,I2) 


TEXT  SYMBOL  _ DESCRIPTION _ 

8z  8z 
bt°'bi 

dz(A) 

bt' 

6z(A) 

8z(b) 

ai1 

6z(B> 

a? 

AY<2)j  Incremental  change  in  Y<2)j 

SeeZETA  297 

E  Young's  modulus 

See  STRESS  102 

Number  of  subdivisions  of  time  step 
Normal  strain  components  on  lower  surface 

Normal  strain  components  on  upper  surface 

E  Young's  modulus 

See  STRESS  195 
See  STRESS  196 

E/d  Surface  force  term  for  EQUIL 

f^2)  Surface  force  term  for  EQUIL2 

Coefficient  of  strain-rate  sensitivity 


Generalized  masses 


FORTRAN  NAME 

TEXT  SYMBOL 

DESCRIPTION 

FMAT 

Material  status  array 

FORCES(J) 

Component  of  externally-applied  force 
per  unit  area  in  j-direction 

FORCEZ(J) 

Convenient  grouping  of  force  components 

g3 

Covariant  basis  vector  component  of  deformed 
shell  in  direction  of  normal 

GKl,J) 

GU 

Covariant  components  of  the  metric 
tensor  of  the  deformed  surface 

GAMMAL(I1,I2) 

Shear  strain  component  on  lower  surface 

GAMMAU(I1,I2) 

Shear  strain  component  on  upper  surface 

GBASai,J) 

J,j 

Cartesian  components  of  the  base  vector 

Gj  in  j-direction 

GBTN 

JjNj 

See  ZETA  225 

GG(I,J) 

G” 

Contravariant  components  of  the  metric  tensor 
of  the  deformed  surface 

GTYPE 

G 

Metric  determinant 

HLAMDA 

X 

Plasticity  parameter 

HM(LA,J) 

Contravariant  components  of  the  relative 
moment-resultant  tensor 

HMKLA,  11,12) 
HM2(LA,I1,I2) 

Storage  of  components  of  HM 

HN(LA,J) 

Contravariant  components  of  the  relative 
stress-resultant  tensor 

HNU 

V 

Poisson’s  ratio 

HNUP 

See  STRESS  104,105 

See  STRESS  103,106 


HNUPP 


FORTRAN  NAME 
HTERM 

I 

ICOUNT 

TCI 

IC2 

IFRACT 

IGAUSS 

TGMAX(U) 

IGMTNOJ) 

IGO 

III 

IJ 

ILAYER 

INORML 

IOUT(K) 

IPLASrr(Il,I2) 

IPSI(I) 

IPS2(I) 

IRY1,IRY2, 

IRY3,IRY4 

1SB 


TEXT  SYMBOL  _ DESCRIPTION _ 

See  STRESS  241 

ik  Cartesian  unit  vector  in  k-direction 

Component  index 

Output  control  counter 

Mesh  indices  selected  by  user  at  which 
specific  output  is  desired 

Failure  model  selector 

IGAUSS  Gauss  point  index 

Gauss  point  index  where  maximum  value  of  stress 
component  occurs 

Gauss  point  index  where  minimum  value  of 
stress  component  occurs 

Selector  for  calculation  or  output  of 
max/m  in  stresses 

Output  control  index 

Number  of  components 

Layer  index 

INORML  Control  number  for  options  regarding 

Y(2)  modification 

Printout  indicator 

Plasticity  indicator 

Coordinates  of  locations  at  which  output  of 
Gauss  point  mixed  tensor  stresses  is  desired 

Indices  corresponding  to  the  limits  of  the 
complete  finite  difference  grid 


ISTRES 


ISTRES 


Sublayer  index 
Plasticity  model  control 


>RTRAN  NAME  TEXT  SYMBOL 


DESCRIPTION 


-EQ 

ISTREZ 

ISUBL 

ISKI) 

IS2(I) 

ITTM(I,J) 

ITIME 

ITIMEF 

ITTMEP 

ITIMM(I,J) 

IV 

IZ 

IZZ 

11 

I1M 

I1MAX(I,J) 

I1MIN(I,J) 

IIP 

12 

I2M 

I2MAX(I,J) 


Plasticity  model  control 
Sublayer  index 

Coordinates  of  locations  at  which  output  of 
geometric  and  stress  variables  is  desired 
(IOUT(12)) 

Time  cycle  of  maximum  value  of  stress  component 
Current  cycle  number 
Final  cycle  number 
ITIME-1 

Time  cycle  of  minimum  value  of  stress  component 
Component  index 
Gauss  layer  counter 
Upper/lower  surface  selector 

11  Mesh  point  index 
IM 

11  location  where  maximum  value  of 
stress  component  occurs 

II  location  where  minimum  value  of 
stress  component  occurs 

11  +  l 

12  Mesh  point  index 
12-1 

12  location  where  maximum  value  of  stress 
component  occurs 


12MIN(1,J) 


12  location  where  minimum  value  of  stress 
component  occurs 


FORTRAN  NAME 

TEXT  SYMBOL 

DESCRIPTION 

I2P 

12+1 

J 

Component  index 

JD1 

Number  of  stress  memory  locations  in  f1 
-direction 

JD2 

Number  of  stress  memory  locations  in  £2  -direction 

JD3 

Number  of  stress  memory  locations  in  £  -direction 

JV 

Component  index 

K 

Bulk  modulus  ”  E/(3(l-2  *>)) 

K 

Component  index 

K1,K2,K3,K4 

Boundary  condition  control  indices  on  the  four 
boundary  lines 

KF 

Number  of  print  options  available 

L 

L 

Number  of  subdivisions  of  time  increment 

LA,LB 

Component  indices 

LC 

Counter  for  time  increment  subdivision 

LEN 

See  STRESS  73 

LI1I2 

Input  to  LMAT 

LL 

Index  of  mesh  points  for  IOUT(13) 

LMAT(I1,I2,IZ) 

Plasticity  activity  arrays 

LS 

Component  index 

LZ 

See  STRESS  444 

MAX 

Integer  controlling  output  heading 

MAUXIL 

Controls  max/min  stress  output 

MIN 

Integer  controlling  output  heading 

4S 


NGAUSL 

NGAUSS(ILAYER) 

NLAYER 

NMESH1 

NMESH2 

NSBL 


Deformed  reference  surface  normal 
Undeformed  reference  surface  normal 
Number  of  Gauss  stations  in  layer 
Number  of  Gauss  stations  in  layer  —  ILAYER 
Number  of  layers 
Number  of  meshes  in  I1  -direction 
Number  of  meshes  in  £2  -direction 
Number  of  sublayers 


fS; 

NSUBL(  ILAYER) 

Number  of  sublayers  in  i-th  layer 

l- 

we 

( 

NUM 

Number  of  mesh  points  at  which  TOUT(13)  is  desired 

l 

P(I1,I2) 

P 

Pressure 

Po 

Peak  value  of  pressure 

P 

PAR 

See  ZETA  304 

PARSQ 

See  ZETA  305 

See  ZETA  376 

lv£: 

PARSQZ 

PGAM33 

yi 

See  ZETA  413 

1 

PPL(I1,I2) 

Pressure  at  next  time  step 

II 

PRSQD1 

See  ZETA  319 

3  * 

PRSQD2 

See  ZETA  320 

Ki: 

r* 

PRSQD3 

See  ZETA  321 

n 

PRSQZZ 

See  ZETA  343 

FORTRAN  NAME  TEXT  SYMBOL 


DESCRIPTION 


QCZ 

QIRCH 

QM 

QPRINT(20) 

QQQ2,QQQ3,QQQ4 

QQ1,QQ2 

QH 

QSHEAR 

QSTRES 

SGAM33 

SIGMA(ISB) 

SIGN 

SN(J, 11,12) 

SQRG 

STRESL(J,I1,I2) 
STRESP(LA,I1 ,12) 
STRJESQ(LA,I1,I2) 

SUMG 

SURFGG 


Logical  variable  for  time  step  subdivision 
Logical  variable  for  Kirchhoff  shdl  theory 
Logical  variable  for  max/min  stress  calculation 
Printout  indicator 

Logical  variables  used  in  STRESS  for  defining 
coefficients  of  constitutive  functions 

Logical  variables  used  in  EQUIL  to  avoid 
calculations  at  boundary  points 

Logical  variable  for  maximum  or  minimum  stress 
selection 

Logical  variable  for  SHEAR  option 
Logical  variable  for  stress  erasure 
r0  Undeformed  shell  reference  surface  position  vector 

SeeZETA  398 

Uniaxial  yield  stress  of  the  ISBth  sublayer 
crY  Static  uniaxial  yield  stress  of  the  material 

±1. 

Nj  Components  of  the  surface  normal 

Vg  Vg 

Generalized  forces  calculated  in  ZETA 

See  ZETA  396 

G33  at  shdl  surface 

U 

T 


TAU(I,J) 


Contra  variant  stress  components  at  Gauss  points 


IRTRAN  NAME  TEXT  SYMBOL 


DESCRIPTION 


FQ 


c 

TAUC 

^rm 

9  m 

Trace  of  the  corrector  stress  tensor 

TAUHI,J) 

i 

Mixed  tensor  stress  components  at  Gauss  points 

TAUM 

m 

'n+l 

Trace  of  the  new  mixed  stress  tensor 

TAUMAX(l,J) 

Maximum  value  of  the  stress  component 

TAUM1N(I,J) 

Minimum  value  of  the  stress  component 

TAUP(LZ,l,J) 

Storage  of  TAUF  stresses 

TAUSPH 

Hydrostatic  stress 

T  AU  S  PL(  LL,  IGAUSS) 

Storage  of  TAUSPH 

TAUSUM 

See  ZETA  400 

TAUT 

T 

rrm 

'  m 

Trace  of  the  trial  stress  tensor 

TAU11(I1,I2,IZ) 

TAU12(I1,I2,IZ) 

TAU13(I1,I2,IZ) 

TAU21(I1,I2,IZ) 

TAU22(I1,I2,IZ) 

TAU23(I1,I2,IZ) 

TAU31(I1,I2,IZ) 

TAU32(I1,I2,IZ) 

TAU33(I1,I2,IZ) 

Storage  of  mixed  tensor  sublayer  stresses 

TC(I,J) 

C 

Ti 

Mixed  tensor  corrector  stress  components 

THIC 

See  ZETA  201 

THICKN 

h 

Shell  thickness 

THICKZ 

See  ZETA  200 

THIKZ 

TIM(U) 


See  ZETA  151 

Time  of  maximum  value  of  stress  component 


YLDFAC 

YLDFAC 

Factor  controlling  subdivision  of  time  step 

YY(J,LA) 

dY(1)j 

eia 

YYU(J,LA) 

y<  i )« j 

YY2(J,LA) 

dY(2)j 

d£a 

Y2(J,I1,I2) 

yi2)i 

Rectangular  cartesian  components  of  D 

Y2D0T2(J,I1,I2) 

Acceleration  of  Y*2^ 

Y3ACELU) 

Rectangular  cartesian  components  of  N 

FORTRAN  NAME  TEXT  SYMBOL 

z  c 

Cg 

ZA(I1,I2) 

ZB(I1,I2) 

ZCEN 

ZCENTR 

zz 

ZZCEN 

>33 

T  i 

ATj 

On 

On+i 


DESCRIPTION _ 

Distance  from  reference  surface 

Distance  of  Gauss  point  from  reference  surface 

£  location  of  interface  of  upper  and  middle  layer 

£  location  at  interface  of  lower  and  middle  layer 

See  ZETA  152 

Value  of  £  at  the  center  of  a  given  layer 
Displacement  field  parameter 
See  ZETA  155 

Covariant  strain  component  in 
through- thickness  direction 

Curvilinear  coordinates  of  particles  on 
the  reference  surface 

Mixed  tensor  trial  stress  increment 
Yield  function 
Quantity  evaluated  at  t^ 

Quantity  evaluated  at  t^ 

Modified  quantity 
Integer  part  of  f  ] 


(m) 

INT  [  ] 


APPENDIX  A 


Listing  of  Significantly  Affected  Subroutines 


SUBROUTINE  AUX1L ( INOEXX ) 


IMPLICIT  LOGICALIQl 

COMMON  /CARTE/  yTEST,YNE».VS»VE 

COMMON  /CARTEL/  V(3, 20. 20). 0(3.20. 201. Y2(3. 20. 201.02(3. 20. 20) 
LEVEL  2.Y.0.Y2.02 

COMMON  /CTIME/  AUX (20) .TIME.0ELT4T.TIMEE, ITIme, tTIMEE.IXUX (20) . 

•  TOUT ( 20 ) «QRR INT ( 20 ) 

COMMON  /CTIMEL/  IRLAST (20.20) .P(20. 20) .PPL (20.20) 

LEVEL  2. IPLAST.P.PPL 

COMMON  /INOEX/  NREAO. NMRITE. NPUNCH. NMESH1 .NMESH2.N1 .N2.NZ.N1M. 6 

•  N1MM.N2MM.I 1. 12, I Z. I1ZERO.I27EBO. IBY1 .IRY2. IRY3.IPYA. ISTR1. 1ST 

•  TSTR3.ISTRA.IC1. IC2. 101 . 102. IPj . IP2. I S 1 . tS2.Ml.A2.K3.Kfc.KBUN. 

•  KZSTOP, KYTEST.I0IR.I1  TEST. I2TEST.KINITL 

COMMON  /OLOGIC/  QAUX (20) .GZETA iOSTPES.OPLAST .OSENS1 .OEOUIL  » 

•  OOIAGN.OINGEO.QlNVEL.OLOAO.QMATPR.OTMIKL.aTEMPE.aSPTEM.OAUXIl. 

•  QAUX I 2.QSPLOA .OIMPUL .OSHARP.QPESO.QIRCM.QSHEaR 

COMMON  /FBAC/  TAUE(3.3I .TAUSPM.NUM.IPSl (10) .IPS2I10) 

I I»INDEXX 

I E ( INOEX  Y  .SE.  A)  1 1 * A 

IE  ( INOEXX  .EO.  11)  H*S 

GOTO(1001.1002,1003.100*.100S).II 


1001  IEdl.EO.lCl.ANO.I2.EO.IC2)  CAU  PBINT(l) 
RETURN 


1002  CALL  PP I  NT ( 2 ) 

AOS  FORMAT  (215. 3E15. 61 

III«IT1ME/I0UT(9)A10UT!9I-1TTME 
IE(HI.NE.O.ANO..NOT.OPRINT(9)  )  GO  TO  A09 
VBITE(NPRITE.AOl)  ITTme 

ao i  format (////bx. ■•circumferential  positions  at  itime»». i5//»  n 

•2".PX."VI",13X."Y2", 13X,«Y3"1 
IE ( .NOT.OIRCM)  MR  I TE( NMRITE. 401  1) 

A 0 1  i  FORMAT(».».6AX.',Y2(l)»,10«.,'Y?(?)i.,i0X,»Y2(3>") 

OO  A 0 2  1 1 • ISTR1 . 1 5TR3 

wRITE(NVRITE.A03)  1 1 , 1 C2 . Y ( 1 . 1 1 , IC2) . Y (2 . 1 1 . IC2 ) . Y ( 3 . 1 1 . 1 C? ) 

IK  (OIRCM)  GO  TO  AO? 

MRITE (NMRITE. 140 3)  Y?(l, II ,IC?) ,Y2(?.I1.IC2) .Y2 (3. II . 1C2 ) 

1  AO 3  FORMAT  (■•«". S9X.3E1S. 6) 

AO?  CONTINUE 

MR  ITE (NMRITE .Al 1 )  ITTME 

All  E0P“AT(////l2*.«CROMN  POSITIONS  AT  ITI"Em». IS//"  II  12". 
*M  ."Yl".  13X."Y?".13X,"V3«) 

IE  (  .NOT.  OIRCM  |  «R  HE  (NMRITE.  40U) 

00  fcl2  I2.ISTR2. ISTBA 

"RITE (NMRITE. 40 3  I  ICI.I2.Y(1.IC1.I2).V(?.IC;.I2).V(3.IC1.I2> 

IE  (OIRCM)  GO  TO  .  1 2 

"RITE (NMRITE. 1  A  03 1  Y?(1,IC1.I?).Y2(?.IC1.I2).Y2(3.IC1.I2) 

Al?  CONTINUE 
4A9  CONTINUE 


AU 

IL 

1 

AU 

IL 

7 

AU 

IL 

3 

AU 

IL 

4 

AU 

IL 

s 

AU 

IL 

6 

AU 

IL 

T 

AU 

IL 

3 

AU 

IL 

9 

AU 

IL 

10 

AU 

IL 

11 

AU 

IL 

1? 

AU 

IL 

13 

AU 

IL 

lfc 

«AU 

IL 

IS 

•  AU 

IL 

16 

AU 

IL 

IT 

AU 

IL 

IS 

AU 

IL 

19 

AU 

IL 

20 

AU 

IL 

21 

AU 

IL 

2? 

AU 

IL 

23 

AU 

IL 

?fc 

AU 

IL 

2S 

AU 

IL 

26 

AU 

IL 

27 

AU 

IL 

29 

AU 

IL 

29 

AU 

IL 

30 

AU 

IL 

31 

AU 

IL 

32 

AU 

IL 

33 

AU 

IL 

34 

AU 

IL 

35 

AU 

IL 

36 

AU 

IL 

37 

AU 

IL 

39 

AU 

IL 

39 

AU 

IL 

*0 

ZAU 

IL 

fcl 

AU 

IL 

4? 

AU 

IL 

43 

AU 

IL 

44 

AU 

IL 

45 

AU 

IL 

46 

AU 

IL 

47 

AU 

IL 

49 

AU 

IL 

49 

AU 

IL 

50 

AU 

IL 

51 

AU 

IL 

52 

AU 

IL 

53 

AU 

IL 

54 

AU 

IL 

55 

AU 

IL 

56 

AU 

IL 

57 

AU 

IL 

59 

AU 

IL 

59 

AU 

IL 

60 

AU 

IL 

61 

1 1  I*!TIM£/IOUT<  10)  MOUT  <  1  0  >  —  T  T I MC 

IE  d  II.NE.O.AND.  .N0T.0PR1NTM  0)  )  SO  TO  420 

NBIT£<NNBITE.466> I TI«E . TIME , TCI , tC2 . Y d . IC1 . IC2) . Y ( 2 . tCl . IC2 I 

•  VO.IC1.IC2) 

4*6  EORHATl"  ITIM£*». IS."  TIME«".£13.6/m  POSITION  OF  OESIRFD 

•  dl»«.t2.».I2«M.I2.'*>  IS  T<1)  *».E13.S.«  Tl?)  *‘>.£13.6. 

*“  Y ( 3  >  «".£13.6> 

IE  (OIBCHI  50  TO  4?0 

*R!TE(N<iRITE.l466)  Y2d.ICl.lC?),Y2f?.ICl.lC2>.Y2d.ICl.IC2> 
14*6  F0RMAT(4SX.»Y?d>«‘'.Ei3.6."  Y? ( 2) «« ,E ) 3 . 6 . "  Y2  (  3)  «« ,E  13  .* ) 
420  CONTINUE 


1003  IE(Il.EO.ICl.iNO.  I?.  E0.IC2I  CALL  PRINTdl 
RETURN 

100*  IEd1.EO.JSl  .ANO.  I2.EQ.IS2)  CALL  PBINTlINOEXXl 
RETURN 

C  MIXED  TENSOR  STRESSES  CHECK 

ions  IIt«ITJNE/IOUT (1 3) *tOUT( 13)»IT IHE 

IEdlI.NE.0  .ANO.  .NOT.  OPBINT(13t  ISOTO  1020 
00  mil)  L •  1  .NUN 

IEdl.EO.IPSl  (L)  .ANO.  I2.E0. IPS21L) 1GOTO  IMS 
miO  CONTINUE 
GOTO  1020 

inis  call  PRINT  1 1 NOE XX 1 
1020  RETURN 
END 


AUXIL  6 2 
AUXIL  S3 
AIJXIL  S* 
AUXIL  SS 
,  AUXIL  SS 
AUXIL  ST 
POINT AUXIL  *3 
AUXIL  SR 
AUXIL  TO 
AUXIL  T 1 
AUXIL  72 
AUXIL  73 
AU»IL  74 
AUXIL  TS 
AUXIL  TS 
AUXtL  77 
AUXIL  TS 
AUXIL  TO 
ALIXIL  SO 
AU«lL  SI 
AU»IL  S? 
AUXIL  S3 
AU*IL  34 
AUXIL  OS 
AUXIL  SS 
AUXIL  »T 
AU»IL  SS 
AU»IL  SO 
AU«tL  RO 
AU»IL  01 
AUXIL  02 


-:v:vv». . »■  v ■££  >K\  ^  Viv/l . 


an  n  r>  ;•)  o  o  r»  r>  r>  o  n  o  o  oooono 


SUBROUTINE  EOUIL  EOUIL  1 

EQUIL  ? 
EOUIL  3 

EVALUATE  THE  DISPLACEMENT  INCREMENTS  ANO  EBOM  THEM  EVALUATE  THE  EQUIL  * 

NEW  POSITIONS  EOUIL  5 

EQUIL  6 

IMPLICIT  lOSICALIO!  EOUIL  T 

EQUIL  9 

COMMON  /CASTE/  TTEST.VNEx.VSAVE  EQUIL  9 

COMMON  /CARTEL/  T ( 3.20 .20 ) .0 (3.20 . 20 ) . 72 ( 3.20 .20 » .02 ( 3 . 20 .20 )  EQUIL  10 

LEVEL  2. Y. 0.72,02  EQUIL  11 

EQUIL  l? 

COMMON  /CTIME/  AU* (20) .TIME, OELTAT. TIMER, ITIME.ITIMEE. I*u*(201 .  EQUIL  13 

•  IOUT ( 20  I ,QPB t NT ( 20 )  EQUIL  14 

COMMON  /CTIMEL/  I  PL AST (20,201  ,P(20,20) ,PPL (20*20)  EOUIL  IS 

LEVEL  2.1PLAST.P.PPL  EOUIL  IS 

EQUIL  IT 

COMMON/CTIMEB/ITIMEC.ITIMEP, DELTAS. OELX.OHR.UNH.hEE  EOUIL  10 

•  , TREES, HTMIK.UEINIS.OEINO.TSTART.YSTABT.VOOTE  EQUIL  19 

•  .ES.8STIV (*1 ,NSTIV(*1  EOUIL  20 

EQUIL  21 

COMMON  /DAMP/  MOANS. QANPE.DEAeT.TDAMP.TOTKP, Cl  EQUIL  22 

EOUIL  23 

COMMON  /EORCB/  EORCEZ ( 3)  EQUIL  2* 

EQUIL  2S 

COMMON  /INOEX/  NBEAO.NMRITE.NPUNCH.NMESHI .NMESH2.N1 ,N?«NZ,N1M,m2M,£QU!L  2S 

•  N1 MM.H2MM. I 1 . I?. 12. I {ZERO. I21EBO. IBV1 . IPT2. IRT3. IRTA. ISTB1.ISTB2 .EQUIL  2T 

•  ISTB3.ISTRA. IC1 .IC2.I01 . 102. IP1 .IP2.IS1 . IS2.K1.K2.K3.KA.0BUN.  EOUIL  2S 

•  K2STOP.KVTEST.I0IB.I1TEST.I2TEST.KTNITL  EQUIL  29 

EQUIL  30 

COMMON  /PUSH/  EORCESI3) .VELOCI3) .RATIO. RATION. 0X1. 0X2. TEMP. OTEMP,  EOUIL  31 

•  ESP ACE. TSPACE .EINCNO.ESTOP. TSTOP. THCOEE  EOUIL  32 

COMMON  /PUSML/  SQRAT (20.20) . SQR A2 (20.20) .EMAS11 (20.20) .  EQUIL  33 

•  EMAS22(?0.20I .EMAS23(20.20) .EMAS33 (20 . 20)  EQUIL  3A 

LEVEL  2.S0RAT.SQRAZ.EMAS11 .EMAS22.EMAS23.EMAS33  EQUIL  3S 


•  ESPACE. TSPACE. EINCNO.ESTOP. TSTOP. THCOEE  EOUIL  32 

COMMON  /PUSML/  SQRAT (20.20) .SQR A 2 (20.20) .EM AS 11 (20.20) *  EQUIL  33 

•  EMAS22(20.20) .EMAS23(20.20) .EMAS33 (20 . 20)  EQUIL  3A 

LEVEL  2.S0RAT.SQRAZ.EMAS11 .EMAS22.EMAS23.EMAS33  EOUIL  3S 

EOUIL  36 

COMMON  /Q8C0N0/  Q1.Q2,Q3.QA.QS.Q6.QT.QR.Q9.Q10.Q11.Q12.Q13.014.Q1SEQUIL  3T 


•  .OEREFl , QERFE2.QES£E3.QEPEEA,0EC0PN,QQ1 ,0Q2 

•  TOAC1.QTOAC2.QTR-C3.QTRACA 

COMMON  /SUPNOM/  SNPQ (3) 

COMMON  /SURNOL/  SNI3.20.20) 

LEVEL  2.SN 

COMMON  /S2/  STRESS (3) 

COMMON  /S2L/  STRESL (3.20.20) .STRES0I2. 20.20) .STRESP(2. 20. 20) 
LEVEL  2. STRESL. STRESQ. STRESS 

COMMON  /TNCOMP/  HMl (2,20,20) .HM2(2. 20.20) .CAPQ1 (2.20.201 . 

•  CAPQ2 (2 . 20*201 .CAPQ3 (2*20.20) .CAP201 (2.20,20) . 

•  CAP202(2.20.20I .CAP2Q3 ( 2, 20.20) 

LEVEL  2. HMl .HM2.CAPOI.CAPQ2.CAPQ3.CAP2Q! .CAP202.CAP203 

COMMON  /VELS/  Y3ACELMI  .Y3DOTC3)  .V3ACPO(4,3) 

COMMON  /VELSL /  VELOI3.20.20I . ACC ( 3 .20 . 20  I 
LEVEL  2.VEL0.ACC 

COMMON  /SUPEOP/  El  (3.20.20) .£2(3.20.20) 

LEVEL  2.E1.E2 


TFPMTaOELTAT/CFLTAP 

QQ1-I1 .EO. IPV1 .0R.I1.E0.IBV3 
0Q2«I2.F0. IRT2.0R.I2.EQ. IPV4 

1 1 »*I 1 ♦ 1 

I1M.I1-1 

I2P«I2.1 

I2M«I2-l 

L1»1SIGN(1  .Ni-2»H  ) 
L2»rSIGN(l.N?-2»I2) 

TERM.) ,/EmaSI 1 ( II. 12) 
TERM2«0ELTAT40ELTAT 
00  2 9P  J« 1 . 3 


.UEOUIL  39 
EQUIL  39 
EOUIL  40 
EQUIL  41 
EOUIL  42 
EOUIL  43 
EQUIL  44 
EQUIL  4S 
EOUIL  46 
EQUIL  47 
EQUIL  49 
EQUIL  49 
EQUIL  SO 
EQUIL  51 
EOUIL  52 
EQUIL  53 
EQUIL  54 
EOUIL  55 
EOUIL  56 
EOUIL  57 
EQUIL  59 
EQUIL  59 
EOUIL  50 
EQUIL  61 
EOUIL  5? 
EQUIL  63 
EQUIL  64 
EQUIL  65 
EQUIL  66 
EQUIL  57 
EOUIL  69 
EOUIL  69 
EOUIL  70 
EOUIL  71 
EQUIL  72 
EOUIL  73 
EQUIL  74 
EQUIL  75 
EOUIL  76 
EQUIL  77 
EQUIL  79 


59 


DISPLACEMENTS  POP  viscous  DAMPING 
DjPbDIJ.I1 .12) 

NEXT  STATEMENTS  AVOID  IMPERFECTIONS  DUE  TO  BOUNO-OPP  EPBOBS 

IP (KI.E0.7.AN0.J.E0.1)  GO  TO  298 
IF<*2.EO.T.ANO.J.EQ.2)  GO  TO  298 
IE<Kl.EQ.5.AN0.11.£Q.2.AN0.J.EQ.l)  GO  TO  298 
IFIK2.EQ.5.AN0.I2.EQ.2.AN0.J.EQ.2)  GO  TO  298 
IFOO.EQ.S.ANO.Il.EQ.NlM.ANO.J.EQ.l)  GO  TO  298 
IP IKA.EQ.5.AN0. I2.EQ.N2M.ANO. J.EO.?)  GO  TO  298 
IP  (Ml .E0.6.AN0.J.EQ.1 1  GO  TO  298 
IPIK2.E0.6.AN0.J.EQ.2)  GO  TO  298 

OEPIVsO. 

OXaOXl 

0XPACTa.5/DX 
IPtOOl)  GO  TO  111 
GO  TO  (102.104.106) .  J 

102  OEBIVbOEBIV.OXFACT4(CAPQ1 (1,I1P,I2)-CAP01 I1.I1M.I2)  1 
GO  TO  119 

104  OEBIVbOEP I V<OXFACT4<CAPQ2< i.IlP. 12) -CAPQ2 < i . 1 1M. 12) ) 

GO  TO  119 

106  0EBIV«0EPIV»0XPACT4<CAPQ3<1.I1P. I2)-CAPQ3( 1.IIM, 12) ) 

GO  TO  119 


FOOIL  79 
E0U1L  80 
E6UIL  81 
EOUtL  92 
EQUIL  83 
EOUIL  84 
EQUIL  85 
EQUIL  85 
EQUIL  87 
EQUIL  98 
EQUIL  99 
EQUIL  90 
EOUIL  91 
EQUIL  92 
EQUIL  93 
EQUIL  94 
EQUIL  95 
EOUtL  96 
EOUIL  97 
FQUIL  98 
EOUtL  99 
EQUIL 1 00 
EQUIL 1 0 1 
EQUlLl 02 
EOUIL 1 03 
EQUlLl 04 
PQUlL 1 05 
EQUIL106 
EQUIL107 
FQUIL108 
EQUlLl 09 
EQUlLl 10 


111  OUMbOXEACT«L1  EQUIL109 

GO  TO  ( i 12. 114. 1 161 .  J  EOUlLllO 

1 12  OEBIVbOEPIV.OUM6(-3.»c*PQ1 <1.11,12) <4.*CAPQ1 < 1 . I 1»L1 . 12) -CAPQ1 t 1 . I EQUlLl 1 1 

•  1 ‘2*L1 . 12 ) )  E0UIL119 

GO  TO  119  EQUIL 1 1 3 

114  OEPIVbOEPIV.OUM6 (.3.*caP02<1 .11. 12) *4.4CAPQ2< 1 , 1 1 »L1 . I21-CAPQ2 (1 , IEOUIL1 14 

•  1 *2#L 1.12) )  EQUIL 1 15 

GO  TO  119  EQUlLl 16 

116  OEPI  VbOEPI  V»OUM4  <-3.*CAPQ3<i,Il,I2).4#*CAPQ311,Il*H.I2)-CAPQ3<l,IE0UlL117 


•  1 *2*L1 . 12) ) 
119  CONTINUE 


OX>OX2 

OXFACT  » • 5/DX 

IE (QQ2)  GO  TO  131 

GO  TO  (122.124.126).  J 

1.22  OEBIV«OEPIV.OXE*CT4<ctPQl  <2,11,  I2P) -CAPO  1  <2.11.12“)) 
GO  TO  139 

124  0EBIVB0ERIV-0XEACT4<CAPQ2<2,li . I2P) -CAPQ2 <2. 1 1  .  I2M)  ) 
GO  TO  139 

126  0ERIVb0EPIV«0XPACT4<CAPQ3<2,I1,I2P)-CAP03(2.I1.I2M) ) 
GO  TO  139 


EQUlLl 18 
EQUlLl 1 9 
EQUlLl 20 
EQUIL 121 
EQUIL 1 22 
FQUIL123 
EQUIL124 
EQUIL 1 25 
EQUlLl 26 
EQUIL 127 
EQUIL128 
E0UIL129 
EQUlLl 30 
EQUlLl 31 
EQUlLl 32 
EQUlLl 33 
EQUIL134 
FQUIL I  35 
EQUlLl 36 
EQUlLl 37 
EQUlLl 38 


ni  0UMb0XPACT»L2  EQUlLl  37 

GO  TO  (132. 134, 136). J  EQUIL138 

132  OERIV»OEBIV<OUM4<-3,#cAPoi <2.11 . 12) .4.4CAPQ1 <2. 1 1 . 12»L2) -CAPQ1 (2.IEQUIL139 

•  I . 1 2»2*L2) )  EQUlLl  40 

GO  TO  139  EOUIL 1 4 1 

134  OEP I VbOERI V.OUM# (-3, *C APQ212.il, I?)<4.6CAPQ2(2«I1*I2<l2) "CAPO 2 (2.IE0UIL142 

•  1 . I2<2*L2) )  E0UIL143 

GO  TO  139  F0UIL144 

136  OPR IVbOEBI V»DUM* (— 3,»CAPQ3 (2.11,12) ♦4«*CAP03( 2. Il,I2<L2)*CAPQ3(2.IEQUlL145 

•  l.I2*2*L2)l  EQUlLl 46 

139  CONTINUE  EQUIL 1 4 7 

FQUIL148 
EQUIL 1 49 

nPRIV«OEBIV.FOPCEZ <J1  FQUIL 1 50 

ACC< J. II . I2)b<OEPIV.E1 (J.li.I?) )4TEBm  EQUIL 151 

0<J.ll.I?)«0<J.Il.I2)4TE9MT.iCC<J.Il,I2)4TEP“2  EQUIL 152 

TIJ.Il . I2)»V<j,[l.I2).0(J.li.I2)  EQUIL 153 

VISCOUS  DAMPING  Cl  EQUIL154 

IP  < TOAMP  .EO.  0.01  GOTO  298  E0UIL155 

0<J.I1,I2I»0<J.I1.I2)-<0<J.I1.I2)»OJP)4C1  EQUIL 156 

298  CONTINUE  EQUIL157 

RETURN  EQUlLl 58 

PNO  EQUIL l 59 


.  .  9FV-. 


u  u  o  o  o 


SUBROUTINE  EOUIL2  EQUIL2  I 

EQUIl?  ? 
EQUIL2  3 

EVALUATE  MIDSURFACE  GEOMETRIC  QUANTITIES (BASE  VECTORS.  SURFACE  EOUIL?  A 

NOPMAL.  metric  TENSOR.  CURVATURE  TENSOR.  ETC.)  EQUIl?  s 

EQUIL?  * 

IMPLICIT  LOGICAL (Q)  EQUIL2  7 

EQUIL?  S 

COMMON  /CARTE/  VTEST.YNEV.VSAVE  EQUIL?  9 

COMMON  /CARTEL/  7(3.20.20) .0(3.20.20>.Y2(3.20.20),02(3.20.20>  EQUIL210 

LEVEL  2, Y .0 • Y2.02  EQUIl?ll 

F0UIL21? 

COMMON  /CTIME/  AUX(?fl) .TIME .OELTAT.TrMEE, TTIME. ITIMEE.IAUAI20) .  E0UIL213 
«  IOUT (20) .QRRINT (201  E0UIL214 

COMMON  /CTIMEL/  IPLAST(20.?0).P(?0.20).PPL(?0.?0)  FOUIL215 

LEVEL  ?. I  BLAST  « P  *  PPL  FOUIL216 

EQUIL717 

COMMON/CTIMEB/ITIMEC.ITIMER.OELTAP.OELX.OmO.UNM.HEF  FQUIL219 

•  .TKfEP.MTHIK.OFINIS.OFINP.TSTABT.YSTART.YOOTF  EQUIL219 

•  ,ES.RSTIV(4I .NSTIV(A)  FQUI1220 

FQUIL2?l 

COMMON  /ORTERM/  MOISPL  !*>  .M0!SDR(4.3>  .PERIOtA.S)  .OISP).  (A. 3)  .Q0ISP3EQUIL22? 
COMMON  /OPTERL/  PBSCPR I  3. ?0 .20)  EQUIL223 

LEVEL  2.PRSCRR  E0UIL224 

F0UIL22S 

COMMON  /INDEX/  NREA0.NNRITE.NPUNCH.NMESH1 .NMESH2.N1 ,N2.NZ,N1M,n?m,E0UII??6 
A  NIMM.N2MM, n . 12. IZ. n  ZERO. I2ZER0. IBY1 . IRY2. I0Y3. IRYA. ISTBJ . ISTB2.E0UIL22T 

•  ISTB3.ISTRA.IC1.IC2.I01.I02.IP1.IP2.IS1  •  IS2  .« I  ,K2  ,K3,  KA.KBIjN,  E0UIL229 

A  KZSTOP.KYTE ST.IOIR.il  TEST. 1 2TEST.K INITL  EQUIL229 

EQUIL230 

COMMON  /OPTION/  mauXIL.MINGEO.MINVEL.MLOAO.MMaTPR.MSPLOA.  E0UIL231 

A  MSPTEM.MTEMPE.MTHIKL.MIMPUL.ISTRES.INORML.ISTREZ  F0UIL232 

E0UIL233 

COMMON  /P0LEGM/SNP0LEIA.3) .VELPOL (4.3) .NUMS ( A ) ,NPL .L0CP01 1  A) .QP0LEFQUIL23A 

•  .OYPOLF (A. 3) .YROLE (4.3) .02P0LE(A,3) .YPPOlE (A. 3) . JPOLE(A)  FQUIL23S 

•  . QPOLET  FQUIL236 

E0UIL237 

COMMON  /PUSH/  FORCES(I)  .VELOCO)  .RATIO. RATION, 0X1. 0X2. TEMP. OTEMP,  E0UIL23B 

•  FSPACE.TSPACE .EINCNO.ES TOP. Ts TOP. THCOEF  EQUIL239 

COMMON  /PUSML/  SQRAT (20.20) .SORAZ (20.20) .FMAS11(20.?0) .  EQUIL2A0 

a  EMAS2?(?0.20).EMAS23(20.20).EMAS33(20.20)  EQUH241 

LEVEL  2.SQRAT.S0RAZ.EMASII.EMAS22.EMAS23.EMAS33  F0Ull?»2 


EVALUATE  MtQSUREACE  GEOMETRIC  QUANTt TIES (BASE  VECTORS.  SURFACE 
NORMAL.  METRIC  TENSOR.  CURVATURE  TENSOR.  ETC.) 


IMPLICIT  LOGICAL (Q) 


COMMON  /CARTE/  YTEST.YNEV.YSAVE 

COMMON  /CARTEL/  Y(3.?0.20) .0 (3.20 . ?0 ) «Y2 (3 .?0.?0 ) .0? ( 3.20.20 ) 
LEVEL  2.Y.0.Y2.02 


COMMON  /CTIME/  AUX(20) .TIME. OELTAT.TIMEE. (TIME. ITIMEE.IAUX(20) . 
*  IOUT (20) .QPRINT (20) 

COMMON  /CTIMEL/  IPLAST (20.20) .P(?0.20) .PPL (20.20) 

LEVEL  2. IPLAST, P. PPL 


COMMON/CTIMEB/ITIMEC.ITIMER.OELTAP.OELX.OhP.UNM.HEF 
a  .TKFEP.HTHIK.QFINIS.QFINP.TSTART.YSTART.YOOTF 
•  ,ES.RSTIV(AI .NSTIV(A) 


COMMON  /OPTION/  mauXIL.MINGEO.MINVEL.MLOAO.MMaTPR.MSPLOA. 
a  MSPTEM.MTEMPE .MTMIKL .MIMPUL . ISTBES . I NOR ML . ISTREZ 


•  FSPACE.TSPACE .EINCNO.ES TOP. Ts TOP. THCOEF  EQUIL239 

COMMON  /PUSML/  SQRAT (20.20) .SORAZ (20.20) .FMAS11(20.20) .  EOUIL240 

a  FMAS2?(?0.20).FMAS23(20.20).FMAS33(20.20)  EQUH241 

LEVEL  2.SQRAT.S0RAZ.FMASII.FMAS22.FMAS23.FMAS33  F0UIL?»2 

E0UIL2A3 

COMMON  /QRC0N0/  Q1 .Q2.Q3.QA.as.06.QT.QR.Q9.Q10.on .Q12.013.Q1A.Q1SEQUIL2AA 

•  .OFBEEl .QFBEE2.QFBEE3.QFREFA.QFC0PN.QQ1 .002  .QEQUIL2»S 

•  TRACI, 0TBAC2.QTRAC3.0TPAC4  F0UIL2A6 

E0UIL2AT 

COMMON  /SUPNOM/  SNPBI3)  E0UIL2AR 

COMMON  /SUPNOl/  SN ( 3 .20 .20 )  E0UIL2A9 

LEVEL  2.SN  EQUIL2S0 

EQUIL2S1 

COMMON  /S2/  STRESE ( 3 )  EQUIL2S? 

COMMON  /S2L/  ST9ESU3.20.20).STPESQ(2.20.?0).STRESPI2.?0.20)  EQUIL253 

LEVEL  2. STBESL . STRESQ .STRESP  E0UIL2SA 

E0UIL2SS 

COMMON  /TENCOM/  YY (3.2) «YYY (3.2,2) .A (3.3) .B (3.3) .AA(3,3) .QR ( 3. 3) .  EQUIL2SS 

•  RM(3,3) ,0A(3«3) .OR (3.3) .G(3.3) *GG(3.3) *0N(3) .00(3.2) ,000(3.2,2) .  F ou I L?S7 

•  OGAM (3, 3) .OGAMMX (3.3) .HN<3.3) .HO (2) ,TAU (3.3) . TAUSPL (3.3)  E0UIL2SS 

•  .002(3.2)  .0002(3.2,2)  ,  YY2  (  3. 2)  •  YYY2  ( 3. 2. 21  .  WU(3.2)  .ORMI3.3)  FQUIL2S9 


COMMON  / S2/  STRESE (3) 

COMMON  /S2L/  ST9ESU3.20.20)  .STRESQ  ( 2. 20 . 20)  .  STRESP  (  2.20 .20 ) 
LEVEL  2. STBESL. STRESQ. STRESP 


COMMON  /TNCOMP /  HMl  (2,20.20)  .MM2(2,20,?'>)  .CAP01  (2.20.20)  . 


E0UIL2SR 

FQUIL2S9 

EQUIL2S0 

F0UIL2M 

E0UIL262 


•  **  .  •  .  *  A^A* 


•  CAPQ2(2.20#20)  .CAP03 (2.20.20)  .CAP2QI  (2.20.20) . 

•  CAP 202(2. 20.20) .CAP203 ( 2. 20 .201 

LEVEL  2.HM1.MM2.CAPQ1  .CAP02.CAP01.CAP2QJ  .CAP2Q2.CAP203 

COMMON  /W6LS/  Y3ACEL(3>.V30OTt3>.Y3ACP0(4.3> 

COMMON  /VELSL/  VEL00.20.20)  .ACC  (  3.2  0.20 ) 

LEVEL  2.VEL0.ACC 

COMMON  /ZACCEL/  SN00T2O.20.20)  . Y20OT2 ( 3 .20 . 20) 

LEVEL  2.SNOOT2.Y200T2 

COMMON  /DELC/  ICOUNT 
COMMON  /OEL/  DELBAR (20.20) 

LEVEL  2. OELBAR 

COMMON  /SUREOR/  El (3.20.20) «F2(3.20. 20) 

LEVEL  2. El. £2 
TERMT«OELTAT/OELTAP 
001*11 .E0.IRY1 .OR.ll .E0.JBY3 
002*12 .EO. IBY2.0R. 12.F0. JUY* 

I1P*11.1 
1 ?P*l 2* 1 
tiM*n-i 

t?M«t2-l 

L 1 « l S I GN ( 1 ,N1 -2*1 1 ) 

L2*tStGN(l.N2-2*I2) 

TERm2*0ElTAT*0ELTAT 
IF (  I  COUNT  .EO.  OlGOTO  *00 
CALL  NaCEL 
00  294  J*l . 3 

NO  Y2  STRAIN  ON  00  SLIOING-Ct AMPEO  EOGE  IN  nOISDR(k.J)* 
IE (001 , ANO .I1.E0.1RY1. A NO .*1  .EQ.9.AN0. MO ISOR(l.J)  .EO.  1 ) 
IF ( QQl . ANO . T I .EO. IRY3.ANO.K3.EO.O.ANO.MOISDR(3«J) .EO. 1 1 
IF  (002.  ANO.  :?.’•> .  I  BY  2. ANO. *2. FQ. 9. ANO. MO  I  SDR ( 2. J)  .EQ.l ) 
IF (0O2.AN0.I?.c0.IBVk.AN0.F4.E0.4.AN0.MOIS0R(*.J) .EO. 1 > 

TO  A VO 1 0  ROUNDOFF  IMPERFECTIONS 
IF  (*1 .EQ.T.ANO.J.EO.1 I  00  TO  300 

IF  IM2.E0.7.AND.J.E0.2)  GO  TO  300 

IF  (*1 .EO.S.ANO • 1 1 .F0.2.AND.J.E0. 1 )  GO  TO  300 
IE1K2.EQ.S.AN0.I2.EQ.2.AN0.J.F0.2)  GO  to  300 
IFOO.EO.S.ANO.l 1.E0.N1M.AN0.J.E0.1)  GO  TO  300 
TF  (F4.F0.S. ANO . 1 2 .EG , N2M .ANO .J.EG.2)  GO  TO  300 
IF  (F1.E0.6.AN0.J.E0.1I  GO  TO  300 

IF  (F2.EQ.A.AN0. J.E0.2)  GO  TO  300 

OFBIVaO. 

Ol*OAI 

OxFACTa.S/OX 

IF (001 )  GO  TO  130 

GO  TO  (120. 124. 1201 .  J 

120  OFBIV*OFOIV.OXEACT*(CAP?Q1  (l.UP.I2)-CAP201  (I.I1M.I2)) 
GO  TO  140 

1  24  OE*IV*OFBI V»DXF  ACT*  (  CAR 202  U  .1  IP.  I  2)  -CAP 202  (1  .  I1**.I?)  > 
GO  TO  140 

124  0fBtV*0EBIV.0XFACT*(CAP2O3(l.IlP.I2)-CAP2O3(l.llM.I2) ) 
GO  TO  140 


130  Ot|M*OXFACT*Ll 

GO  TO  (132.134,1341 .  J 

I  3  2  OFBIV*OEB1V.OUM*  (-3.*CAP201  (  i  .  1 1  ,  12)  -4. *CAP20 1  ( 1  .  1 1  *L  1  . 
»  1 . 1 1  *  2*L 1 -12) ) 


E0UIL243 
EQUIL244 
E0UIL245 
E0UIL244 
E0UIL24T 
FQUIL24B 
EGUIL240 
FOUIL2TO 
E0UIL2T1 
E0UIL2T2 
E0UIL2T3 
F0UIL2T4 
E0UIL2T5 
E0UIL2T4 
E0UIL2TT 
E0U1L2TS 
E0UIL2TB 
FQUIL2B0 
EQUIL2AI 
E0UIL242 
E0UlL2*3 
E0UIL204 
E0UIL2*S 
EOUIL204 
EQUIL2ST 
E0UIL24B 
E0UIL2A0 
E0UIL290 
EQUIL201 
E0UIL202 
1  0IBECTI0NE0UIL293 
GO  TO  390  E0UIL294 
GO  TO  390  EQUIL29S 
GO  TO  390  F0UIL294 
GO  TO  390  E0UIL29T 
E0UIL294 
E0UIL299 
EOUILIOO 
E0UIL101 
EOUIL192 
EOUIL103 
E0UIL1 04 
EOUILIOS 
E0UIL1 04 
EQUIL10T 
EQUIL10B 
E0UIL1 09 
EQUIL110 
EOUIL 1 1 1 
F0UIL1 12 
EOUILll 3 
EQUIL 1 1  * 
EQUIL115 
E0UIL1  14 
E0UIL11T 
E0UIL1 I* 
EOUIL 1 1 9 
EOUIL 1 20 
EOUIL 121 
EOUIL 1 22 
E0UIL1 23 
EOU I L 1 2* 
E0UIL1 2S 
E0UIL1 >4 
I21-CAP701 (EQUIL12T 
FOUR)  2* 


.  •**  •  •  *  %* 


coy  co  o  o 


c 

c 

c 


60  TO  140 

114  OFP I V*0EP I V*0UM#  (-3.4CAP202U.il.  I?)  »A.4CAP?02tl.Il.Ll.l2) 

•  i.ii.?4li.i2)) 

GO  TO  140 

116  n?PIV»0€PIV«OUM*(-3.«C*P2Q3(l.tl.I?)»4.4C»P203U.n*Ll.I2) 

•  1  .IW^Ll.I?)  ) 

140  CONTINUE 


ox«ox? 

0XEACT«.5/0X 

IF(0Q2)  SO  TO  151 

60  TO  (14?. 144, 146).  J 

14?  0£PIV40EBIV«OXF*CT»(C»P?01 <?. II. I?P)-CAP?01 (?. I1.I?N) I 
60  TO  160 

144  Oe»IV*OeoiV«OXFACT»(CAP?Q?l?.!l.I?P)-CAP?Q?(?,tl.I?N! 1 
60  TO  1 6 0 

146  nePIV»0EPlV*0XP»CT»(C»P?03(?.Il.I?P)-CAP?O3(?. Il.IEMl ) 
60  TO  160 


151  0UM*0xFACT*L? 

60  TO  (IS?. 154. 156) .J 

IS?  0fOIV40EPIV*0UM«(-3.»CAP?01 (?. II . I?) »4.»CAP?01 (?. 1 1 . I?*C2>  • 

•  ?.  1 1  ,  !?♦?•(.?)  > 

60  TO  160 

1S4  0EPtV40EPIV40UM*(>3.*CAP?0?(?.Il,I?)»4.*CAP?0?(?.Il.I?.L2)> 

•  ?.n.i?*?»L?>  i 
GO  TO  160 

1«G  nE»IV40eoiV*OU«*(-3.«CAP?03(?.U.l?)»4.»CAP?03(?.Il.I?*L2)- 

•  2.11, I?.?*L2) ) 

160  CONTINUE 

Y  200T2  U.  II.  I?)  *  (OEBl  V*E?  ( J.  II.  I?) -STBESLU.il.  12) -EMAS23U 

•  *V3ACEL ( J) ) /ENAS22 (11.12) 

02U.I1 . 121*02  ( J.  1 1 . 121  •TEPNT»TEPN?4Y?00T2  (J.  1 1.12) 
Y?(J,I1.I?)*Y2(J,I1,I2).02(J.I1,I2) 

3*0  CONTINUE 
?0B  CONTINUE 

TAKE  OUT  NOPNAL  COMPONENT  WHEN  IT  IS  DESIPED  TO  8E  ZE»0 
IEIINOPML.NE.l )  60  TO  300 

OELNOP»SN  (i,Il.I2)*02(l,ll,I?)*SN{?«H,I?)*02(?.Il.I?) 

•  .SN(3. II. 12)402(3. II. 12) 

OELSNl «0ELN0P*SN (1,11,12) 

OELSN?»OELNOB*SN (2,11,1?) 

0ELSN3«0ELN0P*SN (3. 1 1 « 1 2) 

GOTO  410 

400  OELSNl  4-0EL8AB  ( 1 1  ,  12)  «SN  ( 1  .  1 1 . 12) 

0ELSN2.-0EL8AB (II, I?) #SN ( ?. 1 1 . 12) 

OELSN3«-0ELSAR ( 1 1 , I?) #SN (3, 1 1  . 12) 

410  02(1.11.12) «02 ( 1 . 1 1 . I ?) -OELSNl 
0?(2.H.t2)«02(2.Il.T?)-OELSN2 
0?(3.I1.I2)«0?(3.I1.I?)-0ELSN3 
Y200T2(1,I1.I2)«Y?OOT2(1,I1.I?)-TERM240FLSN1 
Y200T2 ( 2, 11.12) 4Y200T2(2, II, I?)-TEPM?40ELSN2 
Y200T2 (3. 1 1 . 12) »Y200T2 (3. 1 1 . I ? I -TEPM?*0ELSN3 
Y2(1.I1.I2)4Y?(1,I1.I?)-0ELSN1 
Y?(?.I1,I2)4Y2(2,I1,I?) -OELSN? 
Y?(3.I1.I2)«Y2(3.I1.I?)-0ELSN3 
300  CONTINUE 
PFTUPN 
ENO 


E0UIL1 ?6 
-CAP?Q2(EQUtLl30 
E0UIL13! 
EQUIL1 3? 
-CAP?Q3(EOUIL133 
E0UIL1 3* 
E0UIL13S 
E0UIL136 
F0UIL13T 
EQUIL13S 
E0UIL136 
E0UIL1 40 
E0UIL141 
E0UIL1A? 
E0UIL143 
E0UIL144 
E0UIL145 
EOUIL 1 46 
EOUIL 1 4T 
EOUIL 1 46 
E0UIL146 
E0UIL1 SO 
E0UIL1S) 
E0U1L1S? 
F0UIL1S3 
EOUIL 1 S4 
-CAP201 (E0UIL1SS 
EOUIL 1 56 
F0UIL1ST 
-CAP?Q2(F0UIL1S8 
EOUIL 1 56 
EOUIL1SO 
-CAPP03 (EOUIL 161 
E0UIL16? 
F0UIL163 

E0UIL164 
II. I?)  E0UIL16S 
E0UIL166 
E0UIL1ST 
FQUIL166 
EOUIL 1 69 
FOUIL1TO 
F0UIL1T1 
E0UIL1T? 
FOUIL 1 73 
E0UIL174 
EOUIL 1 T5 
E0UIL1T6 
EOUIL l 77 
EQUIL1 7® 
FOUIL) 70 
E0UIL1 60 

EOUIL 1 6 1 
EOUIL 1 6? 
FQUIL153 
E0UIL164 
EOUIL 1 65 
FOUIL 1 66 
EQUIL167 
FOUIL 1 66 
EOUIL10B 
KQUIL160 
EOUIL 1 61 
FOUIL 16? 
F0UIL163 
F0UIU64 


£1 

& 

& 

■•■A; 

ft 


,*  • 


k- 


< 


SUBROUTINE  PRINT ( INDEX X) 
IMPLICIT  LOGICAL ( 0 ) 


COMMON  /ALLENE /  TOTAL. TOTk IN. TOTElA . TOTPL A. TOTME X . TOTTFM, INERGY 

•  .TOTVIS.TOTE1 .TOTE2.ECHECK 

•  . OTP { 3 . 1 ) .S0BG.S0BA.SIGHS0.A7.BZ.cz 


COMMON  /CABTE/  YTEST.YNEM.TSAVE 

COMMON  /CARTEL/  V ( 3. 70 .70 ) .0 ( 3 .20 . ?0 > . V? (3 .20 .?0 I . 02 ( 3 . 20 .20  I 
LEVEL  E.T.0.T2.02 


COMMON  /CTImE/  AUA (201 . TIME. OELTAT.TIMEP, ITIME. ITIHEF. IAUX l?OI . 
•  l OUT < 20 ) .OPBINT (201 

common  /CTImEl/  IBLAST(20.20I .PI20.Z0)  .PPLI20.20I 
LEVEL  Z. IPLAST.P.PPL 


COMMON  /IN0EX/  NBEAO.NMPITE.NPUNCH.NMFSNl ,NmeSh?.N1 .NP.NZ.NIM.NPM 

•  N1 mm.nzmm, 1 1 . 12. iz. I IZE oo. IZZEBO. I0V1 . IPr2. IBT3. IPV*.  ISTBl ,  ISTBJ 

•  ISTR3.  ISTB4.ICI.  IC2.I01 ,102. IP! .  IP2  .  I  SI . IS2 .* 1 .*2 ,K3 . *  4 .KRUN. 

•  EZSTOP.EVTEST. IOIB. I 1TEST. I2TEST.FTNITL 


common  /PMVSCN/  EE, HNU. alpha. const. expon.factop. bate. PHO. 
•  mlahOA.COEFE (SI .SIGMA (51 . TN ( 3.31 • TC 13.31 .DELTA (3. 3) 


CQMMON/POLEGM/SNPOLE (*.3 1 .VELPOL (4,31 .NUMS (4) .NPL.LOCPOL (4) .OPOLE 

•  .OVPOLE (4.3) . YPOLE (4.31 .02P0LE (4.3) .T2P0LE (4,3) . JPOLE (4) 

•  .QPOLET 


COMMON  /QBCONO/QI  .02.03. 04,05. Q6. <1T. (JB.0B.01  0.01 1 .0I2.Q13.QI4.01S 

•  .0EOEE1 .0EBEE?.QFBEE1.0FOEE4.0FC0BN.Q(Jl,Q02  . 

•  TBAC1 « 0TPAC2 . 0TBAC3 . OTP AC4 


COMMON  /QlOGIC/  QAUX(20) .Q7ETA.OSTBFS.OPLAST.OSENS1.0EOUIL. 

•  OOIAGN.arNGEO.QINVEL.OLOAn.OMATPO.QTMtKL.OTEMPE.QSPTEM.QAUXll , 

•  QAUX  17,  OSPL  0  A  «  Q  I  MPilL  .OSHARP.  QPESO.QI  PCN.OSHEAB 


COMMON  /SURNOM/  5X08(11 
COMMON  /SUONOL/  SN(3.70,70) 
LEVEL  2.SN 


COMMON  /S?/  STRCSEOI 

COMMON  /S7L/  STPESL I  3,26.20) . STRESO I ?. 20 . 20 ) . STPESP ( 2. 20 . 20 ) 
LEVEL  2. STPESL. STRESO. STPESP 


common  /TCNCOM/  YV (3.2) «  YY Y (3.2,7) ,A  13,3) ,R(3.3> , A  A ( 3 , 3  I .BBI3.3) , 

•  BM( 3. 3) ,0A ( 3.3) .0" (3.3) .G ( 3. 3) ,GG (3, 3) .ON ( 3) .00 ( 3.2) .ODn ( 3.2.71 . 

•  OGAM (3, 3 | .OGAMMX (3,3) ,hn(3.3I ,mQ  (7) ,TAU(3,3) . TAUSPL (3,3) 

•  .002(3.2) .0002(3.2.7) . YY2 ( 3 .2) . YYY? 1 3. 7 . 7) . YYU I  3. 2) .DBM (3.3) 


COMMON  /ThknS/  7.7Z.7CENTR.THTCEN.ABSCIS (4,4) .HEIGHT (4,6) . 
•  NGAUSS(A) , I  GAUSS , NL  AYER. ILAYEP .NSURL ( 4 ) • I SUBL 


COMMON  /TNCOMP/  HM1 (7.70.201 .hM2(7. 70.70) .CAPQ1 (7.70,70) . 

•  CAPQ2 ( 7,70 . 20 ) .CAPO  3 (7. 70. 20) .CAP201 (7,70.20) . 

•  CAP202  (7,20.20  I  «CAP7f)3  ( 2 , 7n  ,?o  • 

LEVEL  7 , hm 1 ,hm?,CAP01 .CAP02.CAP03.CAP201 . CAP7Q2, CAP7Q3 


COMMON  /STPOUT/  ISR.L.LC.0tSCB,TP(3.3) .TN(3. 3) 
COMMON  /GPP  I  NT  /  GTypc.gbaSE  (  3.3)  .OC.BASE  (  3,  3) 
COMMON  /TwPEF/  AGAM33 


POINT 

1 

POINT 

? 

PP  I  ST 

1 

PRINT 

4 

PP  \  NT 

s 

PRINT 

6 

PP  I  ST 

7 

PPIST 

5 

PRINT 

9 

PP  J  NT 

10 

POINT 

U 

PP  1ST 

1? 

POINT 

13 

PP  I  NT 

\4 

PP  I  NT 

IS 

PP  I  ST 

1* 

,PP [NT 

17 

•  PP I  NT 

ia 

PPIST 

19 

PRIST 

?0 

POINT 

?1 

POINT 

?? 

POINT 

n 

PPlNf 

24 

point 

ps 

point 

?6 

point 

27 

point 

point 

?9 

OPPINT 

10 

point 

31 

point 

3? 

PRINT 

13 

POINT 

34 

point 

35 

point 

36 

point 

3T 

POINT 

39 

point 

39 

POINT 

40 

POINT 

4  l 

POINT 

4? 

POINT 

43 

©PINT 

44 

POINT 

4  5 

POINT 

46 

POINT 

47 

POINT 

44 

POINT 

49 

POINT 

SO 

POINT 

SI 

POINT 

s? 

POINT 

S3 

point 

S4 

POINT 

ss 

POINT 

56 

POINT 

57 

POINT 

5« 

POINT 

59 

00  I  NT 

SO 

POINT 

si 

POINT 

SP 

64 


3V>  .3'.3  . 


y  j. 1 


W-W  A'  "•  '•  3-  V  u  1  ‘ 


"•T 


l  nno 


i  noi 

TO  1 


COMMON  /FRAC/  TAUF (3.3) .TAUSPM.NUM. IPS) M0)»IPS2(10> 
COMMON  /MATRIX/  YLOFAC.ANUM 

COMMON  /MA  TP  t  L /  LMAT (20,20. 16) ,FMAT (6.20.201 
LEVFL  2.LMAT.FMAT 

IF  (ITIME  .£0.  -1)  RETURN 
II«IN0FXX 

IF  { I  NOE  XX  .SE.  4)  II»4 
GOTO! 1000.2000.3000.4000) . I! 


POSITION  AND  TINE 

I I ImITIME/ tOUT (l )*I0UT ( 1 1-ITINE 
IFIIII. NE.O.ANO.. NOT. OPRINTO  I  )  GO  TO  T01 
NR  I  TE  ( NMR ITE.1001)  11.12. I  TINE 
FORMAT(»lli»".IS."  I2«".I5." 

CONTINUE 


ITInEn“.I5/) 


PRINT  03 
PRINT  64 
PRINT  OS 
PRINT  00 
PRINT  07 
PRINT  OS 
PRINT  OR 
PRINT  TO 
PRINT  71 
PRINT  72 
PRINT  73 
PRINT  7  4 
PRINT  7S 
PRINT  70 
PRINT  77 
PRINT  78 
PRINT  7R 
PPINT  SO 
POINT  SI 
PRINT  82 
PRINT  S3 
PRINT  84 
PRINT  SS 
PRINT  SO 
PRINT  87 


C  FIRST  ANO  SECONO  Y"S  DERIVATIVES 

C  FIRST  ANO  SECONO  OY"S  DERIVATIVES 

IIIpITIME/IOUT (?)*IOUT (21-ITTME 
IFIIII. NE.O.ANO.. NOT. QPRINTI?) I  GO  TO  702 
WRITE(NWRITE.IOO) 

100  FORMAT (//RX. "FIRST  PARTIAL  DERIVATIVE  OF  Y".20X. "SECONO  PARTIAL  DEPRINT  88 
•PIVATIVE  OF  Y"//12X. "Y Y(J,1)",8X«"YY(J,2)".IRX," YYY(J»1.1)".5X.  PRINT  8R 
•«VVY(J,1.?)«.5X."YYV(J, ?,?)"/)  PRINT  RO 

WRITE (NMR I TE. 102)  YYCl .1  I ,VV(i .2) .YYVI 1 .1.1)  ,YYV  t 1 . 1 , ? 1 . Y YY ( 1 , 2 , 2 1 PP I  NT  31 
MRITF(NMRITE.102)  YYI2.11 .YYI2.2) .YYY (P.1.1) .YYVI2.1 .2) , YYY ( 2. ?. 2) PRI NT  R2 
Vfl ITE I NMR I TE . 1 02 1  YY13.1 I . YY ( 3,2) , VVV ( 3. 1 . 1 1  . Y YY ( 3 . 1 . 2 1 . YYY ( 3 . 2. 2) PP I  NT  R3 
102  FORMAT (7X.2E15.6.12X.3E1S.0)  PRINT  R4 

WRITE(NWHITE.103)  PRINT  R5 

1(13  FOBMAT1//OX. "FIRST  PARTIAL  DERIVATIVE  OF  DELTA  Y»,  1 4X,  "SECOND  PARTPP I  NT  90 
•IAL  DERIVATIVE  OF  DELTA  Y"// l 2X. "OD ( J. 1 1 ». 8X ,"DD I J. 2) " . 1 9» • "000 ( J. PRINT  RT 
*1,1)".5X,"000(J,1«2)",5X,  "OOD  ( J.2.2)  ••/)  PRINT  RS 

WRITE  (NWRITE.  102)  00(1.11.00(1,2)  .DODO  .1.1)  .ODDI  1.  1.2)  .ODD  ( 1 .2.2)  PR  I  NT  RR 
WRITE (NWR ITE. 102)  00(2.11.00(2,2) .000 (2. 1 . 1 ) .000 ( 2 . 1 . 2) .000 (2.2. 2) PRI NT  1 00 
WRITE (NWRITE . 1 02 1  00(3.11 .00 ( 3.2) . DOO ( 3. 1 . 1 ) .000(3.1.2) .000 ( 3. 2. 2) PR  I NTl  01 
IF  (OIRCH)  GO  TO  5101  RRINT102 

WRITE (NWRITE. 5100)  PRINT103 

8)00  FORMATI//RI. "FIRST  “ARTIAL  DERIVATIVE  OF  Y2". 1RX. "SECONO  PARTIAL  DPR  I  NT  1 0  4 
•PRIVATIVE  OF  Y2».//1?x."YY2(j,1)«,7X,»yy?(J.2)",18X,»yyy2(j,1,1)",point105 

•  4X,»Yvy7(J.1.2)».4X,"YYY2(J.2,2)»,/)  PR  I  NT  1 00 

WRITE  (NWRITE.  1021  YV2  (1,11  .YY2I1.2)  .VYY2U  ,1.  1)  «  YYV2  (  1  .  1.2)  ,YYY2(1  PP  I  NT  1  07 

•  .2,2)  PR  I  NT  1  OS 

WRITE (NWR ITE. 1021  YY2I2. 1)  ,YY2 ( 2.2) , YYY2 (2. 1 . 1 ) . YYY2 ( 2 , 1  . 2 ) . YYY2 ( 2PB INTI  OR 

•  .2.2)  PPINT110 

WRITE (NWR ITE. 102)  YY2 ( 3 . I ) . YY? ( 3. 2) . Y YY2 ( 3. 1 . 1 ) , YYY2 (3,1,2)  . YYY2 ( 3PB I  NT  1 11 

•  .2.2)  PR  I  NT  11? 

WB ITF ( NWRITE .5103)  PR  INTI  13 

5103  F0PMAT(//5X. "FIRST  PARTIAL  DERIVATIVE  OF  DELTA  Y2", 14X. "SECOND  PARPR INTI  14 
•TIAL  DERIVATIVE  OF  DELTA  Y2"//12X»"00?(J.l)".7X, "DO? (J»?)",1BX, "OOPP I  NT  1 1 5 
•02 ( J.l .1 )».4X,"000?(J.l .2)".4X,»0002( J, ?.?)"/)  PR  I NTl 10 

WRITE (NWR ITF. 102)  DO? (i,i) .002 ( 1 .2) .0002 ( l • 1 . 1 ) .0DD2 ( 1 . 1 .2 ) .0002 ( 1PB INTI  1 7 

•  .2,?)  PR  I  NT 1 1 S 

WRITE (NWRITE. 1021  002(2,1) ,002(2.2) .0002 ( 2. 1 . 1 ) .0002 ( 2. 1 . 2) .0002 ( 2PR I NTl 1 9 


•  .2.2) 


PRINTI20 


WRITE (NWR ITF. 102)  002(3.11 .002(3.2) .0002(3.1  .1) .0DD2 ( 3 . 1  .2 ) . 0002 ( 3PR INTI  21 


8101 

702 


»  .2.2) 

CONTTNUE 

CONTINUE 

FIRST  ANO  SECONO  METRIC  TENSORS 

FIRST  ANO  SECONO  METRIC  TENSORS"  INCREMENTS 


PRINT!?? 
PP I  NT  1 23 
PR  I  NT  1  24 
PR  I  NTl 25 
PBINT1 20 
BRINT1  27 
"PINT)  28 


1 1  I*IT  IME  /  I0UT  ( 3 )  •10))T  ( 3)  •  ITTME 
IFIIII. NE.O.ANO.. NOT. OPRINTO)  )  GO  TO  703 
WRITE(NWRITE.IOO) 

Too  F0RMATI//14X. "FIRST  METRIC  TENSOR". ?RX . "SECONO  METRIC  TENSOR"/) 

WRITE  I  NWRITE. 9071  A(i.l) .4(1,21.8(1. 11.811.2) 

RO 7  FORMAT (?X,"A( 1 .1 ) "".FI 3. 6."  4(1.2) "".El 3.6. OX ."8(1.1) ■" .E 13.6. 

•  "  8 (1.2)«", El  3. 6) 

WRITE (NWR l TE.R08 )  4(2, l). 4(2. 2). R(2. 11.8(2,2) 

BOS  FORMAT (21. "A (2, 1 ) *".E13.6."  A ( 2, 21 «".E1 3.6.6X."8 (2. 1 ) *".E1 3.6. 

•«  8(2.21 »«.E13. 6) 

WRITE (NWRITE. 10R) 

109  FORMATI//RX, "FIRST  METRIC  TENSOR  INCREMENT". 1RX. "SECONO  METRIC  TENPB I  NT  1 40 
•SOR  INCREMENT"/)  PBINT1A1 

WRITE (NWRITE .9091  OA ( f , 1 1 , DA ( 1 .2) .08 ( 1 . I ) .08 ( 1 .2)  PPINT14? 

ROR  FORMAT!"  OA ( 1 , 1 ) •» . E 1 3 . 6 , "  0 A ( 1 . 2 ) »" .E 1 3 .6 . »  OR ( 1 . 1 ) «".E 1 3.6.PR I  NT  1 4 3 


PP  I  NT  1  23 
PRINT1 30 
PR  INTI  31 
PR  I  NT  13? 
PP INTI  33 

PR  I  NT  1 34 
PRINT! 38 
PRI NT  1 36 
PR  INTI  37 
PR  INTI  38 
PRINTl 39 
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cu  o  coou  uoo 


•  »  0B(1.?I«“.E13.6)  MISTI** 

MBITC(NMBITE.BIO)  0AC2.1) .0A(2. 2) .08(2.1 > .08(2.2)  PRINT145 

«10  FORMAT  ("  0A(2.l)o",Ei3.6,"  0A(2.2)o“.fl3.6.»  08(2.1)»“«E13.6»PBIsn*6 

•  "  08(2. PI*". £13.4)  PRISTl*r 

7ft3  CONTINUE  PR  1ST  1 *0 

PR  INTI *9 

STRESS  RESULTANT  BRIST150 

PBIST151 

IlIotTtME/lOUTISI ‘lOUTlSI -ITTME  PR  1ST  ISP 

IF(III.NE.O.ANO..NOT.OPRINT(5) I  80  TO  70S  PR  I  ST  I  S3 

SPITE (NMBITE. 11*1  PR  I  ST 1 5* 

118  FORMAT!//"  CONTRAVAflIANT  COMPONENTS  OF  THE  RELATIVE  STBESS-RESULTAPR I ST1 55 

•NT  TENSOR"/)  PRINT156 

MRITC(NpRITE.RlT)  hn  ( I ,1 ) ,HN (I .21 .MM (1 .3)  PRINT1S7 

PIT  F0RMAT(1X."HN(1.1)»".E13.6."  HN ( 1 • 2) ■"•El  3. 6* "  HN ( 1 • 3 ) ■" • El  3.0)  BPINT158 
MRITEINSBITE.R18)  NN ( ?, 1 ) .HN 12 ,2) .HN ( 2 ,3 )  PRIST15R 

RIB  FORMATUX,«HNI2.1)»".E13.6.«  HN(2.2) «".E1 3.6."  HN  ( 2. 3)  ■"  «E  13  .  A)  PPTNT160 

70S  CONTINUE  PBINT101 

RETURN 


RESULTANTS 


2000  IF  tlTIHE.EQ.O)  RETURN 

III«ITIME/IOUT (6) *IOUT (6) -ITTME 
IF(III.NE.O.ANO..NOT.QPRINT(6) )  80  TO  706 
ITIMEP»IT!ME-1 

timep«itimep«oeltat 


MBITE  TIME  ANO  ITIME  AT  PRECFEDIS8  STEP 


WBITE(NSRITE.2I  TINE®. ITIME® 

MO I TE  ( SMB ITE.920) 

920  format I//36X. “TENSOR  OEFInEO  IN 
•ONENTS  OF  TME"/33X.»  THE  PETROS 
•ESULTANT  TENS0R"//3X."II 
•  •  12) "/) 

00  120  11*1 .N1 
00  120  I2«1.N2 


PRIST16? 
PRIST163 
PR  I  ST  1 6* 
PR  I  ST  1 65 
PRIST1 56 
PR I  ST  16  7 
PR  I  ST  1 65 
PBINT169 
PR  I  ST  1 70 
PR  INTI  71 
PRIST17P 
PRIST1 73 
PR  I  ST  1 7* 
PRIST1 75 
PR  1  ST  17  5 

».29X ."CONTRA VARIANT  COhPRRIsTI 7T 
A  REPORT". 36X. "RELATIVE  MOMENT  RPP 1  ST  1 75 
I2",32X."0(LA.J.I1  . 12I".*5X.«H(LA.LB.I1PPIST17<) 

PR  I  ST  1 "0 
PR  1ST  I* I 
PR  I  NT  1 5P 


MB ITE( NMBITE .921)  II . 1 2. CAPO  1 (1,11.12) .CAPQ2 ( 1 . 1 1 . 12 1 .CAPQ3(1. II .  PRIST  153 


PRIST15* 
PRIST185 
PRIST186 
PR  I  ST  1 9  7 
PRIST 199 


•  12) .HM1 (1 .11, 12) .HM] (2, I) , I?) 

021  FORMAT  ( 2 1 5  ,  »  CAPQM.1)«".E13.6."  CAPO  (I  .2)  »".£  1  3.6. 

•  "  CAPO  1 1 .3)  ■"•E13.A."  HM1 1  ■'•  ,£  13.6."  hh12»».E13.61 

MRITE (NmBITE. 922)  CAP01 (2. 11 . 12) .CAP02 (2. 1 1 . 1 2 ) . CAP03 ( 2 . II . 12 ) . 

•  HM2(i, ii. I®), HM2(2. 11,121 

922  FORMAT (13X. "CAPO (2*1 (•".E13.6,"  CAPO (2. 2) ■"•£ I  3.6,"  CAP0(2.3)»"  P0IST19R 

•  .El  3.6,"  MM21«".E1 3.6."  mh22.«.E1 3. 6 )  PBIST190 

IF  (OIRCH)  GO  TO  129  PRIST19) 

MRITE (NMRITE.021 )  1 1 . 1 2. CAP  201 ( 1 . II . 12) .CAP202 ( 1 . 1 1 .  1 2) . CAP203 ( 1  . 1  PR  I  NT  1 0  2 

•  1 .12) ,STRESP( 1.11 .12) .STRESP(2,I1 .12)  BPIST193 

021  F09MATI2I5."  CAP 20 (1,1 )■"•£) 3. 6,"  C*P20 ( 1 .2) •".£ 1 3.6."  CAP2Q ( 1 . 3 > PR  I  NT  1 9* 

P (2) ■" .E 1 3 .6)  PP I  ST  1 95 


•*" ,E 1 3.6 , 


P  ( 1 )  ■"  ,E  1  3 .6  ." 


MBITE(NMRITE.822ICAP201 (2 , 1 1 . 12) ,C AP202 ( 2, 1 1 . 12) .CAP203 ( 2. 1 1 • 12 1 .  PRIST  196 

•  STRESOI 1 .11 .12) .STBFSQ(2,I1 .12!  PRINTJ97 

022  FORMAT  I  1 2X ."CAP20 ( 2. 1 ) »".E1 3.6,"  CAP20 ( 2. 2 ) ■" .E 1 3.6.  PPIST199 

•«  CAP20 (2. 3) ■",£  13.6."  0 ( 1  ) ■",£ 1 3 .6."  0 (2) ■"•£ 1 3 .6 1  POINT199 

MRITE (NMRITE.023)  1 1  .  1 2 .  STRE  SL  ( 1 « 11 . 1 2 ) .STRESL 12> 1 1.1?) .STRESL I  3.  PBINT200 

•  11.12)  PRIST201 

023  FORMAT  (215."  STRESL  ( 1  )■"•£!  .3,6,"  STRESL  ( 2)  ■"•  E  1 3. 6.  »  STRE  Si  ( 3 1 POI ST202 


•»"«E1 3.6) 
128  CONTINUE 
766  CONTINUE 


NORMAL 


II  UITIME/IOUTITIMOUT  (71 -ITTME 
IF(!II.NE.O.ANO..NOT.QPBINT(T) )  SO  TO  707 
ITIMEP. ITTME. 1 
TIMEP»lTIMEP*OELTAT 


MBITE  TIME  anO  ITIME  AT  PBECFEOINO  STEP 


MRITE (SPRITE.?)  TIMFP, ITIMEP 
MR  I T  E ( NmBI TE.1602) 

106?  F0RMAT(//17X. "COMPONENTS  OF  th£  SURFACE  S0RMAL"//3X ," 1 1 
•6X . "5N ( 1 1 " . 1 6X ,"SN ( 2) " • 1 flX  »"SN ( 3) "/ ) 

DO  1007  I  1 ■ 1 , N 1 
00  moT  I2-1.N? 


12". 


PRIST203 
PBIST20* 
PRIST205 
PRIST206 
PRIST207 
P0IST269 
PRIST209 
PBIST210 
PR  t  ST  2 1  1 
PRIST21 2 
PR  I  ST? 1 3 
PRIST21* 
PRIST21 5 
RR I  ST? 1 6 
PRIST21 7 
PRIST218 
ORIST219 
DRIST220 
PR  I  ST?  2 1 
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.•  V. 


r'  ■ 


WPITEINWRITE.  10031  1  1,12.  SN< l , 1 1 , 12) ,SN (7. 1 1 . 12) .SN ( 3. 1 1 . 12) 
1003  F0RMAT(2IS.3£i5.6> 

1007  CONTINUE 

IF(.NOT.OR)  GO  TO  70T 
00  1010  t*l.NPL 
K.LOCPOL(I) 

foio  WRITE(NWR!TE.1020i  K,  SNPOLEIX.I) .SNPOLE (K.2) .SNPOLE (K.3) 

1020  FORMAT!"  POLE  ON"/"  SI0E».I2."  ■  ••.3E15.A) 

TOT  CONTINUE 

:  POSITIONS  ANO  INCREMENTS 


PPINT22? 

PBINT223 

PPINT22A 

PRINT22S 

PRINT22S 

PRINT22T 

PBINT228 

PPINT229 

PRIST230 

PBINT731 

PP1NT23? 

PBINT233 

PPINT734 

PBINT235 

PRINT23ft 

PPINT23T 

PBINT23S 


II  I»ITIMC/IOUT 18) •lOUTISI-ITTME  PPINT734 

IE ( III.NE.O • ANO, • NOT  .OPR I  NT ( R ) )  GO  TO  TOR  PRINT23S 

WRITE (NWRITE.7I  TIME.ITIME  PRINT236 

2  FORMAT Cl" .RX."T  IME»».E1 3.0,8*"  ITImE»m.IS/)  PRINT23T 

WRITE (NWRITE.3)  PRINT23B 

3  FORMAT (/48X, 'CARTESIAN  COOBO INATES ' //30 X .' POSI T ION' . 37X .' CHANGE  INPBINT239 

•  POSITION* //3X. • II  I2<  .RX. <Tl  •  ,  13X. < Y2# . 13X, • T3» , 1RX, >01 • . 1 3X , • OPRINT24  0 

•2* , 13X,'03* , 12X. 'P* • 1  OX  , ' IPLAST*'/)  PPINT2*! 

WRITF  (NWRITE.4)  (<U, 12, Til. II, 12), YI2.I1, 12), TO, II, 12), Oil, II. 12)  PR  I NT24  2 

•  .0(2, II, 12), 0(3,11, 12) .Pill. 12), I PL AST (11,12) ,  PPINT243 

•  l ?• IRY2, 1RY41.I1 ■ IRY1 , IPY3)  PRINT244 

4  FORMAT  (2IS.3E15, 6,  '  <  .4£ IS. (SOX .  15 )  PBINT245 

IF  (OIRCH)  GO  TO  1004  PRINT24G 

WRITE (NWRITE. 1053)  PRINT24T 

10S3  FORMAT  (////52X,"Y2  COORDINATES  "//30X ."POSI TION", 37X ."CHANGE  IN  PRINT24R 

•POSITION"//  3X 1 1  I2".ftX."Y2(l)".t0X."Y7(2)».10X."Y2(3)".15x.  PR  I NT749 

•"02(11". 10X."02(2)".10X. "02(3)  IPLAST  •"/)  PRINT250 

WRITE (NWRITE. S) ( ( II . I2,Y2(1 , II , 12) ,V2(2, 11,12) ,V2(3. II ,12) .02(1 . I1PRINT2S1 

•  .12), 02(2. II. 12), 02 (3. II. 12). IPLAST (II. 12) . I2«IRY2. IRV4) . Il»  PRINT752 

•  IRY1.IPY3)  PR  I NT2S3 

5  FORMAT (215. 3E1S. 6. •  < .3E1S.6.3*. 151  PRINT254 

Tflfl4  CONTINUF  PPINT255 

IFINPL.EQ.OI  GO  TO  920  PRINT2SG 

no  927  lal.NPL  PR  I NT257 

K  aLOCPOL ( I )  PPINT258 

WRITE (NWBITE.928) K.VROLE (X. 1 ) .YPOLE (K.p) .YPOLE (K.3) .OyPOLE (F.l) .  PMINT259 

•  OYPOLE («. 21 .OYPOLE (K. 31  PRINT2S0 

928  FORMAT ("OPOLE  ON  S IOE" .12."  HAS  POSITION  ANO  INCREMENTS  *"/  PBINT2S1 

•  101.3515.6."  "*  3E 15.0 1  PBINT267 

IF  (OIRCH)  GO  TO  92T  PR  I NT263 

WRITE ( NWRITE, 19281K.Y2P0LE (K.l I «  Y2P0LE (K.2) .Y2POLEIK.3) ,02POL£(K. 1PRINT264 

•  )  ,02P0LE(K,2)  ,02P0LE(K..l)  PRINT265 

1928  FORMAT!"  POLE  ON  SI0E".I2.".  r'i  POSITION  ANO  INCREMENTS  «"/  PRINT200 

•  10X.3E15.ft."  ".3E15.0I  PR  I NTPftT 

92T  CONTINUE  PR  I NT2ft8 

979  WRITE ( NWRI TE .930 )  PR  I NT269 

930  F0PMAT(/22X,»«  IPLAST  GREATER  THAN  ZERO  INDICATES  PLASTICITY  AT  TMPRINT270 

•IS  TIME  STEP")  PRINT2T1 

T08  CONTINUE  PR  I NT2T2 

C  PR  I NT2T3 

C  WRITE  LMATRIX  PBINT274 

II I»ITIME/IOUT (14)»I0UT(14)-ITIME  PBINT275 

IF ( I  1 1 . NE • 0  .ANO.  .NOT.QPRINT ( 14)  ) GOTO  T30  PRINT2T0 

WPITE(NWRITE.TU)  tihe.ITIME  PRINT2TT 

I7»0  PRINT2T9 

DO  T 1  0  ILAYER«1 .NLAYER  PPINT7T9 

NGAUSL*NGAUSS ( ILAYER)  PRINT280 

00  710  IGAUSS«1 .NGAUSL  PPINT281 

NS5L»NSUBL( ILAYER)  0RINT282 

00  710  ISBwi.NSBL  PRINT283 

IZ«IZ*1  PRINT284 

WRITE (NWRITE.715)  ILAYER . ILAYER, IGAUSS.ISi, ( 12 . I 2» IPY7 , IRY4)  PRINT285 

00  710  IlwIRYl.IRYS  PRINT286 

WRITE (NWRITE. 720)  1 1  , ( LMAT ( 1 1 . 1 2 , 1 7 ) . 1 7* IRY2 . I RY4 )  POINT287 

710  CONTINUE  PRINTERS 

711  FORMAT ( *  1  * , flX • 'TIME*' .El 3.A.8X, *  ITIMEw',15/  PRINT299 

•RX . 1 SURO IVISIONS  OF  TIME  INCREMENT  IN  STRESS'/)  PRINT290 

71 5  FORMAT (//20X, 'L"AT (1 1,12', 12,')  LAYFR',12,'  GAUSS  PT.',I2,  PRINT291 

••  SUBLAYER* . 12//'  II  12*' .40 13/)  PRINT797 

770  FORMAT ( •  '. 14, 5X, 40131  PR  I NT293 

C  WRITE  FMAT  PBINT704 

730  1 1 1*  I  TIME  / 1  OUT  ( 15 ) * IOUT  ( 1 5 ) - ITTMf  PRINTERS 

IFIIII.NE.O  .ANO.  .NOT. SPRINT!  15)  1G0T0  TftO  PRINTERS 

WRITE ( NWR  t  TF , 740 )  T I  ME .ITT  ME  PBINT297 

00  735  ILAYERal , NL  A YE°  PRINTERS 

NGAUSL*NGAUSS ( ILAYER I  PPINT799 

00  735  IGAUSS«l .NGAUSL  PBINT300 
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WRITE(NWRITE.T*5)  ILAYEP.  (GAUSS,  (12.  12* I STR2. 1 STRA )  PBINT301 

00  735  II*ISTP1,ISTR7  9BINT31? 

WRITE (NWBITE. 750)  tl . (FMAT ( I  GAUSS, n • 121 , 12* ISTP2. ISTOA)  PB I  ST  30 3 

735  CONTINUE  PBINT30* 

7*0  FORMAT (« 1 • ,8X. 'TIME*. .£1 3.6.8*. •  I T IM£* • .15/1  PPINT305 

7*5  FORMAT (//?0X, *FmaT ( 1 1 , l?, < , I?. ' 1  GAUSS  PT.'.I2.  PBINT306 

•//•  II  12** .4013}  P9IMT30? 

7*0  FORMAT!'  *. I*. OX. *0*31  PRINT308 

700  CONTINUE  PBINT309 

BETUBN  PBINT310 

PBINT3I 1 

STBESSES  PRINT31? 

PBINT313 

3000  IE  (ITINE.EO.OI  »ETljBN  PBINT31<. 

1 1 I*ITIM£/I0UT (*) »IOUT (*l -ITTNE  PB I  NT  3  IS 

IEIIII.NE. 0 . *N0 . .NOT , OPR  I NT ( * ) I  GO  TO  3003  PP I NT3 1 0 

lEIILAVER.EO.l.AND.IGAUSS.EQ.I  . *N0 .  ISU8L.E0.1)  *BI TE (NWR I TE . 032 1 pb I NT3 1 7 
932  FORMATI'M",  1*X."C0NTRAVABIANT  CONPONENTS  OE  TM£  STBESS  TENSOR'*//  “PPINT31* 
•  L»T£P  GAUSS  STATION  SUBLAYEP  T*USBLII.l)  TAUSBLII.2)  TPB I  NT  3 1 9 
•AUSRLtl.S)"/)  OPINT320 

WRITEINWBITE, 1933)  ILATEO. IGAUSS. ISURL  PBINT32V 

^933  F0RMATI2X.12.BX, 12, UX. I?)  PR1NT3?? 

WBtTElNWPTTE. 19341  !  (TAUSHL (I . J! *J*1 .3) ,1*1 .3)  PPINT323 

193*  EOPNAT (34X.3E15.0I  PBINT32* 

NSBL*NSUSL ( ILAYEP)  PPINTS’S 


IE ( ISUPl .E0.NS8U  WRITE (NwPITE.300?)  ILAYEP . 1  GAUSS . 

•  (  (TAU(I.J) .J*l .3) .1*1.3) 

3002  FORMAT!/**  C0NTPAVAR1ANT  COMP.  OE  TMF  STRESS  TENSOR.  TAU(I.J).  PPB 
•LAYER". 12.'*.  GAUSS  STA  T  .  ".  12.  "  «"/ ( *5*  .  3E1  S.  0)  I 

3003  RETURN 
C 

*000  IE  f ITTNE  .EQ.  0 ) RETURN 

lEdNOFXX  .EQ.  UIGOTO  8300 
III*ITINE/I0l'T(12)  •  IOUT  ( 12)  “IT  IN£ 

IEIIII.NE. 0  .AND.  . NOT .QPRINT (121 )GOT0  *010 
1 1* ( INOFX X-31 

GOTO (*100.5000 .6000. 7000 .800  0.81 00. *20  0  I . 1 1 

c 

*100  WRITE (NWBITE. *001 )  ISB 

*001  FORMAT ( *  SUBROUTINE  STRESS*/*  IS8»*.I3./I 
WRITE (NWBITE .*002)  ( (TNI J. I) .1*1 .31. 1.3) 

*002  FORMAT (•  TN(1.1)«*.E2?.15.'  TNI  1.2) *• .E?2. 15. >  TN ( 1 . 3) * • .E22. 15/ 

•  '  TN (2, 11 **.£??. 15, •  TN<2. ?)*'.£??. 15. •  TNI2.3I** .E22.15/ 

•  *  TNI3.1)**.E22.1S.*  TN(3.?)«*. £22.15. '  TN  1  3. 3 ) *  * . £2? . 1 5 ) 

*010  RETURN 

C 

5000  WPITEINWRITE. 50051  ( (TR(J.I) .I*1.3).J*1*3) 

5005  FORMATI*  TR( 1 . 1 ) *■ ,E22. 15. '  TR (1 ,21  *• ,E?2 . 1 5.  •  TR (1,31*'. £22.15/ 

•  '  TR(2, 1) *' ,£22.15. '  TR (2,2)** ,£??, 15. *  TP ( ?. 3) * ' ,E2? . 15/ 

•  •  TR 13. !)»*•£ 22. 15.*  TR 1 ?. 2) *  * . £22. 1 5. •  TR ( 3. 31 «* ,E22. 15) 

WRITE (NWRITE. 5006!  I  SB .L .L C .f 2 • S IGMSO 

5006  FORMATC  ••«  I  SB*'  .13,'  L»*.I3.'  LC»*.I3«*  CZ»* *E?2. 15 « *  SI&mSQ 
••.E22.15) 

return 

c 

600  0  WRITE  (NWRITE. 60 05 1  A2.R2.0ISCR 

6005  FORMATC  A7*'.E22.15.*  B  7*  ' .  E?2. 1 5  .  »  0  ISCP»*  .£  22 . 1 5 ) 

WRITE (NWRITE. 60 06)  ( ( TC ( J . 1 1 . I  *  1  *  3 ) . J* 1 • 3 ) 

6006  FORMATC  TC  ( 1 . 1 )  *  * .  £  22 . 1  5 .  '  TC  ( 1 .2)  «'  .£22. 15.  '  TC  ( 1 . 3)  *  •  .£22. 1 5/ 

•  •  TC(2.l)*'.E22,l5.*  TC (2.2) *• .£22. 15. •  TC ( 2 • 3 > ■' .£22. Is/ 

•  '  TC (3, 1 ) »• ,E22. 15. '  TC(3,21**.E22.15.»  TC ( 3 . 3 ) * » .£22 . 1 5 1 

RETURN 

C 

7000  WRITE (NWRITF. 700 1 )  HLAMOA 

7001  F0PMA7C  "LAMOA*  *  .E22. 15) 

WRITE (NWRITF ,7006!  ( ( TM ( J, I ) . I  *  1 . 3 ) . J* 1 • 3 1 

7006  FORMAT!'  TM ( 1 , 1 ) *  * ,£22 . 1 5 . •  TM ( 1 , 2  I  *  * . £22 . 1 5 . •  T“ (1 . 3 1  * • .£ 22 . 1 5/ 

•  *  TM(2, 1  I ,E22. 15. •  TM(?.2) »• .£2?. 15. '  TM < 2, 31  *  * . £22. 1 5/ 

•  '  TM(3, I )■*.£??. 15. •  TM(3,2) ** ,£22. 15. •  T«< 3. 3) * • .£2? . 1 5 1 

RETURN 

C 

*000  WRITF(NWRITE.*001>  lTTME.TlMF.IGAUSS.il . I? 

*001  FORMAT!'  I T I M£  a • , I  * . f  T I M£ ■ ' . £  IS , * , •  I  GAUSS* ' . I* • *  II**.  I**'  12*' 
•1*1 

wRITF (NWRITF .*002) 

*002  FOPMATC  SUBROUTINE  7FTAM 

WRITE (NWRITF .*0031  ( (G8ASE ( J . ( I . I* 1 . 3 1 . J*1 *3 ) 

*003  FQRM4T  (  '  GB4SF  1 1  . 1 )  •'  .E22.  15.  •  GBASE  1  1 .2)  ■  •  »E22. 15.  *  GBASEU.3I*' 


PRINT326 
OP1NT327 
POINTS?? 
PP  I N  T3  29 
BO1NT330 
OOINT331 
PP 1NT332 
PPINT333 
PP  T  NT 33* 
PPINT33S 
POINT336 
PBINT337 
PRINT33* 
POINT339 
BPINT340 
POINTS*! 
POI NT3*? 
pp ! NT3*3 
PP I NT3** 
PR INT3*5 
PPINT3*6 
PP 1  NT  3*  7 
PBI NT3»* 

do INT3*9 

PBINT350 
PBINT351 
■PP 1  NT  35? 
PD INT3S7 
oor.T35* 
PP I  NT  355 
PDINT356 
DPTNT3S7 
POINT353 
DDINT3S9 
DPINT360 
POTNT361 
PP I N  T  36  2 
PB I  NT  36  3 
pp! NT3S* 
POINT36S 
PP I N  T  366 
PR  I  NT  36  7 
PP I N  T  3  6  * 
PPTNT369 
PPINT370 
OP1NT171 
PO  T  VT  3  7  3 
,PQ I 

OO  J  ST  j  7  4 

OO  J  S.T1  TC, 
DU \ ST  3  TS 
DO  I S  T  3  7  7 

•  PW * S  T  3  7* 


6b 


•E22.is/  pptnt379 

••  GB4SE (2. I > ■• .E22. I*. •  GR4SE ( 2. 2) » • ,F22. IS. •  GR4SF ( 2 . 3 ) * • .£22. 1 5/BO  I NT350 
••  GB4SEO.ll**. £22. IS.*  GP4SE  (  3  .?>■*. E22.  1  5.  •  GB4SE (3. 3) *  * .E22 . 15 > PP ! NT 3B 1 
BRITE (NWRJTE. 80041  ( ( OG« 4SE (J. I ) . I » 1 . 3 ) . J* 1  *  3  I  BBINT38? 

5004  EOR«4T(*  0GR4SE (1 . 1 ) «• ,f 22. 15, »  DGP4SE 1 1 . 2 1 « • . £22 . 1 5 . •  06R4SE ( 1 . 3 ) Po I  NT  35 1 

•*••£22.15/  BBINT344 

••  0GP4S£(2.1 )**.E22.is.*  0G8  4SE  £  2  •  2)  •  • .E22.15. •  0GR4SE ( 2 . 3 ) * • • E22 . Bo  I NT3 55 
*15/  P01NT335 

••  0GR4SE  (  3. 1  )  »<  .£22. 15.  •  (1GR4SE  (  3.2)  *•  .£22.  15.  •  D6P4SE  (  3 . 3  )  *  '  .  F22  .BB  I  NT  3“  7 
•15)  BB 1  NT  34  5 

W9ITE(NWflITE.flOOS)  I (G( J. I ) . 1*1 . 3) • J*1 .31  BBINT35B 

5005  F0PN4TI*  G  II  .  1  )  *•  .£22.15.  •  G(1  .21**. £27. 15.*  G  I  1  .  3)  »<  .£22. 15/  PBINT390 

•  •  G (?. 1 ) * • .£22. 15, •  G 12.2) .£22. 15. •  G (2.3 1 • • .£ 22. 15/  PBINT391 

•  •  GO.  1  )  *•  .£22. 15.  •  GO. 2)**. £22. is..  G  (  3. 3)  «•  .£22. 1 51  PB1NT3P2 

WRITFINWPITE.BOO*)  GTVPf  bb I NT3P3 

800*  F0»M4T(*  GTTPE**  .£22.15)  PPINT394 

•BITE (NWSITE. BOOT)  <  ( GG < J • I  I  • I * 1 . 3 )  . J* 1 . 3 1  PPINT39S 

BOOT  £09M4T(*  GG  (  1  ,  1  )  »•  .£22. 15.  •  GG  1 1 . 2  )  »  •  .  E22  .  1  5  ,  >  GG  (  1  .  3  1  «  •  .  £22 . 1  5/  op  I  NT39S 

•  •  GG<2. 1 ) *• .£22. IS. •  GG ( 2  •  2 )  *  * . £22. IS  ,  •  GG ( 2. 3 ) * • .£22. 15 /  BBINT39T 

•  '  GGO.  1 1  ••  .£22. 15.  •  GG  0.2)  »•  .E22. 15.  •  GG  O.  3)  .£22. 15)  PRINT39B 

WRITE  (NWR 1 TF  .8008)  (  <0G4N(  J.  t  )  .  1*1  .31  .  J*t  .3)  PPINT399 

•  00*  E0PM4TI*  0G4MI 1  ,  1  )  ••  ,£22. 15.  •  0G4M  ( 1 . 2 )  *•  ,E22. 15.  •  PG4M  ( 1 , 3 1  « •  .E22PB I  NT400 
•.15/  PO1NT401 

••  0G4M(?. i ) .. .£22. 15. •  OGAM  12, 2) ** .£22. 18. •  0G4N ( 2, 3 1 « • - £22. 15/  PBINT402 

••  0G4M< 3. I ) *• .£22. 15. •  0G4MO.2)  *•  .£22. 15.  •  DG4M ( 3 . 3 ) * • .£22 . 15  I  PBINT403 

wQITF(NWfltTE.B009l  ( (OGAMMX ( J. It . !«1 .31 .J*l .3)  POINT404 

POOP  F0PM4T ( •  0G4MMX ( 1. 1) a* . £22.15, '  OG 4BBX (1 . 2 ) * • , £22 . 1 5 . •  0G4MMX ( 1 . 3 ) PP 1 N T* OS 
•*•.£22.15/  PBINT40G 

•  •  0G4MWX (2. 1 1  X. ,£22. 15. '  0G4MMX (2.2) *• .£22. 15. •  GG4m*x ( 2 , 3 ) * • .£22 .PP ! NT*0 T 

*15/  PBTNT408 

•  •  DG4MMX(3.l  >*..£22.15.*  0G4MMX  (3.2)  *•  .£22.15.  1  DGxmmx  ( 3 . 3  )  *  • .  £22. po  ]  NT4  09 


•15) 

PFTUPN 

P100  WRITE ( NWR 1TE.P101)  DG4MMX (3.31 

8101  FORMAT ( *  (REVISED)  OGAMMX0.3I«*.E??.1S) 

RETURN 

B20  0  WBITE(NWRITE.B201)  4G4H33 
B201  FORMAT!*  4G4H33*' .£22. 15) 

PCTIIBN 

MIXED  TENSOR  STRESSES 

B3  fl  0  WRITE (NWB1TE.B320)  T I wf , I  TIME . 1 1 . T 2 . IG4USS . T4USPH 
WBtTF ( NWS  I TE. 8325 1  ( (TAUE (K.J) .J*l .3) .*») .31 

PFTUPN 


PP r  NT4 1 0 
OPINT41 1 

PP 1 NT4 1 2 
PRINT413 
PBTNT414 
PRINT415 
PBINT415 
PPINT41T 
PBINT418 
PRINT419 
POINT420 
Pfl  1  N  T4  2  1 
PBINT42? 
OPINT423 
PBI N  T  4  2  4 
PBJNT425 
BPINT424 
PP I NT4  27 

PBINT42B 
PP JNT429 


B320  FORMAT ( *  time*' .E15.A.3X, < ITTwE  *'.I5.'  II*'. 13.’  I?*'. 13.  PPINT427 

••  G4USS  BT.  a*  « 1 3. •  T4USRM  *»«£l0.i)  PQINT420 

B325  F0OM4T  (58X  , 'MUEO  TENSOR  STRESSES*/  PPINT420 

•39X.I  T4UF ( 1 . 1 ) *• .Flo. I . •  T4UF ( 1 .2) •• .Fl 0. 1 . •  T4UF (1.3)**, £10.1/  PPINT430 
•3RX  «  *  T4UF  (2. 1  )  »•  .Fl  0.1 .  •  T4iiF(2,2)*«,£10.1.  •  T  4(Jf  (2. 31**. £10.1/  PQINTOI 
•  IPX.*  T4U£ ( 3. 1 ) »• .£10. 1 . '  T*U£ (3.2) *• .£10. 1 . '  T4u£ ( 3 . 3 1  * ' . £  1 0 . I  I  0907432 
ENO  PBINT433 


SUBROUTINE  N  (  T  AU  .1  CO 

IWHICI!  LOGICAL  ui 


CU^ON  / INOE  */  NRfc  Aj  ,.x  #8  [  TE.NPuNCM.NNtShl.N^ES^.M^.N^.MONM. 

•  si"*',h^n.n.i^iz.iii£9Cii^f»Qii(miNr^t  lariJKY^iSTki, 

•  lSTQ*,ISTR«»«Kl.ZL^lCl,KW*lr>l*l^*i5l*IW«*l*'<<?fK3,A4’«*;<LS. 

•  k/stor  .« yte s r .  iui*.  jitest,  i?t£St»kim  tl 

c  3**0*  /  r**NS/  ; ,  i  i  •  2C£s7R»TniC<N*A0Stisi6.M.«fiu*ri**M. 

•  sG#USS<^*»ICAU»S.NLAYt*».UAYEk*sSu8H^>*ISUBL 

CO*"CN  /**  T  !  *E  /  aua ( *  C J  «  T  I  *k  •GELTAT.TlNtf.l  T  ,  II  i  "ff  ,  l  Aox ( 20 )  * 

•  ;cLrcjo«-*>‘MMttv» 

CG**Cn  /C  T  I  -t-  l  /  IRLASM  2C.2C  I  .2 I  2U.20  )  «?Pl  i  20,  2C  ) 

L  6  V  E  L  2.  I  PL  AS  T  .  P  .P^u 

J I -E  N$  J  Cn  T  au  <  3  ,  il  ,  Tau-a<  U  ,  )j  .TauONU,  3i  ,  i  MAX(  ]»3l»li“l*»<3t3)t 

•  I27'A*«3.j).i2MN(J.)),IG',PM3.3>.lG“IMJ,Jl.iTIMJ.j)  *  M  •-  (  3.31 

•  .  I  T  I-PI  3.  3  I  .T  MM  J  *  3) 

OAIA  CH/.E  ALSfc  .  /  MAX /•  “AM/MIN/*  "  t  N  *  / 

G0T0(b*35». IGJ 

S  (MCN)GQTC  i'O 
CO  1C  I  - l . 3 
00  10  J-  l  .  3 
T  AUf*  AX  (  I  .  J  )  *0*0 
TAUMNl  l  .  J)-O.C 
I  MAY(  I  .  j  1  -0 

1  IMN«  1  *J>"C 

I2*AX<  I  .  J  )  -C 
I2MN4  1  ,  J  I  ■  0 
IGPfNC I • J  >*0 

10  CONTINLE 
QW - . TOlE  . 

20  DU  30  1-1.3 
no  30  J-1.3 

01  l-TAU(  I  .  J)  -GT.  TAU^AX(I,J| 

IP (  •  NQ T  .  Cl  l  JOGTU  2*> 

TAUPAX  (  1  .  j  )  ■  T  A  O  (  I  ,  J  I 
I 1NAX ( I . J ) • l l 
UMAX  <  I  »  J  )  »  I  2 
I  T  Ml  I  .  J  I  -  I  T  I-fc 
TMII.  JI-TI-6 
( G N A X ( I  .  JI-  ICALSS 
2*>  Cll-TAL(I.J)  .Li.  TAUfllMl.JI 
I F  < .NOT.  U  1 1  )Gu  T  C  30 
TAUFINI I . J) «TAL i i , J ) 

UMMliJIMi 
t2"  (NU  «JI  *  12 
ITI“NU*J)-ITME 
ri^(I,JI*T|."E 
I GN  I N  4  I  .  J  )  -  I  GAoS  S 
3u  CONTINUE 
BE  TUPN 


-A  X" I N  1 
* AX  -  I N  2 
1  A  »  p  [  N  3 
MAX^IV  <• 

P  A  *  *  I  N  *> 
PA  *P  I N  b 
-A  *  -  I  N  ? 

«  A  x* I N  4 
P  A  *  P  I  N  ? 

-  A  i  P  I  N  l  0 
P  A  X  »  I  N  l  1 

-  A  *  P  i  N  1  .? 
p  A  *  P  l  N  l  3 

-  A  X  P  I  N  1  P 
**  A  X**  i  N  l  5 
-A  XP  1  N  It 
-A  XM  I N  1  7 

-  A  XP I N  t  d 

-  A  X  P  I  N  M 
P  A  X  P  I  N  2  3 
- AXP  I N2  1 
- axp 1 N22 
-A  x* I N2  J 

-  A  <P  [ N2  N 

-  AXP  I  N  2  *> 
P A  X  *  IN  2  b 

p  Ax- I N2  7 
p  A  x  p  I  s  2  s 
p  A  *  *  I  N  2  ^ 

-  A  x-  1  N  30 
-AX-  1 N  31 
PAX- I N 32 
p  t x“ I N  3  3 
PAX-  I N  3<* 
p  A  X  p  I  N  3  ** 
-A  X  -  I N  3b 

-  A  x  -  I  N  3  7 

P  A  X  P  I  Nib 

-A  *  -  I \  3  * 

-  A  x“  i  NP'> 
p  A  X  p  I  S  *•  l 
p  a  X  -  i  \  P  2 

PAX  “  |  N“  3 
P  A  t  p  l  N  <.* 

•  A  y  “  1  N*  *> 
-Ax-  (  \*b 
p  A  y  P  I  N  H  ? 

p  A  X  p  I  N  *»  P 

-  A  X  p  I  N  «.  0 
PA  r-  I  NSif 
p  A  x  P  ;  N  p  1 
“  A  x  p  I  N  *>  2 

p  A  X  -  I  N  *3  ) 

-  A  X  *  1  N  S  t> 
p  A  *  p  l  S  N  b 

-  A  X  -  I  N  *>  b 

*  4  »  “  I  S  •)  ." 
p  A  X  -  [  N  S  •% 


j  *?  rf*  1  TEl  N-Rl  TF  .1 JC)  pax 
NO  DU  pS  I  -  I  .  3 
r*ii  x,  r,  » «  \  .  x 

tad  I  TE  (  NhH  ITE  *1 10)  i  •  J  I 

»  •  I  C»4*  (  l  ,  i  1  ,  l  ,  J  ,  r  4U»<4  l  <  I  .  J  > 

contincf 

taRITEINtadirE.lJOl  1IN 
00  50  1-1,3 
DO  50  J-1.3 
tad  I  TE  <NtaR  I  TE  ,  L  10  3 
•  .  I  C»IM  1  .  J  )  ,  I  .  J  ,T4U“  I  ,N<  [  ,  J  I 
50  CONIINCE 

»e  turn 

IOC  F0dta4Tt/'  •,4‘.,<l”Ud  “!*E0  tEtaSOd  SIdE55E5'/l 
l  10  FORMAT!1  TI-t-'.Ell.ft.'  I  U"E  -  *  .  I  5  .  *  f  1  -  *  -  1  ^  *  U  -  *  .  !  /  . 
•'  C4USS  EAUF  I  ’  .  i  l  ,  •  .  ’  ,  U  ,  •  t  ■*  ,E  1 1 .5  I 


“A » -  I \fi  ; 
p  A  x  P  !  *4  t  I 

«*  4  »  -  ’  *:  * 

P  A  X  M  I  Nfe  3 
p  A  X  *  I  V  b  P 
“  A  X  p  |  S  t>  •> 
pax  -  1 Nbb 
pax-is^; 

-  &  X  -  I Sbrt 
" A  X  -  INN) 

-  Ax  -  I  N  7D 

-  A  x  *  |  N  7  l 
p  a  X  p  i  n  7  2 
•  A  X  p  I  N  7  3 

-  A  x  -  I  N  7  U 

-  A  X- 1 N  7S 
“ A  »  p  l  N  7b 

-  A  x  p  i  N  7  7 


c 


c 


c 


c 


c 


c 


c 


e 


] 


?* 


10 

i* 


is 


SUBHOOT T  NE  SFORCE  (CS.TMtCKN.YY.YY2. I  1 . I2.SURFGG) 

SURF  ACE  FORCES  (ONLY  FOR  SHEAR  OPTION) 

COMNON  /CARTE/  YTEST.YNEN.YSAVF 

COMMON  /cartel/  Y (3.20.201 .0 (3.20 .20) .Y2 (3.20. 20) .0? (1.20.20  I 
LEVEL  2.Y.0.Y2.02 

COMMON  /PUSH/  FORCES  (1)  .VELOCO)  . RAT  10. RATION, OX  1 .0*2. TEMP, OTEMP , 

fspace.tspace.fincno.fstop.tstop.thcoef 

common  /PUSHL/  SO»AT<20.20) .SORAZI20.20) .FMAS1  1 (20.20) . 

FMAS22(?0.20) .FMAS 23(20.20) .FMAS33 f 20 .20) 

LEVEL  2.S0PAT.SQRA2,FMASn  .FM4S22iFMAS23.FMAS33 

common  /SUBNOM/  SNPR(1) 
common  /SUHNOL/  SNO.20.20) 

LEVEL  2.SN 

COMMON  /SJPFOR/  El (3.20.20) •F2I3.20.20) 

LEVEL  2. El .E2 

01  MENS  I  ON  02 (3) .CSI3.P) .58  ASF (3.3) .5 (3. 3) .GG ( 3. 3) • YY ( 3. 2) . VY? ( 3. 2 ) 
•SUMEI3) 


SFOOCt  1 

sforce  R 

SFCOCE  3 
SFORCE  » 
SFORCF 
SFORCE 
SFORCE 
SFO»CE 
SFORCE 
SFORCE 10 
SFOcCE 1 1 
SFORCE l 2 
SFCaCE 1 3 
SF0RCE1* 
SFO°CE 1 5 
SF0RCE1S 
SFORCE  1  T 
SF0BCE1S 
SF0RCE1R 
SF0RCE20 
SF0RCE21 
SFORCE’? 


00  3  J*  l  «  3 
El (J. II. 121*0.0 
E2IJ.I1. 121*0.0 
CONTINUE 

00  SO  122*1 .2 
22*0 

CALL  FDASE  (02.3) 

02(31*1.0 

GOTO  (2S.30I.I2Z 

22*-.5*TmJCFN 

S I GN* 1.0 

CALL  lOAOL 

GOTO  3* 

22».S*TMICKN 
SIGN.-i .0 
CALL  LOAO 
00  3S  J*  l  .  3 

GRASE (l.J)*YYij,l) -Z2*CS ( J. 1) 

GRASE (2. J) *YY  < J.2) -Z7‘CS  < J.2) 

GRASE  (  3.  J)  *SN  ( J.  11,12) 

GRASE ( 1 , J ) *GRASE ( l . J) .02 ( 1 ) »Y?( J. II . 12) »22»yy? [ j, 1 | 
GRASE  (  2.  J I  *GBASE  (2.JI  «02(2I*Y2  < J. II . 12) ♦ 2Z‘VY2 < j.2 ) 
GRASE ( 3 . J ) *GH ASE I  3. Jl *02 ( 3) *Y? ( J. II ,12) 

CONTINUE 

lill  •  j  )  *GBASE  (1.11  *GBASE  ( 1  .  1 )  .GBASFI  1,2)  *GRASE  (1.2) 

’  .GRASE ( 1 .3) *G8ASE (1 .3) 

G( 1 .2) *G8ASE (1.11 ‘GRASE (2.1 ) .GBASE ( 1.2) *GR ASE (2.2) 

'  .GBASE ( 1 .3) *GBASE (2.3) 

G( 1 . 3) *G8ASE (1,1) *GRASE (3.1) .GRASE 1 1 .2) *GRASE (3.2) 

'  .GRASE ( 1 . 3) ‘GBASE (3.3) 

6(2. 1 ) *G< 1 .2) 

G(2.2) »G8ASE (2. 1 )*GRASF (2. 1 ) .GBASE (2.2) .GRASE (2.2) 

'  .GRASE ( 2. 3) * GRASS (2.3) 

G( 2. 31 *G8ASE (2.1) ‘GRASE ( 3. 1) .GBASE <2. 2) .GRASE (3.2) 

'  .GBASE ( 2.3) ‘GBASE ( 3. 1) 

6*3.1) *8 (I. 3) 

5(3.21*6(2.3) 

5(3.3) *08 ASF (l.li .GBASE (3.1) .GRASF (3.2) ‘GBASE (3.2) 


SFORCE  23 
SF0RCE2* 
SFORCE  ’5 
SFORCE?*. 
SF0RCE2T 
SF0RCE28 
SFORCE ’R 
SF0RCE30 
SF0RCE31 
SFORCE  3? 
SF0RCE33 
SF0RCE3* 
SFC-0CE35 
SFO&CE  3*. 
SEC  RCE  3  7 
SF0RCE3B 
SFORCE  33 
SFORCEAO 
SF0RCFA1 
SFORCE  A? 
SFORCF  (.3 
SFORCE** 
SFODCE*S 
SFORCE  *0 
SFC=CE*T 
SFORCE  *8 
SFORCE  *  R 
SF0=CE50 
SF0RCES1 
SFORCES? 
SFCRCFS3 
SFORCES* 
SF0RCF55 
SFORCES* 
SFGBCES7 
SFORCES? 
SFORCES* 
SFDRCESO 
SFO’CESl 
SFO=CES? 


71 


•  .GBASE ( 3.3) •G84SE ( 3. 11 
G.GTYPE 

GTYPE.G  ( 1  . 1 |4<G(?.2I •G(3.1)-Gi2.3)  ••?)  •6(1.JI»(GIU3I*5IMI- 

•  C,  ( 1  ,2)  #G(  1.1)  I  .G(  1 . 3)  »  (G  (  1 . 2)  »G  (2.3) -6  ( 1  .  31  «G  (2 .2)  > 

SQRG.SQRT (GTVPE) 

GTYPE*1 ./GTVPE 

GG(i.ll«(G(2.2l»6(1.1l-G(2.3l»*2)»GTYPE 
GG ( 1 .21. <G (] .3) *6(2. 1) -G ( 1 .2) »G<3.3> >»GTYPE 
GG ( ?. t  >  *GG (1.2) 

GG ( 1 .3) «<G( 1 .2) •0(2.11 -6(1 .3) *6(2.2) ) *GTYPE 
GG ( 3. 1 ) «GG (1.3) 

Gtt (>.?)■ (6 I) .1  I *613. 31 -6 11 .3) **2) *GTYP2 
f.G  (2. 3)  »  (G  ( 1 .21  »G  <  1 . 31 -G  (  1 . 1 )  »G  (  3.2)  )  *GTVP£ 

GG ( 1. 2) *GG (2.3) 

G6t 3. 3) «(G( 1. I > *G( 2.2) -G ( 1 .2) ••2)»GTYPE 
IP ( I  Z?  .20.  2)  SUR2GG.GG<3.3) 

CALL  ERASE  (SUME.3) 

GO  40  K«1 .3 
00  40  1.1.3 

SUME (K).SUHE (K ) .GG (3. I l«GBASE (  I  •  K  ) 

40  CONTINUE 

00  45  J.  1 . 3 

12 (SN( J. 1 1 , 121  .EG.  0.0) GOTO  45 
PRESS*- (FORCES (JI/SN ( J.I 1 . 12) ) 

TENP»S0RG4SUNE I J) 4PR2SS 
El ( J. II . 12). El ( J. II . I2).SIGN«TEMP 
E2( J.I 1.12 1* E2IJ.I1. 1 2) .SIGN.7Z4TEMP 
45  CONTINUE 
50  CONTINUE 
PETURN 
ENO 


SFC9CES3 
S20BCE44 
SF09CE55 
SF0OCE45 
S202CE5T 
SF0OCE55 
SE0BCE5R 
SF0OCET0 
SF0BCET1 
SF0BCET2 
SF0SCET3 
SEC  BCE  T  4 
SFOSCETS 
SF0PCE75 
SFOBCETT 
SF0BCET5 
SF0BCET9 
SEOBCEAO 
SF0BCE41 
SF0»CE9» 
SF0BCE43 
SF0BCE54 
SFOBCESS 
SF0OCE56 

SF0PCE3T 

SFOBCEBP 

SF0BCE49 

SF0BCE90 

S2DBCEU 

S209CE92 

S200CF93 


SUBROUTINE  STRESS  (  JO  1  ,  J02 , JC3 . TAu 1 1 . r AU 12 , TAU 1 3 .  STRESS  1 

•  TAU21,TAU22,TA023,rAu31,TAU32,TAu331  STRESS  2 

STRESS  3 

EVALUATE  STRESS  INCREMENTS  AND  STRESSES  STRESS  <1 

stress  s 

IMPLICIT  LOCICAL1CI  STRESS  6 

stress  7 

COMMON/  ALLENE  /TOTAL. TOTKIN.TCTELA.TQTPLA.TOTkEX.TOTTEM.INERGY  STRESS  6 

•  .  TOTV IS.TOTE I .  TC  T  E2  .ECHECK  STRESS  9 

•  .OTMI3.3l.SORC.SCRA.SlCMSa.A2.32.C2  STRESSIO 

STRESSll 

COMMON  /C  T I  ME /  AUXl 20)  ,T IME .OELTAT.T IMEF . I T I  ME . IT IMEF . I A0X1  20) .  STRESS12 

•  I0UTI20).CRRINT(20)  STRESSES 

COMMON  /CTIMEL/  IPLAST  (20.20  .21  20.20)  .RRLI20.20)  STRESSl* 

LEVEL  2. IPLAST. P./RL  STRESSES 

STRESSE6 

COMMON/ C TIMER/ IT IMEC.ITIMER.CELTAP.OE LX. OHR.UNh, HE E  STRESS17 

•  .  TKEEP  .MTHIK,  Cf  IMS  ,0F  INP.  TSTaRT  ,r  STAR  T  .  rOOTF  STRESS18 

•  .ES.8STIV(*),nSTIV(*)  STRESSIO 

STRESS20 

COMMON  /  I  NO  EX/  NREAO.NI.R  I  TE  .  NPUNCH  ,  NME  SHI  .NMESH2.NE.N2  .N2.MM.N2M.STRESS21 

•  NEMM.N2MM.I1.I2.12.I12ER0.122ERO.IRYL.IRT2.IRT3.IRTN.ISTRE.ISTR2.STRESS22 

•  ISTR3.ISTRN*1CE.1C2«1CE.1C2.1P1.IP2.ISE.IS2.RE.R2.R3,KN,kRcn,  STRESS23 

•  K2ST0P,KYT£ST,I0IR.UTEST.12TEST,<INITL  STRESS2N 

STRESS2S 

COMMON  /OPTION/  mauXIL.MINCEO.MINVEl.MLOAD.hmaTRR.MSPlCA,  STRESS26 

•  MSPTEM.MTEMPE.MTHIKL.PIMPUL.1STRES.INCRML.ISTRE2  STRESS27 

COMMON  /CPTFRA/  IFRACT.CFRACT  STRESS28 

COMMON  /MAI  I /  NCCNT,NRITE,NTRAIN.MRIT£,NDELP.ETA01.ETAC2.NSTRN  STRE'S2<) 

«  .FACTOM.MO.FACTON.NO.FK IN  STRE. ,30 

STRESS3E 

COMMON  /PHYSCN/  EE ,HNU, alpha .CONST. EXPCN, Factor .RATE ,RHC,  STRESS32 

•  HLAMDA.COEFF ( 5 ) ,S I CMA l S ) , Tr ( 3, 3  I , TC I  3, 3  I , DELTA ( 3. 3 )  STRESS33 

STRE  SS3* 

COMMON  /PUSH/  FORCES! 3 )  .VELCCI 3  ). RAT  10, RATI OM.CXE ,0X2 . TEMP .CTEMP,  STRESS3S 

•  FSPACE .TSPACE ,F INCNO.FSTOP .TSTCP. THCOEF  STRESS36 

COMMON  / P  u  S  HL /  SCRATI20.20)  .S0RA2 120,20  ) , FM ASl 1 1  20 ,20 )  ,  STRESS37 

«  FMAS22(20.20),FmaS23(20.2CI.FMAS33I20.20)  STRESS38 

LEVEL  2.SCRAT.SCRA2 .FMAS11.FMAS22.FMAS23.FMAS33  STRESS3N 


COMMON  /OPTION/  MAUXIL.MINCEO.MINVEL.MLOAD.MMATRR.MSPLCA, 

♦  MSPTEM.MTEMPE.MTHIKL.PIMPUL.ISTRES.INCRML.ISTRE2 
COMMON  /CPTFRA/  IFRACT.CFRACT 

COMMON  /MAI  I /  NCCNT.NRITE.NTRAIN.MRITE.NDELP.ETA01.ETAC2.NSTRN 

♦  .FACTOM.MO.FACTON.NO.FK IN 

COMMON  /PHYSCN/  EE ,HNU, ALPHA, CONST. EXPCN, FACTOR .RATE ,RHC, 

♦  HLAMDA,COEFF!5),SIGMA!5),Tr(3,3),TCt3,3l,DELTA(3,3) 


•  FSPACE. TSPACE. FINCNO.FSTOP. TSTCP. THCOEF  STRESS36 

COMMON  / P  u  S  HL /  SCRATI20.20)  .S0RA2 120,20  ) , FM ASl 1 1  20 ,20 )  ,  STRESS37 

«  FMAS22120,20),FmaS23120.2CI.FMAS33I20.20)  STRESS38 

LEVEL  2.SCRAT.SCRA2 .FMAS11.FMAS22.FMAS23.FMAS33  STRESS3N 

STPESS*0 

COMMON  /CLOCIC/  0AUX12C1.C2ETA.CSTRES.CPLAST.WSENS1.CECUIL.  STRESS*) 

•  OOIACN.OINCEO.OINVEL.CLOAC.OMATPR.OTHIKL.OTEMPE.OSPTEM.CAUXll,  STRESS* 2 

•  OAUX12,CSPLOA,CIMPUL.CSHARP,CP£SO,CIRCH.OSHEAR  STRESSN3 

STRESS** 

COMMON  /TENCOM/  YYI3. 2). YYY13, 2.2). A13.3). o(3«3). AA13, 3). 6613.3),  STRESS'S 

•  8M(3,3).0AI3.3).C8(3.3).C(3.3),CG(3,3).0M3).CC13,21.CCC(3,2.2).  STRESS** 

•  DC  AM(  3,3  )  .OGAM  MX  1  3,31  ,HM3*3),HO(2).TAU13.3)*TAuS0L(3.3)  STRESS*7 

•  ,0021 3,2) .00021 3,2,21 ,YY 21 3,2), YYY 21 3,2.2) .yycI3,2),08mi3,3)  STRESS*8 


•  0GAM(3,3).0GAMMX(3,3I  .HM3,3),HO(2).TAUI3,3).TAuSBH3.3)  STRESS*7 

•  ,0021 3,2) .00021 3,2,21 ,YY 21 3, 2), YYY 21 3,2.2) ,yycI3,2),06m13,3)  STRESS*8 

STRESS** 

COMMON  / THK  NS /  2 . 22 .2CENTR . THICKN. A8SC IS  I  6.6 1 . »E IGHT (6.6 ) .  STRESSSO 

•  NCAUSS 1 *) . IGACSS.NLAYER, ILAYER .NSU8L l *  1 . ISUBL  STRESSSl 

COMMON  /CUt/  OGAM  33  STRE5SS2 

COMMON  /GG/  SURFCC  STRESSS3 

ST»E  SSS* 

COMMON  /OELC/  I C  CUNT  STRESSSS 

STRE  SSS6 

COMMON  /STROUT/  ISB.L.LC  .OISCR.TRl 3. 3). TNI  3.3)  STRESSS7 

STRE  S  SSd 

COMMON  /FRAC/  TAUF  13.3  I.TAUSPH.NUM.  IPS  1 1 10)  .  IPS21 10  1  ST»£SSS<) 

STRE  S  S  60 

COMMON  /STRPLO/  TAUP16C.3.3).  TAUSPL11C.6)  STRESSSl 

COMMON  /MATRIX/  YLOFAC.ANUM  STRESS62 

COMMON  /MATRIL/  LMAT120.20. It) .FMAT 16.20.20  I  STRE5S63 

LEVEL  2.LMAT.FMAT  STRESS6* 

STRESS'S 
STRE  SS66 

01  ME  NS  ION  TAUlll J01  .  J02. J03 I ,T AU121 JDl, JC2. J03I .TAU 1 31  JCI. JC2. JO  3 1 STRE SS6  7 

•  , TAU2 11  JO) . J02 .J03I.TAU22IJC1.302. JO  3) .TAU231 JCI . J02.JC3)  STRESSES 

•  . TAU3 11 JOI«J02,J03).TAU32(jCL,J02,JD3).Tau33(JC1.JD2.JD3)  STRESS  60 


IF  (OSTRES)  GO  TO  20 
LEN-JOl* J02»JD3 
CALL  ERASE1TAU11.LENI 
CALL  ERASE1TAU12.LENI 
CALL  ERASE  I TAU13.LENI 
CALL  ERASEI TAU21 .LEN) 
CALL  ERASE1TAU22.LENI 
CALL  ERASEITAU23.LEN) 
CALL  ERASEI TAU31.LENI 


STRESSES 
STRESSES 
STRESS70 
5  T  RE  5  S  7 1 
5TRESS72 
STRE  S5  7  3 
STRE  SS  7* 
STRESS7S 
STRESS76 
STRESS77 
STRE  S578 
STRE  S  S  76 
STRESSSO 


CALL  ERASEITAU)2,LEM 
CALL  ERASEI TAU13.LEN) 

CSTRES-.TRUE. 

C 

20  CONTINUE 

0002*. FALSE. 

0003*. FALSE. 

0004*. FALSE . 

IF  I  ISTRE2.EQ.2)  0002*. TRUE. 

IF  ( ISTRE2.E0.2.CR. ISTRE2.E0.3I  0QQ3*.TRuE. 

IFIISTRE2  .EC.  4 t 0004- . TRUE . 

C 

OGTE«P*-U.  »MNU  I  PALPHAPCTEPP/HNU 

CGAN*A-CCA)INXt  l ,  l 1  *0  CAP  NX  l  2 ,2  I pOGAMPX  I  3  ,3)»OCTE"P*3. 

1 F ( 0003  .OR.  0CC4I  aGAPPA-CCANPXI 1 , l )*CCAMNX(2 .2  I •UGTE PP»3 . 
IFIHNU  .EO.  0.51  CCANNA*3.4(0GTEPP*ALPMA*0TEPP 1 
IFIHNU  .£0.  0.51  HNUP-HNL73. 

C 

ES-EE 

IF  IBSTIVI ILAYER) .NE.O. I  CALL  ESTIFF 
FACTOR*!. 

EEP*ES/Cl.*HNUI 
HNUPP*  L . / 3. 

IFIHNU  . NE •  0.5)  HNUP-HNU/ I  1.-2 . PHNU) 

IF  I  COO  3  .OR.  0C04I  HNUP-HNU/ I l.-HNU) 

IFIC003  .OR.  UCC4  1  HNUPP" II .-2.4HNU)/ I  3.* ( l.-HNU I  I 
C 

IF  I  CONST. NE .0. 1  CALL  SENS1T 

C 

CALL  ERASE1TAU.9) 

CGAN33-0.0 
NSBL-NSUBLI  ILATER) 

C 

c 

Call  erase ioth i 
Call  erase  itauf,9i 

C  THE  FCLLOxING  T33.T33PL  CALCULATION  IS  CNly  valID  FOR  INITIALLY 

C  CONSTANT  THICKNESS  PLATE 

C 

IF!. NOT.  C004ICCTO  25 

T33— PI  I1,I21*SUi>FGG*I  L  .  *  2  .  ♦// TH  ICKM  /  2  . 

T13PL  — PPLI  1 1.  12  l*SURFCCP|  1.p2.p2/THICKM/2. 
OTAU33-G(3,3)*IT3)PL-r13l 
25  CONTINUE 

IF!  ICOUNT  .EO.  C I  GO  TO  tCIO 
IFII1.E0.IS1  .ANC.  12. £0*1521  GOTO  1000 
GOTO  1C10 

1000  I  I  I -IT  I  NET  I  OUT  I  121 • I  OUT  I  12  I- I  TINE 

I F  I  I  1 1  .NE.  0  .AND.  ..NOT. OPR  INTI  121  1GOTO  1010 
NR  I  TE  I  <1,100  7  1  Pill, 121, PPL!  1 1 ,  1  2  I  ,  SURF  GG 
1007  FORHATI*  PI  1 1, 12  I -' ,E22. 15, '  PPL!  1 1 ,  1 2 1  *  *  ,  E 22 . 1 5 . 

••  SURFACE  GGI3,31**,E22.15) 

1010  CONTINUE 

C  OETERNINE  STRESSES  FOR  EACH  SUBLAYER 

00  250  ISB-l.NSBL 
0C2*. FALSE. 

LUI2-0 
L-l 
EL  «  l  . 

SICPSO-SIGNAI I  SB  IPS IGHAI ISP  1 »F ACTOR »F AC  TOR 
I  2  •  I  2  ♦  1 

C 

CALL  ERASE  I  TN , 5  1 
TNI  1,11-TAUllI  I  1  .  12,  12  I 
TNI  1 ,2 1  * TAu 12 1  I  1. 12 ,  121 
TNI  2.  I  1-TAU21I  I  1 ,  12, 12  I 
TNI2.2)*TAU22! 11,12,121 
IFIOCO21COT0  3050 
TNI 1.31-TAU13I II , 12  ,  12  1 
TNI2.3I-TAU23I I l, 12, 12  I 
TNI  3, II-TAU31 I  I  1  ,  12,  12  1 
TNI  3,2 1-TAU32I 11,12.12) 

IF  I CC0  3  1GOTO  3C50 

TNI  3, JI.TAU33I 1 1.12 , 12  1 

IF1ITIPE.GT.1  .AND.  TNI  3,3) . E 0. 0 . 0 1 T 3 3*0 . 0 
I  F I C0Q4 ) TNI  3  ,  3  I  *  T  3  3 
3050  CONTINUE 

c 

C  CALCULATE  -  TMI.JI  AT  EACH  SUBLAYER  ANO  SU* 


STRE  SS8  l 
STRE  SS82 
STPESS83 
STRESSB4 
STRESS85 
STRESS86 
STRESS97 
STRE  SS88 
STRESS89 
STRESS90 
STRESS91 
STRE  SS92 
STRE  SS93 
STRE  S  SOL 
ST  PE  S  S  95 
STRESS  36 
STRE  SS97 
STRESS')!) 
STRE  SS99 
STRE  S 100 
STRE  S101 
STRE  S 102 
STRES103 
STRES104 
STRES105 
STRE  S106 
ST»E  S 107 
STRES108 
STRES109 
STRE  S  1 10 
STREStll 
STRE  S 1 12 
STRESU3 
STRE  S  114 
STRES115 
STRE  SI  16 
STRE  S 1 1 7 
STRE  S  1 18 
STRE  S l 19 
STRE  S 120 
STRES121 
STRES122 
STRES173 
STRES124 
STRE3125 
S  TRE  S 126 
STRE  S 12  7 
STRES128 
STRE  S 129 
STRES130 
STRE  S l 3 l 
STRE  S132 
STRE  SI  33 
S  TRE  S 1 34 
S  TRE  S 1 35 
STRES136 
STRES137 
S  TRE  S  l  38 
STRE  S 1 39 
S  TRE  S 1 40 
STRES141 
S  TRE  S 1 42 
STRES143 
STRES144 
STRES145 
STRE  S  146 
STRES  147 
STRE  S  148 
STRE  S149 
STRES150 
STRESISl 
STRFS152 
STRES153 
STRES154 
STRES155 
STRE  S IKfc 
ST8ES1S7 
STRE  S158 

ST»E  S150 


74 


CTNIl.ll- 
0TMI.2I- 
CTNI2.1I- 
0TM2.il- 
IF  ICOC2I 
0TM1.3). 
0TM2.3I- 
CTKI3.1  I- 
QTM3.21- 
0  TM  3*  3  1- 
50 50  CONTI  NoE 


'OTNll, 
‘CTNI 1, 
■0TM2  . 
■0TNI2. 

CO  TO 
l0T  Nil, 
'OTM2  . 
0T*<3. 
0  T  f-  (  3  , 
'CTM3. 


II  — TNI  1 
2 l-TNI  l 
1I-TNI2 
21-TNI2 
5050 
31-TNI 1 
3  I  -  T  N  4  2 
1I-TN13 
2I-IM3 
3 l-TNI 3 


«  1 1  -CCEFF I  l  Sd  I 
,  2 1 -CCEF  F I  I  SB  I 
-1MCCEFFI  ISB) 
,2 l-CCEFF I l SB  I 

.  3  I  -CGEFF I ISB I 
,31-CCEEFI ISB) 
. I l-CCEFF I  I  SB  I 
,21-CCEFFI ISB) 
.  31-CCEFF  (  I  So  I 


Tn(la.l8I  ARE  THE  NIxEC  INITIAL  STRESSES  FOR  IMS  PARTICULAR 
TER 

SU8CIV1CE0  TINE  STEP  STARTS  MERE.  IF  -ECUIREO 
55  CONTINUE 

TNI  1 ,1 l-TNI 1 .1  I 

TNI  1 ,2). INI  1.2  I 

TN II, 31- TN4I.3I 

TNI2.ll-TNI2.il 

TNI2.2I-TM2.2I 

TNI2.3I-TM2.3I 

TNl3.lI-TNI3.il 

TNI  3,2 l-TNI 3,2  I 

TNI3.3I-TM3.3I 

TAUN.TM  l.l  I  -  TNI2.2  I  -  TN|  3,3  I 

LC-l 

IF  I  I  COL  N  T  .E3.  II  CALL  X  L  X  I L I  5 | 

CALCULATE  “IXEC  TkIAL  ELASTIC  STRESSES  TRILA.L3I 

50  CONTINUE 
ML  A«OA-C . 

CALL  ERASEITR.91 

ETERNl-EEP/EL 

ETERN2-hnuP«0GAN-a»ETERNI 

TR  I  1.1  I -TN|  l ,1 | -ETERNl-OCANNXl 1, l ) -ETERN2 
TRIl,2l-TNIl,2l-ETE8Nl*0CA«Nxli,2l 
TRI2.1I-TNI2, ll-ETE8Nl-CGANNXl2.il 
TR I  2.2  I -TNI  2, 2  I ♦ETERNI»CCANNX12,2) -ETERN2 
IF  I00C2I  CO  TC  6050 

TRI1.3)-TNI1,3)-ETERN1«0CANNXI1.3) 

TR I  2,3 I-TNI2 ,3  I -ETERNl •CGANNXI2.3) 

TRl3,ll-TNI3,ll -ETERNl-OCANNXI  3,  l ) 
TRI3.2l-TNI3,2)-ETERNl»0CANNX(3,2) 

IF  l  COO  3 ICOTO  bC50 

TRI3, 31-Tni 3, 3|-ETcRNl*CCANNXI3,3|.ETfeftN2 
IFI.NOT.  C0C5 I CC TO  6G5C 
TR I  l ,  1 1 -TR I l . I  I -mnuP-0  T  Au3  3/El 
TRI2,2I-TRI2,2I*mnuP»OTAU33/EL 
TRI3, 31-TNI 3,31-CTAU 33/EL 
6050  CONTINUE 

TAUT-TRIl.ll-TRI2.2l-TRI3.3l 
YIELD  PM 

C2-TRI1, II* TR II, 11-2. -ITRI1.2I-TR12. ll-TRIl.3l-TRI3.il 

•  ♦TflI2,3|-TRI3,21l-TRl2.21-TRI2,2l-TRI3.31-TR|3,3l 

*  -I TAUT--2-2 .-S ICNSC 1/ 3. 

IF  IL  .CT.  LCICCTC  55 
IFIICOUNT  • E  0 •  l  I  CALL  AUXILI5I 

TEST  YIELO  CONDITION 

55  IFICT.LT. 0  .ANC.  U.EO. LCICCTC  220 
IF  IC2.LT. 0.1  CO  TO  121 
I  F I CC2 I  GOTO  60 

L-INTI YLCFAC-ISCRT I  I 1.5-C2-SICNSOI/SIGNSC1-I.OI l-l 

el-l 

QC2-.TRUE. 

IFIL  .CT.  L  I  CO  TC  50 
6C  CONTINUE 

PLASTIC  BEHAVIOR 


- 1  I F  C  CORRECTOR  3  OINENSIONAL  STRESSES  TCILA.lBI  as*  INITIAL 
OE  V  I  A  TOR  I C  STRESS  (FOR  THIS  INCREMENT  OR  SUB- I NCRE -Fn 7  I  oASEC 
TNILA.LBI 


STRE  S 160 
STRES161 
STRES162 
STRE  S163 
STRES16* 
S  T  RE  S l 65 
STRES166 
STRE  S 167 
STRESlbB 
STRE  S169 
STRE  SI  70 
ST9ES171 
STRES172 
SL8LASTRES173 
STRES175 
STRESL75 
STRES176 
STRES177 
5TRES175 
S  TRE  S  l  79 
STRES110 
STRE S 131 
STRE  S 182 
ST-ES133 
STRE  S  115 
STRES135 
SUESlob 
STRE  S l- 7 
ST»E  S  138 
STRES189 
STRES190 
STRE  S 19 1 
STRE  S 192 
STRES193 
ST»E  S 195 
STRES195 
S  T  o  E  S 1 96 
STRE  S 197 
STRE  S198 
STRE  S 199 
STRE  S200 
STRES201 
STRE  S202 
STRES203 
STRES205 
STRE  S205 
STRE  S206 
STRE  5207 
STRES208 
STRES209 
STRE  S2  10 
STRES211 
STRE  S212 
STRE  S2 13 
STRES215 
STRE  S215 
STRES216 
STRES217 
STRES218 
STRE  S219 
STRES220 
STRE  ,221 
STRES222 
STRE  S2  2  3 
STRE  S225 
STRES225 
STRE  S 226 
STRE5227 
STRFS228 
STRE  S229 
S 'RE  S230 
STRE  5  23  1 
STRES232 
STRE  S 2  3  3 
S  TRE  S2  3* 
STRE  S2  35 
S  TPF  S  2  36 
STRE  S2  37 
EC  ON  STRES238 
STRES239 


HTE«».HNUPP»TAUN 

CAUL  ERASE  (TC.9I 

TCII.I»«TN«1,11-hT£RN 

IE  (COO  A)  TCUiil'IClM  I-MNlP*TNI  3,  3) 

TC(l,2)«TN(l,2> 

TC(2.1l-TX(2.ll 

TCt2,2).TNi2,2>-NTERN 

IF  (OOO  A  I  TC(2,2l-Tu(2,2  l-AMP>I«(  3, 3  I 

IF  (C0C2I  CO  TO  7050 

TC(l»3)-TNIl,33 

TCI2.3I.TN12, 3) 

TCI3,ll-TNt  3,11 
TC(3,2).TNI3.2 3 
IF ( 000 3  .OR.  OGCAtcaro  7050 
TC(3,3 1-TNI3.3I-HTERN 
7050  CONTINUE 

TAOC»TC(l.n»TC!2.2l»TC(3,3» 

IF  (ISTRES  .NE.l)  CO  TO  70 
TANCENT  APPROACH 
TAUN3.TAUN/3. 

hlanoa.o. 

0CANAT-0CANNA/3.-ALPHAP0TEPP 
00  65  LA-L.3 
CO  65  L8-1.3 

65  hLAX0A.NLAX0a.|TN(la,L8I-TAUN3*DELTA(LA.L8>> 

•  *(0GANNK(U8  ,LA>-DGAXAT»CELTA(LB,LAI  ) 

HLAP0A-HLAP0A«EEP»l.5/(  S  ICPSOAEI.  t 

IF  (HLAPOAI  15C.150.L21 
70  CONTINUE 

COMPUTE  »NC  CHECK  VALUES  CF  A2.B2  AND  CISCRININANT 

A2*TC(l.l)*TCU.ll*2.*(TC(l,2)*TC(2.1t«TC(1.3l*TC!3.ll 

•  ♦TC(2,3l»TC(l,2l)»TC(2,2l*TC(2.8l»TCl3.1IPrC(3.3l-(TAUCA»2)/3. 
fl2-TC(l.ll«TR(l,l)*TC(1.2)*TR(2.1»»TC11.3»»TR(3.1) 

•  ♦TC(2.1»«TR(l,2)*TCI2.2l*TR(2.2>*TC(2.3)»TS(3.21 

•  ♦TC(3,l)»TR(l,3)»rC[3,2l»TR(2,3MTC13.3)»TR(3,3)-TAuC*TAur/3. 
0ISCR-82«B2-A2«C2 

IF (L  .GT.  L  1 1 1  2  }  LI  1 12  *L 
I F  (  ICOUNT  .EO.  UCALL  AUK  IL  (  6  1 

c 

C  TEST  a2 

C  IF  A /  IS  NEGATIVE  -  PRINT  ERROR  5ESSAGE 

C  IF  A2  IS  ZERO  -  SuE-INCR£N£NT 

C  IF  A 2  IS  POSITIVE  -  CONTINUE 

C 

IFIA2I  80. 150, ICO 
80  «R  I  TEINaRI  TE.,')C  I 

NO  F0RNAT11H  .ak.iamaa  NEGATIVE  AT  » 

GO  TO  180 
C 

C  test  DISCRIPINANT 

c 

c  IF  0 1  SCR  IS  NEGATIVE  -  SU8- INCREMENT 

C  OTHER. ISE  CONTINUE 

100  IFICISCR.lt. 0. I  CO  TO  150 

c 

C  TEST  92 

L 

C  IF  82  IS  NEGATIVE  OR  2E«0  -  SU8- INCREMENT 

C  OTHERWISE  CONTINUE 

C 

IF (82.LE.O. 1  GO  T  C  150 

C  CONFUTE  hlANOA  ANO  ELASTCPLASTIC  STRESSES 

HLAH02«  92*S0R  T (Cl  SCR  I 
HLAH0A-C2/MLAND2 
121  CONTINUE 

T"(l.l»-TR(l,l I -HU  A  NO A*  TC(  l  .  1  ) 

T«( 1 .2 l*TR( 1,2 1-HUANDAPTC ( l .21 
TN(2.ll«TR(2,l I -NLA  NO A*  TC ( 2,11 
TN(2.2»-TR(2.2I-HLAN0a«TC(2.21 
IF  ( COG  2  >  CO  TC  9050 
T»(l,3).TR(l,31-“LAN3A«TC(1.3) 

T*(2,3»*TR(2, 3  1-HL  A.ndavTC  (  2.11 
TN( 3,1 >-TR( 3,1 l-HLAHOA»TC(  3.1) 

TN(3.2I-TR(  3.2  )-HLA*<OA» TCI  3. 21 
If(0003IGOTO  8050 
TN1 1,31 -TR ( 3,3 1-mlanoa*  TC (3,31 
IFIQ00A1  TN( 3, 3 )-TR( 3, 3 > 

8050  CONTINUE 


STRES2-1 

STRES2*2 

STRES2A) 

STRE  S2aa 

SIRE  S2*5 

STRES2AO 

STRf  S2A  7 

STRES2A8 

S  TR  £  S  2  an 

STRE  S2*0 

S  T  R  E  S  2  5 1 

STRES252 

STRES253 

STRES25A 

STRES255 

STRES256 

STRES257 

STRE  S 2 5 8 

STRE  S259 

STRES260 

5TRES261 

STRES2b2 

STRE  S263 

STRES26A 

STRES265 

STRES266 

STRE  S267 

STRE  S25d 

STRES26R 

STRES270 

STRE  S271 

STRES272 

STRE  S273 

STRE  S 2 7 A 

STRES275 

STRE  S276 

STRES277 

STRES278 

STRES27<) 

STRES280 

STRE  S  28 1 

STRES2B2 

STRES253 

STR£  S28A 

STRES285 

STRES286 

STRES287 

STRE  S  2  9  5 

STRES299 

STRE  S260 

STRE  5231 

STRE  S2R2 

STRES263 

STRES29A 

STRES205 
STRE  S2R6 
STRES2R7 
STRES2R8 
S7RE520  7 
3  TRE  S300 
STRE  5301 
STRES302 
ST»E  S303 
STRE  S10A 
S  TRE  S  305 
STRES306 
S  TRE  S  30  7 
STRE  S  30 9 
STRE  S30N 
ST»E  S110 
ST»E  .311 
5  TRE  S  M2 
STOES’l  3 
STRE53IA 
STRE  S315 
STRE  S316 
STRE  S3  l  7 
STRE  5318 
STRE  S  3  19 
STRE  S320 


TAUH.THU,ll.TNI2,2l»THU.3l 
I F  <  I  COUNT  .EO.  1ICALL  tUXILWI 
CHECK  THE  SUB-INCRE“ENT  NUH8ER 

IF  (LC.EC.L IGCTO  210 
LC  *LC ♦ 1 
CO  TO  SO 

HAKE  SUB-INCREHENTS  5"*IL£R 


CHECK  H  AX  l HUH  NU  HBE  R  OE  ALLOWABLE  S UB- I NCRE HE N T S 


STRES321 
STRES322 
STRE  S  32  3 
STRES326 
STRE  S  3  25 
STRES326 
STRES327 
STRE S  32  8 
STRES329 
STRES330 
ST»ES33l 
STRES332 
STRES333 
STRES33A 
STRES335 
STRES336 
STRES337 
STRE  S338 
S  TR  E  S  3  3*7 


IFIL-IOOI  65, *.5.160  STRES336 

0  wRITE(NWRlTEtl7C)  STRES337 

3  EORHATUH  ,*.X,36HSTR£SS  CALCULATION  UNS  AT  I  SF  AC  TORT  AT  >  STRES338 

3  WRITECHWRITE.1901  I T  IHE  ■  1 1 .  12 . I L ATE R  .  I S e ,L ,LC  STRES339 

3  EORHATUH  .9X.9HTIME  S  TEE  .  1 1 . 5X  .  3H  1 1-,  1 2 ,5X,  3H  1 2»,  1 2 ,8  X  ,6HL  ATE  R- .  I  S  TRE  S  3*0 
♦2.9HSUBLAYER., 12,5X,2HL*.13.5X,3HLC..I31  STRE S3  6  1 

wRITE(NWRIT£,2C0l  A2 , 8 2 . C2 , C  I  SCR ,  STRES362 

•  ( <  THU  V,  JV  1  .  J  V- 1 ,  3  )  .  I  v*  1 .  3  I  ,  ( <  TR  i  I  V.  JV  1  ,  STRES363 

♦  JVl,3»Uy*l,3l,((TC(  IV,  JV),JV.L,3ltI<-l.31  STRESS*.*, 

3  F0RHATU0X,"A2  -".E 15 .6 . 10X,"B2  •  •*.  e  1 5.6. 10X  ,  "C2  ■  "  ,  1 15 . 6S  TRE  S  3  <.5 

•.10x,"0ISCR  -".E15.6  STRES366 

*/10X,"HIX£C  INITIAL  STRESSES". 16X,“-  in  •" . 3E 15.6 75hx , 3E 15 .6/  5TRES3‘7 
•56X,3£15.6/10X,"HIXED  TRIAL  ELASTIC  ST  RE SSE S “ .  10X ,  TR  STRES368 

•3E15.6/56X,  3E15 .6/ 56X,  3E  15.E/10X  .“HI  XEC  CORRECTOR  3-C  I  HENS  I  CNAl  STST»ES3*.9 


•RESSES*  TC  ■“.3E15.6/56X.3E15.6/56X.3E15.61  STRE  S  350 

CALL  0IACNCI5I  STRES351 

STRES352 

HAVE  REACHEO  PLASTIC  SOLUTION  ST»ES353 

STRES35*. 

3  CONTINUE  STRES355 

IPLASTUl,I2|.tPLASTUl,l2l»l  STRES356 

CO  TO  222  STRES357 

3  CONTINUE  S  TRE  S  358 

STRE  S  359 

ELASTIC  TH  equals  trial  TR  PER  SUBLAYER  S  TRE  S  360 

STRES361 

THU,ll«TR( 1,1)  STRES362 

THU.21-TRT1.2)  S  TRf  S  36  3 

THT2, ll.TRT2.il  STRES36A 

THI2.2I.TR12.2I  S  TRE  S  365 

IF  (Q0C2)  GO  TC  9080  STRES36N 

THU,3t.TRU,3|  STRES367 

THI2,3|.TRI2,3I  S  TRE  S  368 

THU.ll.TRU.il  S  TRE  S  369 

THU.2l.TRU.2t  STRES370 

THT3.3t.TRU.31  STRES371 

3  CONTINUE  STRES372 

STRE  S373 

7  CONTINUE  ST»E83?*. 

CGAHHXT  3,3  l«  ( (  rH(3,3|-TM3.TI  1-hnu*  IT  TH(  l .  1  l-TN  U  ,  l  I  I  ♦(  TH  (  2 , 2 1- TNT  STRE  S  3  75 
•  2 .2  t  I  1  !/ES»!l.  .HNut*(  2.  MM  3. 3  l-TNU,  l  I-TM2.2  1  l*HLA"CA/  t  3.  «ES  I  ST JE  S  3  76 
IF  !  I  COUNT  .EC.  1  t  CALL  ALXIL  U1  S  T  =  E  S  3  7  7 

5TRES378 

STORE  HIXEO  TENSOR  SUBLAYER  STRESSES  STRES379 

I F  T I  COUNT  .EO.  0  t  CO  TO  2500  S  TR  E  S  3  80 

STRES381 

TAUllTIl.I2.I21.THTl.il  STRE  S  382 

TA012T  11.12.  I2I-THU,  21  STRE  S  383 

TAU21lll.l2.l2fTHT2.il  STRES38A 

TAU22TIl.I2.I7l.THT2.21  STRES3B5 

IF (QQQ2 1C0TC  2500  S  TRE  S  386 

TAU13T I l .12, 12 1*TH( 1,31  STRES387 

TAU23Ul.I2.I7l.THI2.3l  STRES388 

TAU31II1.I2.12I.THT3.11  STRES389 

TAU32TI1.I2.I2I.THI3.21  STRES390 

IFT0CQ3IG0T0  25CO  STRES391 

TAU33T 1 1.12, 12 I'THll.J I  ST»E  S3<»2 

3  CONTINUE  *T»t  S,',J 


TAUF(2.l)*TAUE  (2.1)*C0EEE(  tSBMTttl2.il 
TAUF(2.2)*TAUEI2,21*C0£FFUS8I*T"(2.2) 

IF  < CQQ2 ) GOTO  2501 

TAUFIl,3l*TAUF(l.3)*CaEEEUSBl*Ttt(l,3) 

T  AuF ( 3, 1 l*T  AuF (3,i)*CQEFFU5B)*TP(3,ll 
TAUE(2.3)*TAUE(2,3)»C0EEF(  ISB)*TP<2.  31 
TAUF  1 3,2  * «T AUF  13.21  ♦COEFF  use  3,21 

IE(0OO3)G0T0  2501 

TAUF (3,3 I* TAUF ( 3,3) *C0EFF( ISS)«Ttt«3.3l 
2501  CONTINUE 

CGAP33*CCAP33*C0EFF(ISBI*0CAPP*(3.3) 

ISUBL* (SB 

IF  (  l  COUNT  ,EQ,  l 1 C  ALL  AUAILI31 
C 

250  CONTINUE 

c  CONTRAVAR I  ANT  TOTAL  STRESSES 

TAUI 1, 1 )*TAUF  (  1  , 11*00(1,  II • T  AUF M .2  I *001 2 • 1  I  *T AUF I l  .  3  I *CC (  3 .  1  1 
TAUU.2  1-TAUFl  1 ,  1 )*GG(  1 .2  »*  T  AUF  I t .2 1*0012  ,2  )  ♦  T  AuF  U  .  3  )  •  GO  I  3.2) 
TAUI  2, 1  »*TAUF  I  2  •  II  *GG1  l,l)*TAUF(2,2)*GGI2,l)*TA0FI2«  3 1  *CC (  3.1) 
TAUI2.2 1-TAUF 12,1 »*CCI 1.21 ♦ TAUF I  2 ,2 1 «CC (2 .2 1 ♦ T AuF ( 2 , 3 1 *00 C  3 ,2  I 
TAU(X,21*.5»ITAUIl.2»*TAUI2.1ll 
TAUI2.11*TAU(1,2J 
I F ( QCQ2 I  GOTO  2600 

TAUI1.3MTAUFCI.  11*0011, 3l*TAuF(l.2)«(.Gt2,3l*TAUFtl.3)*CG(3.3) 
TAU I  3 , 1  I »TAUF ( 3 , 1 1  ACC  I  1 , l » ♦ TAUF ( 3 ,2 )*GC (2 . * ) *T  AUF! 3, 3 1 *00 l 3 , l I 

TAU ( 2 , 3  > *TAUF ( 2 , 1 1 *00(1 «  3 1 *T AuF ( 2 ,2  I *00 ( 2 , J l ♦ T  AuF ( 2 , 3  I *GG l 3 , 3  t 
TAU(3.2l*TAUF(3,l)*GG(1.2l*TAUFI3.2)»GGI2.2)»IAUFI3.3)*CG(J,2l 

TAU(3.3I*TAUF (3, 11*00(1.1) ♦ TAUF 13.2 1*00(2, 3)*TAoF (3. 31*0013, 31 
TAUtl,3»*.5*(TAU(l,3!*TAU13.l)» 

TAu ( 3. 1  l-T AU( 1 . 3 1 
TAUI2, 3I*.5*(TAU(2, 31* TAUI 3.2) I 
TAU  13.2 I-TAUI2, 31 
2600  CONTINUE 

IF(  ICOUNT  .EO.  C100TC  2)0 
ANUP-1HF 

IF (  IPLASTC I  1  .  I  2  I  .OT.  0)  ANUP-1HR 
T AuS  PH* ( TAUF ( 1 ,1 1 »TAuF (2 ,2 ) *TAUF ( 3 . 3  I )/ 3,0 
IF(»AUXIL  .EO.  1)  CALL  *AXAIN  I  TAUI-  ,  1) 

CALL  AUKIL  (ID 
00  260  LL*l.NUr 


STRE  SAOl 
ST»E  SA02 
STRE$*03 
STRE  SAG* 
STRESA05 
STRESA06 
STRESaO? 
S  TRE  S 5 0 8 
STRE  S*OS 
STRES*10 
STRE  S  * 1 1 
STRE  S*12 
ST*E  3*1  3 
STRE  5*1* 
STRE  S  *  1 5 
S  T  »  t  S  *  1  6 
STRE  S*l  7 
STRE S  *  1 3 
STRE  5  *  Is 
STRE  S*20 
STRE  S*21 
STRE S  *2  2 
STRE  S*2  3 
ST»tS*2* 
STRE  S  *  25 
STRES*2b 
STRE  S  A 2  7 
ST»E5*28 
STRESA2S 
STRE  S*30 
STRE  S*3  l 
STRESA32 
STRESA33 
STRE S  »  3  A 
STRE  5*  35 
STRE  S  A  36 
ST  =  ES*3? 
STRE  S  A  3  3 
STRE  S*  30 
STRE  S**0 


IF(  11  .NE.  IPSKLLI  ICOTC  260 
IF(I2  .NE.  IRS2(LL))GUTC  260 
TAuSPL(LL.ICauSS)*TauS?f 

L2-ICAUSS*NCAUSS(  (layer  I  *(LL-l>*I(-A»6R 
00  255  1*1,3 
00  255  J» l , 3 
TAUPIL2.I , J)*TAUF( I, J) 

255  CONTINUE 
260  CONTINUE 
C 

FPATIICAUS5,ll,l2)*ANuP 
270  CONTINUE 

C  CALCULATE  CHANCE  IN  TP  PER  GAUSS  STATION  AS  TP  -  TN 

C 

CTP(I,I)*GT*(l,l).TAU(l,l)*C(i,l)*TAU(l. 21*0(2, 1)*TAU(1, 31*0(3,1) 
0TPI1, 2) *0TP(1, 2)*TAU(l.l)*C(l.2l*EAu(l. 21*0(2, 2 l*TAUt l , 31*0(1, 21 
DTPI2,II*0TP(2,1I*TAUI2,II*C(1,1)*Eau(2.2)*0(2«1)*TAUI2,3)*0(3.1) 
0TP(2,2)*0T*(2,2)*TAU(2.1I*C(1,2)*Tau(2,2)*C(2,2)*TAu(2,3)*C(3.2) 

IE  (00C2)  GO  TO  SO  70 

OT"(I,3)*OTP(l,3l*TAU(l,  1I«CI1,  1 1  ♦  T  AU(  l  ,  21*0(2,  3  I *T  AUI 1  ,  31*0(3.3) 
0TP(2,3)*0Tp|2,3I*TAU(2,1)*CI1.3)*Tau(2,2)*CI2.3)*TAU(2,3)*C(3,3) 
OTN(3,l)*OTP(3,l)*TAUl3.l)*Cll,l)*TAU(3«2)*C(2.l)*TAU(3,3)*C(3.1) 
0TP(3,2)*0T“(3,2J*TAU(3,l)*C(i,2)*TAU(3,2i*CI2.2l*TAU(3,3)*C(3.2) 
OTP(3,3l*OTPt3,3l*TAU<3,l)*C(l.3>*TAut3.2)*GI2.)l»rAu(3,3)»GIJ.l) 
SO 70  CONTINUE 
RETURN 
END 


STRESaaI 
ST»E  S A * 2 
STRESaa) 
S  TRE  S  AAA 
STRE  S**5 
STRE  S * *6 
STRESAA7 
STRE  Saab 
STRESaas 
STRES*50 
STRE  S  *5  l 
STRE  5*52 
STRE  SA53 
STRE  5*5* 
STRE  S*5S 
STRE  S *56 
STRE  S  *  5 7 
STRE  S*5 j 
S  T  o  E S*5S 
ST5ESA60 
STRE  S*6  1 
STRE  SA62 
STRES*63 
STRESA6A 
S  T»E  5*65 
STRE  S  *66 
ST«E  S*6  7 


78 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  JET*  JET* 

JET* 

EVALUATE  STRESSES  THROUGH  THE  THICKNESS  ANO  EV»LU»TE  JET* 

contravariant  copponents  of  the  relative  hones!  resultant  tenscr  jet* 

*NO  OF  THE  RELATIVE  STRESS  RESULTANT  TENSOR  JETa 

JE  T* 

IPPLICIT  LOGICAL (Q)  JET* 

JET* 

COPPON/  ALLENE  /TOTAL. TOTKIN,TOT£LA,TQTPLA.TOT*EX,  TOTTEp,  INERGT  JETa 

•  .  TOTVIS.TOTEl.TOTEJ  .ECHECK  /eta 

•  ,0TP|3,3I .SQRG.SCRA.SIGHSQ.AJ.BJ.CJ  JET* 

JETa 

COPPON  /APTRAC/  KSPOTS.KTt * > ,XDIRI *, 3) .PPTYPE I  *  I ,PT TYPE  I  *  I , XSPCT ,  JETa 

•  KSIOE.  STIFMA,  31,  TORSAIAl  /eta 

COPPON  /APTRAL/  TA1 I  3  ,801 , TA Jt 3.80 1 ,M«l 1(90  I  .MP22 ISO )  JET* 

level  2,TR1,TA2,hp11,«P22  /ETa 

JETa 

COUPON  /CARTE/  YTEST.YnER , YS AvE  /ETa 

COPPON  /CARTEL/  y(3. 20. 20),  CI3. 20, 203.Y2I  3.20.201  ,021  3 .20,20  JETa 

LEVEL  2.Y.C.Y2.C2  /ETa 

/  ETa 

COPPON  /CTIPE/  *U*< 20) . T I»£ ,OEL TAT , TIPEF , I T IPE , I T I PEF . I Au« < 20  I .  JET* 

•  IQUT(2C),0PRINT(20)  JETa 

COPPON  /CTIPEL/  IPLAST«20,20),P(20.20I.PPlI20.201  JETa 

LEVEL  2,  IPLAST.P.PPL  JETA 

JETa 

COPPON/CT HKN/  2AI20  .20) .281 20,20) ,D2A1 120.20) .02*2120 .20) .  JETa 

•  0JB1I2C.20 ) ,0202 (20,20  I ,88AR  ISO )  JETa 

LEVEL  2, 2A*2B«C2A1,Q2A2, 0281,0282, 3BAR  JETa 

jni 

COPPON  / INDEX/  NREAO.NvR I T£ ,NPUNCH,NAESH1 ,NPESp2,NI,N2,N/,n1P,n2p,2£Ta 

•  Nl PP.N2PP.il, 12, IJ.ILJ£Ra.I2JER0.IRYL,IRY2.iRY3.I<YA,ISTRl,IST«2./ETi 

•  ISTR3*ISTRa,1CL.1C2,IC1,IC2,IPI,1P2*1SI,IS2.R1,R2.<3,Ra,kRln,  JETa 

•  k2STQP,KYTEST,!CIR,I1TEST,I2TEST,KIN1TL  JETa 

JETa 

COPPON  /N£  v 1/  CCISP.POISP.I5Ti.IST2. 1ST). 1ST*. CTRaC.«T6aC,  JETa 

•  T  AUPAX,TAUPIN,GAP«AX,GAPPIN.ELP»X,ELPIN,  UPAXS,  I2“AXS  ,  I  GAP  AS.  JE  TA 

•  I L APRS, I  IP  INS. I2PINS. IGAPIS.lLA'IS. I  IP  A  XX, IJPAxx ,1 JPAXX. IIPINX.  JET* 

•  I2PINX,I2P!NX,IIPPA*,I2PP*X,IJPPAX.IIP»IN,I2P»IN,IJ*-|N,UI(2C).  JET* 

«  OTRACT.CTRACP.1PHVI20I.PPHYS  JET* 

JETa 

COPPON  /OPTION/  P  Aux  It. ,  P  I NC  EC.PINVEL.PLOAO.PPATPR.PSPLCA,  JETa 

•  PSPTEP.PTEPPE.PTHIKL.PIPPUL.ISTRES.INCRPL.ISTREJ  JETa 

JETa 

COPPON  /PHYSCN/  EE ,HNU. ALPHA, CONST, EXPCN, FACTOR. RATE ,«HC.  JETa 

«  MLAPOA.COEFFI ‘ ) ,s ICPAI 5  I , !P« 3.3) , TCI  3, 3) .OELTai 3 ,3)  JETa 

JETA 

COPPON  /PUSH/  FORCE  5 ( 3  ).  VELCCI 3 ) .RAT  10, RAT  10" .Cx l ,0X2 . TE *P ,CTE 'P,  JETa 

•  FSPACE ,TS?*CE .F INCNO.FSTOP.TSrCP.THCUEF  JETA 

COPPCN  /  PuShl/  SCRAT|20«20),SORAJ(20.20).FPASIU20.20),  JETa 

•  FPA522(20,201,FPAS23<2C,2C),F*A533(2Q,20>  JETa 

LEVEL  2  .SCRAT.SCR AJ  ,F«AS  l  l.F*AS22 .FPAS2  3  ,F*AS3 3  JETA 

JETA 

COPPON  / CLCGIC/  CAUXI2C  I  ,OJETA,OSTRE S.CPLAST . JSENSl ,0E Cu IL.  JET* 

•  QO  IAGN.OINCEO.C  INVEL.CLQAC  .CPA  TP  R  ,  C  TH  I  «  L  ,  U  TE  P  PE  .CSPTEP.CAUXU,  JETA 

•  QAUXU.CSPLOA.CIPPUL.CSfaRP.OPESO.OIRCH.CShEAR  JETa 

JETA 

CO“PON  /SuRNOP/  SNPRI3I  JETa 

COPPON  /SURNOL/  SN  <  3,20,201  JET* 

LEVEL  2 , S N  JETA 

/ETA 

COP*ON  / S 2 /  STRESS  I  3 )  JETa 

COPPON  / S2L /  STRESL (  3.20 ,20  ) ,STRE SGI  2 .20,20 ) , S t RESP l J ,20.20  )  JETa 

LEVEL  2, STRESL. S TRE SC , 5 TRE 5 P  JETA 

JETa 

COPPON  /TENCOP/  YY I  3,2 ) , YYYI 3,2.2 ) , Al 3 , 3 ) ,81 3. 3  I  ,AA(3 . 3  )  ,881 3, 3  )  ,  JETa 

•  8PI3, 3), 0*13,31,0813, 3), GI3, 3). GGI3.il, ON  131, COI 3 .2  I .CCC I i ,2 ,2  I ,  JETa 

•  OGAPI  3, 3  I  .OCAPPX  I  3  ,  3  I  ,HM  3 , 3)  ,HC(  J  )  ,  TAUI  3  ,  3  )  ,  TALS8L  13  ,  3  I  JET* 

•  ,0021  3,21.00021  3.2.2  1  ,YY2l3.21.YYY2li,2,2).YY'jn,2).0BP(3.3l  JETA 

JE  TA 

COPPCN  / THP  NS /  2 ,22 ,2CENTR, THICKN.ABSC l S (6.6) . »E  I  CMT (6.6  I  ,  JETA 

•  NGAUSSI *) , IGAUSS, PLAYER , IL AYER .NSUBL t *1 . I SUBl  JETa 

JETA 

COPPON  JTNCOPP/  HP  II J,  20 ,20  I ,HPJ I J ,20,20 ) .CAPGl I  2 .20  .JO  I ,  JET* 

•  CAPQ2(2,20,20),CAPQi(2.20,20l,CAP2Cl(2.20«20>.  JETa 

•  CAP202I2,20,2C) , LAP  20  3  I  2 , 20 . 20  I  JETa 

LEVEL  2,HPl,HP2»CAPQi,CAPC2,CAPC3,CAP2Cl«CAP2u2.CAP203  JETA 

JETA 
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3 
6 

5 

6 
7 

e 

9 

10 
1 1 
12 
13 
16 
1* 
16 
l  7 
18 

19 

20 
21 
22 
23 
2* 

25 

26 

27 

28 

29 

30 

31 

32 

33 

36 
35 
3b 

37 

38 

39 
60 
61 
62 
63 

6  R* 

65 
6b 
6  7 
68 
69 

50 

51 

52 

53 
5«» 
55 
5h 

57 

58 

59 

60 
61 
62 
63 

66 

65 

66 

6  7 
68 

69 

70 

71 

72 

73 

76 
75 

7  6 

77 
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COHTGN  /VARTHK/  VTMlM  6*20*20  I .CENTRVt* ,2G, 201  26T&  ?9 

LEVEL  2*  VT>f  v  7ETA  80 
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i  ^  v*v»v* 


WAll  [  *  W  mi1  --1 1-1  >. '  •.'■  *  U  v.k  j.it.m  g 


Bia 


vfts 


:jS 

& 


►•l 


COPPON  /VI S 1/  IVIS.AV15I3I ,CVIS( 31 ,OSTRll,OSTR l2.QSTR13,OSTRZl. 

♦  DSTR22,OSTR23,OSTR3l.D5fR32,OSTR33 

COMMON  ZVISILZ  VSTRIL(3,20,20I.VSTR1213,20.20I.VSTR13«3.20.2CI, 

♦  VSTR22(3,20.2C1,VSTR23(3,20.20> 

LEVEL  2»VSTRll*VSTRl2»VSTR13,vSTR22,VSTR23 

COPPON  /STRINT/  EPSLII20.20),EPSL2C20.20),GRppal(20,20I. 

♦  EPSU1(20.20)«EPSU2(20.20I,GAPP*U(20,20> 

LEVEL  2.EPSL1.EP5L2  .GAPPAL .£ PSU l ,£ P Su2 .GAPPAu 

COUPON  /PA  II/  NCONT ,NR I TE .STRAIN, PR ITE .NCEL? .ET A01.ETAC2 .NS I«N 

♦  ,  P ACTOP  ,PQ , FAC  TON. NQ , f K I N 
COPPCN  /CNE/  OCAP33 
COPPON  /THREE/  A  CAP  3  3 
COPPON  /TEPP/  P GAP 3  3 

COPPON  /OELC/  (COUNT 
COPPON  /CEL/  0EL8ARI20.2C) 

LEVEL  Z.CELBAR 
COPPON  /GG/SURFCC 

COPPON  /CPRlNTv  CTYPE .C8ASE < 3. 3  I .CUBASE ( 3 . 3 > 

COPPON  /ERA C/  TAUF l  3.31  .TAUSPM.NUP. IPSll 10J . IPS2I 10  I 
DIPENSICN  C2(3I,CSI3,2I.Cx(3,3>,hp(2,2> 

PAR  >0. 

FOR  ORTHOGONAL  ccorcinates.  to  »ECuC£  error 

IF ( K 1 .NE . 7 . ANC.K2 .he .? )  GO  TO  10 
OA  1 1,2 )-0. 

0A(2,1)-C. 

OBI  1,2  1-0. 

0B(2,l)-0. 

10  CONTINUE 

CALL  EfiASE(H-,3l 
STRESLIl. 11,121-0. 

STRESLI2. 11,121-0. 

STRESLI 3,(1.121-0. 

CAPC1I 1,11,(21-0. 

CAP01 12 , 1 1 , I  2 ) -C. 

C  A  PC2  t 1,(1,121-0. 

C  APC2  <  2 , 11,12  1*0. 

C  A  P  C  3 ( 1,11,(21-0. 

CAP03»  2 , I  1, (2) -0. 

CAP2GK1, II, 121-0. 

CAP2Qlt2,H,(2)-0. 

CAP2Q21 1,(1,121-0. 

CAP2Q2(2,Il.I2)-0. 

CAP 203  I 1, ( 1, 12 1-0. 

CAP2Q312, 1 1 , 12  I -0. 

5TRESPI1.U.I2I-C, 

STRESP(2,Il.I2)-0. 

3TRES0Cl.il, 121-0. 

STR  ES0C2 , 1 1 . (2 1-0. 

C2CU-C. 

02(21-0. 

02(31-0. 

22-0. 

08(3,31-0. 

IPLAST<  11,(21-0 
- ( C-SURF  ACE  OFFSET 

FOR  KIRCHHQFF  GR  SHEAR  SHELL  CG**SIGER  CEL-0  SINCE,  F»OP 
The  POINT  2-0  IS  The  half  THICKNESS  OF  THE  Shell  anO  IS 
IN  The  OENSIIY 


.CR.GSHEARI  GO  to  2205 


•1,3 


VARIATION 
CEL-O. C 
IF  (OIRCH 
THIK2-C. 

00  220*  ILATER 
CALL  THIKLAC1I 
220*.  THIK2-THIK2»TMIC«N 
ZCEN-Th|kZ*.8 
IL  AVER-  1 
CALL  THlKLAC  11 
22CEN-THICKN-2B ( 11,121 
0EL-2CEN-/2CEN 
2203  CONTINUE 

CS<l,ll-eP(l,ll*VY(l,ll.8P<2.H-vr<l,2) 

CS(l.21-BP(l,2l-YY(l,ll.0P(2.2l-»Y(l,2l 

CS(2,ll-RP(l,ll»YYI2,ll,BPI2.ll»YY(2,2l 

CS(2,21-BP(1,2I*yy(2,11,8P(2.21»Yyi2,2I 

CS(3,ll-BP(l,ll»YY(3.ll-RP(2,ll»YY13,2l 


2  E  T  A 
/ETA 
2  E  T  A 
2E  TA 
2ETA 
2ETA 
2  ETA 
2ETA 
/  ETA 
/E  TA 
2  ETA 
2E  TA 
/  E  T  A 
2E  TA 
2E  TA 
ZE  TA 
ZETA 
ZETA 
ZETA 
ZETA 
ZE  TA 
ZETA 
ZETA 
Z  E  Ta 
ZETA 
ZETA 
ZETA 
ZE  TA 
ZETA 
ZETA 
ZETA 
ZE  TA 
ZETA 
ZE  TA 

ZETA 

ZETA 

ZETA 

ZETA 

ZETA 

ZETA 

Z  eta 
ZE  Ta 

ZETA 

;e  ta 
/‘ta 

ZETA 
ZETA 
ZE  TA 
ZETA 
ZETA 
ZETA 
ZETA 
ZETA 
ZETA 
ZETA 
ZETA 
ZETA 
ZETA 
ZETA 
ZETA 
ZETA 
ZETA 

TmIkla,  ZETA 
NOT  -tIGHEZETA 
ZETA 

ZE  ta 
ZE  Ta 
ZE  TA 

ZE  ta 

ZETA 

ZE  Ta 
ZE  Ta 
ZE  Ta 
ZETA 
ZETA 
ZE  TA 
ZETA 
ZETA 
ZE  TA 
ZETA 
ZETA 
ZETA 


•fl 

92 

91 

83 

93 
96 

92 
38 
96 
60 
61 
62 
63 
63 
63 
6b 
67 

6e 

66 

100 

101 

102 

103 

103 

105 
10b 
107 

109 

106 

110 
111 
1  1  2 
113 
l  1  3 

115 
lib 

117 

118 

116 
120 
121 
12  2 

12  3 
123 

125 

126 

127 

128 
126 

130 

131 
1  32 
133 
1  33 
1  35 

136 

137 
118 
136 
130 

13  1 

132 

133 

1  33 
136 
I  -6 
13  7 
1  38 
136 
160 
16  1 
t‘2 
1  6  3 

163 
136 
1  66 

167 

168 
136 
160 
161 
162 


Lr.V». 
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CS(3,2l-BN(l,2)-YY(3,il.BB(2.2)-YY(3,Z) 

ZETA 

163 

call  erase  ihn.ri 

ZE  TA 

16* 

Kltllh  U.  12  1-0.0 

ZETA 

165 

HB1(2,  l  l ,  1 2  )  -0 . 0 

zfta 

166 

MB2ll.il, 121-0.0 

ZE  TA 

167 

HB2I2. 11.121-0.0 

Zf  Ta 

168 

12-0 

ZETA 

169 

IF1IC0LNT  .£0.  1)  CAUL  £ N£ 9 G Y ( 2 , l . 0 1 

ZETA 

170 

I  F  I  I  NOR  PL .E  C.2 )  0041*43*0. 

ZETA 

171 

IF( IN0RBL.EC.2)  THICK/* 0. 

ZETA 

I  72 

DGB33-C.O 

ZETA 

173 

SUBG-0.0 

ZETA 

1  7* 

TMIC-0.0 

ZETA 

176 

SGAB33-0.0 

ZETA 

1  76 

CALCULATIONS  FCR  ALL  LAYERS 

ZETA 

1  77 

00  169  ILAYER- 1 .NLAYES 

ZETA 

1  7  8 

ISTRE/- ISTRES 

ZETA 

179 

IF  (01RCH)  CC  TO  63 

ZETA 

100 

IF  (NLAYER.EO.il  00  TO  si 

ZETA 

181 

IFIILAYER-21  60.61,62 

ZETA 

18  2 

60 

ZZ- ZB  1 11.12  1 

ZETA 

103 

OZ  1  l  l-CZRl ( 1 1.  12  1 

ZETA 

10* 

OZ  12 ) -C  Z82 I  11.12) 

ZETA 

185 

OZ ( 31-0. 

ZE  TA 

1  06 

00  TO  63 

ZE  Ta 

1  0  7 

61 

COST 1N0E 

ZETA 

188 

1FI  .NOT  .QIRCH.ANC. .NOT .OSMEAR.ANO. I STRE  S.  £  0 . 2  )  ISTREZ-3 

ZETA 

189 

OZ (11-0. 

ZETA 

190 

0ZI2I-0. 

ZETA 

191 

02(31*1. 

ZETA 

192 

CO  TO  63 

ZETA 

193 

62 

ZZ-ZAII1.I2) 

ZETA 

19* 

OZ  I  ll-CZAK  1 1,  12  1 

ZETA 

195 

02(21-02 A2I 11.12) 

ZETA 

196 

02(31-0. 

ZETA 

197 

63 

CONTINUE 

ZETA 

198 

CALL  THIKLA(l) 

ZETA 

199 

IF ( IN0RBL.E0.2I  TM1CKZ-TMICKZ*THICKN 

ZETA 

200 

TMIC-ThlC-THICKN 

ZETA 

201 

NCAUSL-NCAUSS(  ILAYER  ) 

ZETA 

202 

Surface  TRaxon 

ZETA 

203 

CALL  SFORCE  (CS.TMICKN.yy.YY2, 11, 12, SuRFOO) 

ZETA 

20* 

CALCULATIONS  FCR  EACH  CAUSS  POINT 

ZETA 

205 

CO  169  IOAUSS-1.NOAUSL 

ZETA 

206 

Z  AS  BEASuREO  FRCB  BIO-SURFACE 

ZETA 

207 

Z-.S-TMICKN-ABSCISI I GAU  5  5 , NO AUSL ) -ZCENTR 

ZETA 

208 

z  as  beasureo  frob  reference  surface 

ZETA 

209 

Z-Z-OEL 

ZETA 

210 

IF  ( ILAYER. £0.2. CR. ( .NOT. CIRCH.ANO.NLAYER. EO.l 1 1  22 - Z 

ZETA 

211 

IF  (CSHEARI  ZZ - Z 

ZETA 

212 

IF  (BTEBPE  .NE.  01  CALL  TEMPER 

ZETA 

213 

IF  (BTEBPE  .NE.  0  .OR.  ICAUSS  .EQ.  1)  CALL  baTPRO 

ZETA 

21* 

if(btebpe.ne.o  .or.  icauss.ec.u  call  enercyis.i.oi 

ZETA 

215 

CBTN-0.0 

ZETA 

216 

CALCULATIONS  for  all  components 

ZETA 

217 

00  9  J-1,3 

ZETA 

218 

CBASE( 3 , JI-SN( J, I l, 121 

ZETA 

219 

IFI.NOT.CSMEAR.ANU.INORBL.EC.Z)  C8aSE(3.JI-CBASE(J.JI-Y2(j,I1, 

>12)  ZETA 

220 

IF  (CIRCH1  CO  TC  9891 

ZETA 

221 

IF (  .NOT .CSHE AR . ANO. INORBL .EC.Z 1  CO  TO  9891 

ZETA 

222 

GBaSE(3.JI-G8ASE(3.JI  ♦CZ(31»Y2(J,Il,t2) 

ZETA 

223 

9891 

CONTINUE 

ZETA 

22* 

CBTN-CeTN-CBASE ( 3. J 1 *SN«  J. 1 1 .12 ) 

ZETA 

225 

CBASf(l.J)-YY(J,l)-Z*CS(J,l) 

ZETA 

226 

CBASE(2.JI-YY(J.21-Z»CS(J,2) 

ZETA 

227 

IFICIRCMI  CO  TO  9 

ZETA 

228 

GBASE(l.JI-G8ASeii,J>»CZ!ll-Y2(.i,ll,l2l»22-YY2<J,ll 

L  ETA 

229 

CBASE( 2 . J>-GBAS£ (2, JI-CZI 2>-Y2( J. U . 12  >  »22-YY2(  J.2  1 

ZE  TA 

2  30 

9 

CONTINUE 

Z  E  TA 

231 

CALL  ERASE ( CG8ASE ,9 1 

ZE  ta 

232 

OCBASE ( 3, 1 1 -ONI  1 ) 

/  ETA 

233 

CGBASE (3.21  -  ON ( 2 1 

ZETA 

23* 

OCBASE ( 3 . 31 -0N( 3  1 

ZETA 

235 

IFICIRCMI  CO  TC  806 

ZETA 

236 

OCBASE (3. 1I-0CBASE<3, 11-02(3) *0211. II. 121 

ZETA 

237 

CCaASE(3,2l-0C8ASE(),2).CZ(3l*02(2.Il.I2l 

ZE  TA 

238 

0C8ASE(3,3l-OCaASEI3.3).0Z(3l«O2(),ll,I2) 

ZETA 

2  39 

886 

CONTINUE 

ZETA 

2*0 

81 


oo  aas  la»i,2 
oo  aas  k«i.3 
00  887  LS*l,2 

887  DC8A5E<LA,K).0CBAS£tLA.K|,<C£LTA(LS.lA)-2»8P(LS.LAI)*0DIK.LS>- 

•  2A08*<L5.LAI«irYlK ,LS l-OOIK ,LS  » t 
If  tOIfcCH)  CO  TO  Sad 

DGBASE(LA,K|,0C8ASE(LA,K|.22*QC2I«.LA).C21LAI*02(K,[1.!2) 

888  CONTINLE 

G<l.l>-C8AS£ll,ll»GaASE<l.ll,CeAS£<t,2)*GBASE(I.2l 

•  »C8AS£ (I . 3 )»CBASE ( 1  .  3  I 

ct i »7)*cbase  a .1 i«g8ase i; ,i >  *ceASE 1 1  .£ i AcaAse i a  .7 > 

•  •GBASE(1.3I»GBASE12.31 

G<l.3)-C8ASE(l,ll«CaASE<3.1l  *G8ASE  <1.2I»GbaSE11,2I 

•  *G8ASE  <1.31 »CBASF  <3.31 
C<2.1)>G< 1.2) 

Gt2.2)-CBASEI2.1)*GaASE<2.1)»GaASE<2.2)»CaASE<2.2) 

•  aCBASE <2 . 3 ) »G8 ASE < 2. 3  I 

G<2.3I>G8ASE<2i 1 1 *G8 AS E <  3 . 1 1 ‘G8 1SE  < 2 . 2  » «G8 AS E < 3.2  1 

•  ‘GBASE <2,3 l»GBAS£< 3.31 
G<3.l»*G<1.3) 

G<3.2)-C<2.3) 

G<3.3l-C8ASE<3.ll«C8AS£(3.1> ‘G8ASE  <3.2)a&oaSE(3.2> 

•  ♦G8AS£<3,3I*G8ASE< 3,31 
CALL  ERASE ( OGAP  ,9  I 

I  J-3 

IF(OIRCH)  [j-2 
00  '102  I  *  1 «  I  J 
00  102  J-l.lJ 
00  102  x«1.3 

102  OGA(*<I,JI-CGAH<|,J1,.5«<GBASEI1.K>»OG8ASE<J.K|,G8ASE<J.|‘I*OGBASE< 

•  I .K I-CC8ASE < I ,K l«OGBASE < J.K  I  I 
G-CTYPE 

GTYPt-G<l,ll»lG(2.2>»C<3.3)-CI2,3l»*2>,Cll,2)»lC(I.3>»G(2,31- 

•  G<l,2)»C(3,3IJ,G(l,3l«<G(l.2)*Gi2.3)-C(1.3)»G(2.2>) 
SORG-SCRTIGTYPE ) 

GTYPE- l./GTYPE 

GG<l,ll»<GI2,2>«G<3.3)-G(2.3»**2)*CTYP£ 
GGIl,2)-<G<l,3»«G<2.31-C<l,2)»G<3.3l)*GTYPE 
GG I  2 , 1 1 "GG (  1,21 

GG<l,l)-<G<l,21*C<2.3l-C<l,3)»G<2.2n*GTYPfc 
GC<  3.1t*GG< 1,31 

GG<2,2>*<C<l.U»G<3,31-C<l,3»»«21»CTYPE 

CG<2,3)*<G<l,2I«C<l,JI-C<l,ll»G<3.2)i*GTYPE 

GG<3.2)"GC(2.3I 

GG< 3.3I-<G< l,l»»C»2,2*-G< l , 2 V««2 » ‘OTYPE 


OGAPPX  <  1,  l  )»GGAP(  1 ,  l  >»GG<  l,  l  •♦OCAlM  1  .2  )»CG<  2,1  I  ♦OGAi't  l  .3  IPCCI  3.  I  1 
DGAP«X  <l,2>"DGAP(2.ll"CC(l,l)»0CA812.2l*C0<2,l).QCA»(2.3l"CG(3.ll 
OGAPXXI 1 , 3)»CGAP<  3, 1 1  »CG<  l.l  lAOCAIAt  3,2  lACGI  2 .1  I  aCGAYI  3  .3  I*GG  I  3. 1  I 
OGArPX<2.1l»0GAr<l,l)»GG<1.2l,0CAf,<1.2)»GGI2.2),0GA''<l,3)»CG<3.2l 
OCAP*X<2.2)'OGAI"<2.1I*GG<1,2  1»0GA8  <2,2  1  "Gc  12,21*0&a-<2.3I*GG<>.2I 
DCA8"X(2,3)"0GA8I 3  , 1 1 "GG I 1,2 1 ♦OGAM) 3.2I*GG(2,2>,0GAp( 3, 3) «GC< 3.21 
OCA«*'X(3,l)"OCA»(l,ll»GCIl,31.0GAf*<l,2)»CGt2,3)"CGA'<l,3)»GC(j.3> 
OGAPPX(3,2l"OGA"<2,ll*GC<l«3)»OGAB<2.2)*GGI2.3l‘OGAP|2,3)"CC<j.3> 
OGAP»>X(3,31*OGAA<<3»ll*GC<1.3),CGAP<3,2IAGGI2,3)*OGA(<<3«31»CG<3.3l 
033S-0GAPPX  <3,31 
I  F  <  I  COLIN  T  .EO.  1ICALL  ALIXILIE) 

THIS  IS  A  CALL  TO  SL8RGUTI6E  YlSCuS 
If  <AVI5<  [LAYER  I  .NE.  0.1  CALL  DIPTAU<3.C) 

THIS  IS  A  CALL  TO  SLBSCUTINE  STRESS 
CALL  0 l " T  AU  <  2 ,C 1 

If 1PPHYS.NE .0  . ANO .  ICCLNT.EC. I  I  CALL  PHYSICI2 I 
?AR-THICKN«.S»»EIGMT< I G  A  u  5  S , NGAU  SL 1 
PAR  SC»PAR«SCRG 

HN<l,l|,MN<l,ll»pARSa*<TA<j<l.ll»GaASEH,l)»TAUIl,2l*GoASE<2.l) 

•  ♦  T  avj <  1 ,3  )"G9aSE  <3.1)1 

mN<1,2I*e,MI,2IaP»RSG*<Tau<1,1I»G8ASE<1.2I»Tau<1.2)»G8ASE<2.2) 

•  » I AU< 1 .3 1 *G?ASE I  3.2 1  I 

HN<l,3)»MM1.31»PA'<SG«<TAU<l,l)»GaASE<1.31»TAu<1.2l»G.3ASE<2.3) 

•  «■  T AU<  I , 3 I»G8aS£ (3,311 

fN12,i>-KN12.1).?ARSG«ITAO<2.ll«GaASE<l,l|,TAij<2.2l*Gb»SE<2.ll 

•  »T AU  <2 , 3  I • G2  A  5  t <3,111 

HN<2,2l*»i»<2,2IAPAKSi»<TAu<2.ll»G8AS£<l,2l,rau<2,2i»GBASE<2.2) 

•  »TAU<2 . 3  I *G8ASE (3.2  11 

MN(2.3l-f.|Y<7,3l,PARSO»(TAO<2.ll»GaA5E<1.3|,TAul7.2l*G8ASE(2.3l 

•  ,T  AU<  2, 3  I AGPASf  <  3  .  3  I  I 
If  (OIRCH)  GO  TO  2  165 
PRSCOl"PARSQ"02<  11 
PRSC02*PARS  C"0  2 ( 2 1 


2  E  T A  241 
2  E  T  A  2  a  2 
2 E T A  243 
2  E  T  A  2‘* 

2  E  T  A  245 
2  E  (  A  2  A  b 
2  E  T  A  2  A  7 
2  E  T  A  2  A  8 
2ETA  249 
2  E  T  A  2S0 
2  E  T  A  2  S  1 
2  E  T  A  2S2 
2  E  T A  2S3 
2ETA  25 A 
2ETA  255 
2 E  Ta  256 
2  E  T  A  257 
2 E T A  258 
2ETA  259 
2E  TA  260 
2 E TA  261 
2  E  T A  262 
2  E  T  A  263 
2  E  T A  26A 
2 E  T  A  265 
2ETA  266 
2ETA  267 
2  E  T A  268 
2  £  T  A  269 
2  E  T  A  270 
2  E  T  A  271 
2  E  T  A  272 
2  E  T  A  273 
26Ta  274 
2 E  T A  276 
2  E  T  A  276 
2  E  T  A  277 
2 E  T A  279 
2 E  T A  279 
2E  TA  290 
2  E  T  A  281 
2  E  T  A  282 
2 E  T A  293 
2 E  T A  284 
2  E  T  A  295 
2  E  T  A  2  86 
2  E  T  A  297 
2  E  T A  288 
2 E T A  299 
2ETA  290 
2  ?  T  A  791 
2 E  T  A  292 
2  E  T  A  293 
2  E  Ta  294 
2  F  T A  295 
2  E  T  A  296 
2  E  T  A  2 97 
2  E  T  A  299 
2ETA  299 
2  E  T  A  3  00 
2  E  T  A  301 
2  E  T  A  302 
2  f  T  A  303 
2ETA  304 
2  E  T  A  305 
2ETA  306 
2  F  T  A  307 
2  E  T  A  3  0  9 
2  F  T  A  309 
2  c  T  A  310 
2  E  T  A  311 
2  r  T  A  312 
2F  Ta  3  13 
2  ETA  314 
2  E  T  A  315 
/ETA  316 
2  E  T  A  317 
2  E  7A  3  18 
2  E  T  A  319 
2  E  T A  320 


PRSQ03-PARSC*02I3) 

SIRS SLll. II. 121-SIRE SL<1. 11.  12 »♦ 

«  PR  $00 1* 

•  ITAUU,1>AGBASE(1,1)»TAL11,2|AG8A5E(2.1).T4U11.3)AC8ASEI3,1)I» 

•  PRSC02* 

•  IT AU 12. 1)*G BA SEtl,l).T4Ut2,2>AC8ASEt2, II. TAUI2, 31 *G BASEI3.il). 

•  PRS003* 

•  ITAUt3,l>*GBASEtl,l|.T4LI3.21*GBASE!2.1).TAUt3.3>*GBAS£l3.1)> 
STRE3LI2.il. I2l-3rRESLI2.Il.  12  I* 

•  PRS001* 

«  ITAUll,l)AG8AS£il,21.IAUll,2)AG8ASEI2,2).TAUll,3)AG8ASE13,2)|. 

•  PRSCD2* 

•  (TAUl2,l|AG8AS£U.2).TAUl2.2l*C8ASEI2.2l.rAut2,3)«CBAS£(3.2>)« 

•  PRSCD3* 

•  ITAUI3,l)*G8ASEtl«2I.TAU(3«2)*GBASEI2,2l-TAUt3,3l*GSASE(3,21) 
SIRE  SL  13, (I.I2I-STRESL  (3. I1.I2L 

•  PR  SOD  1 • 

•  (TAUIl.l)«GSASE<l.3|.rAU(1.2)AGBASE<2,3).TAull.3|AC8ASEt3.3>). 

•  PRSC02* 

•  ITAUI2.1>aCBASE11,3I»TAUI2.2)aGBASEI2.3).TauI2.3)»GBASEI3,3I)» 

•  PRSG03* 

•  ITAUt3.1IAG8ASE<l,3I.T4Ut3.2l*G8ASEt2,3>.rAut3.3)*GaASEI3,3)l 
PRS022-PARSC*22 

CAP2QI I  1 , l l, 12 I-CAP2Q1 1 1 . I  1. 121 .PRSC72«<  TAull, 1) *GBAS£ 1 l. 1 ) 

•  *TAUIl,2)«GBASEI2.l>. Tal« 1. 31»G8 AS E ( 3. 1 1  I 

CAP20U2, 1 1  ,  12  1-CAP2QU2, 1 1.  12I.PRSC22*!  TAU12. 1 1  ‘GBASE  1 1 ,1) 

•  .TAUI2,2)*GBASE(2.1).TAUI2.3MGBASEI3,1I) 

CAP202T  l.ll«I2)-CAP202<l.ll,!2)»PRSC22*<TAU(l.l>AGaA5£ll,2) 

•  .TAU(l,2l*C8ASE(2.2).rAUU.3l*GBASEI3.2)) 
CAP2G2l2«ll«12l-CAP2Q2t2.Il«12)»PRSC22*ITAL12,l)*CdASEIl,2) 

•  .rAU(2.2)*G8ASE<2.2|.TAU(2.3IPGBASEI3.21) 
CAP2a3(l,Il.l2l-CAP2a3tl.ll.i2).PRSC22AtTAutl,l)*G8ASEIl,3> 

•  .TAL(1.2)AGBASE(2.3I.TAU(1,3I*GBASE(3,31I 
CAP2Q3l2«Il,l2)-CAP203l2,Il,I2)»PRSC22*tTAut2,l)*GaASEll,3) 

•  .TAU(2,2|AGBASE(2,  31.TAU(2.31*CBASE(3.3I  I 
2  165  CONTINUE 

IF  (QSfEAR)  GO  TO  2  OS  1 

CXIt.ll-GBASElt.llAYYUll.l 1 .GBASE I 1.21AYYU12.1I 

•  «CBASE ( 1 . 31 *YYU< 3. 11 

CX(l,2I.GBASE(2.l I  *YYU I  1 . 1 1 .CBASE l2.2)*TYUt2,l) 

•  «GBASEI2.3)ayyli 3.11 

CX ( 1, 31- GBASE ( 2 . ll«YYUIl,ll .CBASE l3.2IAYYUI2.il 

•  .GBASE13.31AYYLI3.il 

CX 12, 11 -CBASE II. 11  ay yu 11,2). CBASE I  1.21 *Y YU  I  2,2  I 

•  .CBASE 1 1,3 IAYYL1 3,2 1 

CX(2.21-GBASE(2.1|ayyu(1,2I.CBASEI2,2Iayyu12.2I 

•  »GPA5E 12 , 3 IAYYL 1 3,2  I 

CX(2,3l«CaA5E13,llAYYull,2|.GBASE13,2IAYYL(2,2l 
a  »CBA5E13,31AYYLI3,2! 

CXI3, 11- C8  ASEll.llASNII.il. I2I»  C8  ASE11,2)*S»<12.[1. 121 
a  *CBASE11,3IASNI 3,11,12) 

CXI3,2)-CBASE12,1|aSn(  1,11,121  .GBASE  I2.21ASM2.11.12) 

•  ‘GBASE  1  2 . 3  I  *SM  3,  l  1 ,  12  I 

CX13. 31 -C8ASEI3, DASMI.il. 12).  CBASE  (3. 21ASN12. II. 12) 

•  »CBA5E  13.31  *SM  3,11,12) 

PARSC2 -PARSCA2 

H- 1 1 , l ) -HA ( 1 . 1 ) .P ARSL2 • 1 T AL 1 1 . I )ACX 1 l , l ) »  t AUt l .2  I  »C < 11  ,2  ) 

•  »TAU( 1 .3 )*CX1 t . 3 ) I 

MM  1,2  l-HM  1,2).PAR5  02  A1TAUII,1)ACX12,1).TAUI1»2)ACX12,2I 

•  » TAUt 1 . 3  I »CX 12 . 3 ) I 

H1M2,1).|-M2,1).PARS02A1TAU(2.1  I*CXU.1).TAU(2,2)»CXU,2) 

•  »  T  All  12.3)*CX  11,3)1 

HN 1  2,2  >*HM  2.2  >»P  ARSQ2  •!  TAU 1 2, 1  )*CX  ( 2,  l  )*TAul  2 .21*0X12,2  I 

•  .TAUt2,3)*CX12,3) ) 

STRESP1 1»  11. 12 1-STRESP1 1. 1 1. I2).P*RS0*«  TAU13.1 )»CX1 1,1 1 

•  »TAU(3,2)*CXI 1 ,2 |.TAU13,3)ACXt 1. 3) 1 

STRESPI2. ll.I2)-STRESP»2. II. 12 1 .PARSQA ( T AU 1 3.1 )*CX!2,l 1 

•  .TAU13,2)ACX12,2|.TAU13.3|ACX(2,3II 
STRES0«1,H,I2)-STRES01 1.1 1,  I2I.PARS02A1  TAUI  l.l  1ACX13, 1) 

•  ‘TAUI1 ,2) »C XI 3.2 1 »TAUI 1 , 3  I *C XI  3. 3) ) 

STRES0I2, ll.I21-STRESQ12.il. 12). P»BSa2A(TAL(2,l) »C X13.1) 

•  .TAUI2,2I*CX13,2).TAU12,3)«CX13.3)) 

2059  CONTINUE 

0CB33-CCP33.033S*PAR 

SUPC-SUPC.CC13.31APAR 

IE  1  I  COUNT  ,EQ.  1)  CALL  ENERGY! 3, PAR ) 

5GAP33-SGAN33.CGAB33APAR 
If  1 IN0RPL.NE.2 I  GO  TO  169 
TAUSUN-CTNI 1. 1  I .DTN12.2 I.0TP1 3.3) 


2 E  T  A  321 
2E  Ta  322 
2 E  TA  323 
2 E  T  A  32s 
2ET  A  325 
2ETA  326 
2  E  T  A  327 
2  E  T a  328 
2  E  T A  329 
2  E  T  A  330 
2  E  T  4  3  3  1 
2  6  Ta  332 
2  E  T A  333 
2 E  Ta  33s 
2  F  T  a  335 
2 E  T A  336 
2  E  T  A  337 
2  E  TA  338 
2  fc  T  A  339 
2ETA  3  SO 
2  E  TA  3AI 
2  E  TA  3S2 
2ETa  3S  3 
2 E  Ta  3SS 
2 E  T A  3s5 
2  E  T  A  3S6 
2ETA  3  S  7 
2  E  T 4  3S8 
2ET4  349 
2E  Ta  350 
2  E  T  A  351 
2 E  T A  352 
2  E  T  A  353 
2 E  T A  35s 
2  E  T A  355 
2  E  Ta  356 
2 E  T A  357 
2ETa  358 
2  E  T  A  359 
2E  Ta  360 
2  E  T  A  361 
2  E  T  A  362 
2 E  T A  363 
2E  Ta  36s 
2  E  T A  365 
2 E  T A  366 
2 E  T A  367 
2  E  7  A  368 
2  E  Ta  369 
2  E  T  A  370 
2  E  T  A  371 
2 E  T A  372 
2  E  T a  373 
2  E  T  A  3  7s 
2 E  T A  375 
2ETA  376 
2  E  T  A  377 
2ETA  378 
2 E  T A  379 
2 E T A  380 
2  E  T  A  381 
2  E  T  A  382 
2ETA  3B3 
2  E  TA  38s 
2ETA  385 
2  E  T A  386 
2ETA  387 
2 E  T  A  388 
2ETA  389 
2 E  T  A  390 
2  E  T  A  391 
2 E  T A  392 
2  E T A  393 
2  E  T A  394 
2  E  T A  395 
2  E  T  A  396 
2  E  T  A  397 
2  E  T  A  398 
2  E  T  A  399 
2  E  T  A  sOO 


DGAf*A3-QGA*A3*  «  (  (  1.-2.  *HNU  WEE  )  • TAUSU* 

♦  -OCAHHXU  .U-CGANNXI2 ,2  >*O.*ALPHA*0TE“P  WPAR 
169  CONTINUE 

0GM3-0GH33/THIC 
AV6G33-SunG/THlC 
AGA»»33-SGAf‘33/TH  IC 
I F  < ICOUNT  .go.  I ) GOTO  180 
0G0C-2.0*(CGM33-AGArt33)/AV6G33 
DEL  BAR C  I  1  ,  12  )— GB  rs^SQRTt  GBTN*»2-OGOG1 
180  CONTINUE 

SURFACE  STRAINS  FOR  PLOTTING 
l F ( I  COUNT  .60.  QIGGTO  190 

IF  (  Il.EQ.ISl  . ANC •  I2.fcG.IS2)  ? G A n 3 3 -PC A *  3 3 ♦ AGA n 3 3 
I F  (  I  COUNT  .60.  1)  CALL  AOXILCIGI 

EPStlt  I  1*12  J-EPSLl  t  1 1*  I  2  >  *0.  SMCAt  1  ,  1 )  ♦THIC*D8U  ,  L)  ) 
6PSL2m,I2)-EPSL2ClUI2)*0.5#I0A(2,2)*TNlC*0B(2.2J) 
GANNALl  li«l2)*GAPHALl 1 1,12)*0.6*{0A( i,2l*THlC*CB( 1,2 )  I 
EPSU1C I  1  *12  >*EP$U1  (I  l  •  12  MO.  5*  (  CA4  1 ,  U-TH  IC  *03  l  l  ,  1)  ) 
6PSU2(Il.I2)*6PSU2(Il»I2)*0.S*IDA(2,2)-THlC»0b(2,2)> 
GANNAUC  1 1  ,  12  )-CAfNAU(  li,I2)»0.5*(DA<l,2)-THlC*CB(1.2)> 
190  CONTINUE 

IF (  INGRPL.NE  .2  I  GO  TO  200 
AVERAGING  AT  EACH  POINT 
0GANA3-0CANA3/ TH  ICK2 
SETTING  Y2  ANO  C2 
IFtCIROMGO  TO  460 

D2( 1,1 1 « 121-02 ( l, I l « 12 ) «0GAPA3*SN( l, l l, 12) 

02(2,11, 12)*02(2»11,I2WDGA*A3*SN<2»11, 12) 
02<3,ll,t2)-D2(3,il,I2)*0GANA3*SN(3,l;,[2) 

*50  Y2(ltIl«I2>l'Y2(ltll«12)+CGAPA3*$N(  1 ,  1 1 «  I  2  ) 
T2(2«I1.I2)"T2<2*I1«I2)»0GAPA3*SN(2*I1*I2) 
Y2<3.I1,I2)-V2(3*U*I2M0GAPA3aSN(3*U,I2) 

200  HNl  (  1.11,12)  -HM  1  ,  1) 
hn 1I2.II.I2  >»«*(  1,2) 

HN2  (1,11,12)  *MP  (2.1) 

HH2(2,I1,I2)-HP(2,2) 

ASSUMING  CRTHOCCnau  C0CRCINAT6S,  TO  «6CUCE  ERRCP 
IF  (K l.AE . 7. AND. *2. N6  .  7  )  RETURN 


2  6  T A  *01 
JET  A  402 
2  E  T  A  A  0  3 
2 E  T  A  40* 
2  E  T  A  406 
2  E  T  A  *06 
2 F  T A  40? 
2  E  T A  408 
2  E  T  A  409 
2ETA  *10 
2  E  f  A  4  L  1 
Z  6  T  A  ^12 
2  6  T  A  *13 
2ETA  41*, 
2  E  T  A  4  l  S 
2  E  T  A  4lt> 
2  E  T  A  4  l  7 
2E  TA  <*18 
2  E  T  A  4  19 
2  E  T  A  420 
2  E  T  A  421 
2  E  T  A  422 
2f  TA  *2  3 
2FTA  424 
2  F  T  a  *,26 
2  E  T  A  426 
2  E  Ta  4 2  1 
2  E  T  A  42a 
2  E  T  A  47=) 
2  E  T  A  430 
2  E  T  A  4U 
2  F  T  A  4J2 
2 ETa  433 
2  F  T  A  4). 
2  E  T  A  436 
2E  Ta  4  36 
2  E  T  A  *37 
2E  TA  4)8 
2  F  T  A  4)9 
2  E  T  A  4  4  0 


HHl  (2,  I  l,  12  )-0. 
HN2(l,  11,121-0. 
RE  TURN 
END 


2  6  T  A  441 
2  E  T  A  442 
2  E  T  A  443 

2  6  T A  444 
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